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FOREWORD 


The work for this report was per- 
formed under contract number NAS 9-14203. 
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assistance of the Contract Technical Monitors, 
Messrs. R. T. Savely and E. R. Schiesser of 
the Math Physics Branch of the Mission Planning 
and Analysis Division of NASA/ JSC. Assistance 
and cooperation were also received from Messrs. 
T. J. Blucker and R. £. Eckelkamp. 

This is the second of two reports 
submitted under the contract and represents 
work performed between September 1975 and 
July 1976. The results of work completed 
prior to September 1975 were presented in 
Ref. 1. 
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ABSTRACT 


A detailed analysis of Space Shuttle 
navigation for each of the major mission phases 
is presented. A covariance analysis program for 
prelaunch IMU calibration and alignment for the 
Orbital Flight Tests (OFT) is described and a 
partial error budget is presented. The ascent, 
orbital operations and deorbit maneuver study 
considered GPS-aided inertial navigation in the 
Phase III GPS (1984+) time frame. The entry and 
landing study evaluated navigation perfoirmance 
for the OFT baseline system. Detailed error 
budgets and sensitivity analyses are provided 
for both the ascent and entry studies. 
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1 . 


INTRODUCTION 


The navigation system for the Space Shuttle will 
employ a variety of navigation aids in providing accurate 
navigation for the entire mission scenario from launch to 
landing. Previous studies by TASC have focussed on the de- 
velopment of error models for these navigation aids, the 
formulation of navigation software algorithms based upon the 
error models, and the generation of detailed performance 
projections for candidate navigation system implementations. 
As indicated in Fig. 1-1, these efforts have encompassed 
Space Shuttle navigation during each of the major mission 
phases. TASC has continued these activ'ities in the current 
study with emphasis on evaluation of potential changes in the 
navigation system, including both software changes and new 
navigation aids, which could significantly improve Space 
Shuttle navigation performance. 

The current study has considered three facets of 
navigation system performance: 

• Prelaunch IMU calibration and alignment 

• GPS -aided* navigation for the ascent, 
orbit and deorbit mission phases 

• Entry, preland and landing navigation 

The cal/align and entry studies addressed navigation sys- 
tem performance for the current Orbital Flight Test (OFT) 
baseline system. The GPS study evaluated the improvement 

NAVSTAR Global Positioning System 
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Figure 1-1 Summary of TASC Studies Relative to 
Space Shuttle Navigation (1973-1976) 


in navigation accuracy which could be attained by using GPS 
as the primary navigation aid for the Space Shuttle. Each of 
these studies is discussed in more detail in the following 
section. 


1.1 OVERVIEW 

The guidance, navigation and control function for 
the Space Shuttle will be provided by the Operational Navi- 
gation Program (ONP). Those elements of the ONP which directly 
affect nav'igation performance are the: 

• IMU Subsystem Operation Program (SOP) 

• Navigation software 

• Redundancy management software 
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The I MU SOP performs the basic ground function of calibrat- y 
ing and aligning the IMUs. During flight, it interrogated 
the IMUs and compensates the measured inertial velociti^ 
for the calibrated instrument errors. The inputs to tl^e navi- 
gat ion software are the compensated inertial velocities as 
well as measurements obtained from the external navigation 
aids. The navigation software processes these inputs to 
generate estimates of vehicle position and velocity. The 
redundancy management software -^--rforms basic consistency 
checks on data generated at the subsystem level, e.g., IMU 
outputs, sensor outputs, position estimates, etc. The studies 
presented in this report deal only with the IMU SOP and the 
navigation software. The scope of the studies is outlined 
in the following subsections. 


1.1.1 Calibration and Alignment Analysis 

In January 1975, TASC initiated a verification and 
validation (V&V) program for the I'lU soft -ctre algorithm. The 
program was formulated in three phases (see Fig. 1.1-1) — 
analytical, functional and performance. As each phase was com- 
pleted, a series of recommended algorithm modifications were 
delivered to NASA and, following approval by NASA, incor- 
porated into the IMU software baseline. A summary of the 
first two phases, and part of the third phase, was presented 
in Ref. 1. A summary of subsequent progress on the perfor- 
mance phase is included in this report. 

The performance phase of the V&V program applies 
three different techniques: 


1 ^ 
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Figure 1.1-1 TASC Approach to Calibration and Alignment 

Software Verification and Validation 


• Comparison of KT-73 calibration and 

alignment test data with the Shuttle 
performance specifications 

• Monte carlo analysis 

• Covariance analysis 

Since the KT-73 IBIU is generically related to the Space Shuttle 
nro, P.rd since the ca?.ibration programs are also similar, the 
KT-73 Data Bank (Ref. 19) developed and maintained by TASC for 
the U.S. Air Force is a valuable tool for ev’^aluating the abil- 
ity of the IMU software to calibrate and align the Space Shuttle 
IMUs In an operational environment . The Data Bank also pro- 
vides a useful data source for updating the IMU error models 
and in many ways will complement the IMU test data currently 
being generated by NASA. 

The monte carlo and covariance simulations provide 
complementary vantage points for computer-based analysis of 
the IMU software performance. The monte carlo simulation is 
a high fidelity simuJation developed to ensure that the IMU 
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software algorithm can control the IHUs and execute the cal/ 
align function. The simulation is extremely useful for 
generating coarse overall performance projections and for 
evaluating the effect of system nonlinearities on the soft- 
ware performance. Unfortunately, detailed sensitivity ana- 
lyses or accurate overall performance projections can be 
prohibitively expensive with monte carlo simulations. 

The covariance analysis simulation is based on a 
linearized IMU error model. Subject to the validity of this 
linearization, the covariance analysis simulation will yield 
a more accurate estimate of overall performance than is 
attainable with the monte carlo simulation and will also pro- 
vide a detailed breakdown of the contribution of individual 
error sources to the calibration and alignment performance. 
The degree to which the nonlinearities alter the covariance 
analysis results can be inferred by comparison with the 
monte carlo results. Together, the two approaches will per- 
mit an accurate assessment of calibration and alignment per- 
formance and will provide a sound basis for reconunending and 
evaluating software modifications which would improve overall 
performance. 

The results of the KT-73 comparison and the monte 
carlo analysis were presented in Ref. 1. Development of the 
covariance analysis program was completed as part of the 
current effort. However, because of recent changes in the 
IMU software algorithm (Ref. 18), a complete error budget has 
not been generated. The error budget will be included in 
a subsequent report. 

1.1.2 GPS- Aided Navigation 

GPS is a satellite navigation system designed to 
provide accurate position and velocity fixes on a global 
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scale. The GPS network will become available on a limited 
basis by the end of 1977 (Phase I). The complete network 
will be operational in 1984 (Phase III). Because of its 
accuracy and its ability to provide navigation coverage 
during virtually all phases of a Space Shuttle mission,* 

NASA is currently evaluating the possibility of including 
GPS as a navigation aid. The objective of the GPS-aided 
navigation study was to provide NASA with an assessment of 
potential navigation accuracy in the operational GPS (post- 
1984) time frame. The study was restricted to those mission 
phases (ascent, orbit, deorbit) where no navigation aids are 
included in the current baseline, although GPS could also re- 
place the current entry and preland navigation aids. 

GPS is designed to satisfy the navigation require- 
ments for both military and appropriately equipped civilian 
users. Seven classes of user requirements have been defined 
by the GPS Joint Program Office. Each class has its own 
unique requirements in terms of the type of user equipment 
needed and the navigation data required. The seven user 
classes are: 


Class 

A - 

High accuracy users such as 
strategic and reconnaisance 
aircraft 

Class 

B - 

Highly dynamic users such as 
fighter aircraft which require 
high accuracy updates 

Class 

C - 

Low-cost , low-accuracy users 
such as transport aircraft 

Class 

D - 

Land and sea-based surface 
vehicles 


♦No coverage will be available during radio-blackout. Also, 
the navigation accuracy may be marginal for landing. 


!< 
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• Class E — Manpacks 

• Class F — Submarines 

• Class M — Missiles 

Two receiver types, the X and Y sets, are designed for users 

with high accuracy requirement . A third receiver type, the 
Z set, is intended for users with lower accuracy requirements. 

The intent of the study was to evaluate GPS-aided 
navigation performance with the best available hardware. As 
a consequence, the Space Shuttle was assumed to be equipped 
as a Class B user. This included a multichannel X-receiver 
and processing of the high accuracy (P-code) satellite trans- 
missions. However, it did not include a cesium time standard 
onboard the Space Shuttle. If the Shuttle Master Timing 
Unit (MTU) were replaced by a cesium time standard, navigation 
accuracy would be somewhat better than predicted in this 
study . 


The GPS study had three basic objectives: 

• Develop a navigation filter for GPS- 
aided navigation 

• Formulate candidate GPS/Shuttle mea- 
surement scenarios 

• Assess GPS-aided navigation 

The emphasis in the navigation filter development was to re- 
tain the basic structure of the entry navigation filter 
(Ref. 45), with state definitions and parameter values suit- 
ably modified. The filter design as well as the performance 
analysis relied heavily upon GPS error models and GPS simu- 
lation programs developed by TASC in previous efforts (Refs. 
1 and 5) . 
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1.1.3 Entry, Preland and Landing Navigation 

The Space Shuttle will employ a number of different 
navigation aids during the entry, preland and landing mission 
phases. In addition to the radio navigation aids indicated 
in Fig. 1.1-2, altitude measurements from a drag pseudo- 
measurement, a baro altimeter and a radar altimeter will be 
utilized at various points during the mission. 


R-U200 



MLS DME and 

azimuth 


Figure 1.1-2 Baseline Orbiter Radio Navigation Aids 

An accurate navigation performance analysis requires 
error models for the IMU and the navigation aids which ade- 
quately reflect instrument performance. Thus a major part 
of the entry, preland and landing study was the improve- 
ment of previously developed error models. Particular atten- 
tion was given to the drag pseudo-measurement and radar 
altimeter error models, although recent performance data 
from the MLS manufacturer was also incorporated into the MLS 
error models. 
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1.2 ORGANIZATION OF THE REPORT 

Chapter 2 siunmarizes the results of the GPS-aided 
navigation study. A candidate GPS measurement schedule is 
presented along with a description of a navigation filter 
suitable for processing GPS measurements during ascent , or- 
bital operations and the deorbit maneuver. Navigation per- 
formance projections are presented for each of the mission 
phases for which the filter was designed. Performance sen- 
sitivity to the dominant error sources is also evaluated. 

A description of the covariance analysis program 
to be used to analyze IMU calibration and alignment perfor- 
mance is provided in Chapter 3 . Although a complete error 
budget is not given, calibration uncertainties Jue to certain 
key error sources are estimated. These estimates are com- 
pared with previous sensitivity analyses based upon TASC’s 
monte carlo simulation. 

The results of the entry, preland and landing study 
for the OFT-1 mission are contained in Chapter 4 . Detailed 
error budgets are provided for the end of radio blackout , MLS 
acquisition and touchdown. The projected navigation perfor- 
mance is compared with results from previous studies and 
recommendations for improving performance are made. 

A brief summary of results and conclusions is pro- 
vided in Chapter 5 . 
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2 . 


GPS-AIDED NAVIGATION 


The navigation error analysis for the Space Shuttle 
ascent phase provides an evaluation of navigation accuracy 
at orbit insertion if the baseline nav'igation system were 
modified to include the NAVSTAR Global Positioning System 
(GPS). The assumption in this study is that P-code (high 
accuracy) measurements from a GPS X-receiver are to be 
used periodically during ascent to update the navigation 
filter state estimates. The results are then extended to 
permit a preliminary evaluation of GPS-aided navigation per- 
formance during both the orbital and deorbit maneuver mis- 
sion phases. 

The methodology used for evaluation of the GPS-aided 
Space Shuttle navigation performance is summarized in Fig. 

2-1 . The first step in the process is the GPS segment 
simulation. This step is required to: 

• Develop a time history of satellite 
orbital positions and velocities 

• Generate the satellite position, 
velocity and clock error covariances 

The GPS segment is performed using the TASC Navigation Error 
Sensitivity Analysis (NESA) program. 

The second st-ep is represented by the FLTCOV module 
in Fig. 2-1. The satellite ephemeris data, along with 
the Shuttle trajectory and measurement schedule, are used 
to generate a recursive solution of the navigation filter 
covariance propagation and update equations. This step is 


2-1 
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Figure 2-1 Methodology for GPS/Shuttle Performance Projections 
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executed only once, with the filter gains being recorded on 
magnetic tape. 

The third step (the SYSTEM module) is based on the 
system truth model and consists of the recursive solution of 
the linear system covariance equations. These equations are 
solved repeatedly to produce an error budget , the same gain 
file being used each time. The satellite error covariance 
generated in the NESA step is used to generate one line of 
the error budget and models for other error sources (IMU 
errors, BJTU errors, etc.) are used to generate subsequent 
lines. An overall system performance projection is then cal- 
culated based on the detailed error-source-by-error-source 
breakdown . 

The results of the ascent phase covariance analysis 
are summarized in this chapter. Section 2.1 describes the GPS 
configuration used in the NESA step and Section 2.2 discusses 
the relationship of the Orbiter trajectory and measurement 
schedule to the GPS satellite positions. The Shuttle navi- 
gation filter covariance generated in the filter gain step 
is discussed in Section 2.3. Section 2.4 summarizes the truth 
model used in the SYSTEM module. The performance assessment 
generated by that module is presented in Section 2.5 and a 
sensitivity analysis is presented in Section 2.6. A first-cut 
estimate of GPS-aided navigation performance during the orbital 
and deorbit maneuver phases is presented in Section 2.7. 


2.1 GPS CONFIGURATION 

GPS is a navigation satellite system under develop- 
ment by a multi-service Joint Program Office within the U.S. 
Air Force Space and Missile Systems Organization. The GPS 
satellites will each broadcast signals containing coded 
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Information giving the time the signal was transmitted and the 
satellite position and velocity at that time. A user with an 
L-band antenna, a GPS receiver, and a clock synchronized with 
the satellite clocks will be able to compute the range to a 
satellite by measuring the time delay between signal trans- 
mission and reception, and he will be able to compute the 
range-rate to a satellite by measuring the signal’s doppler 
shift (see Fig. 2.1-1). Since clock synchronization errors 
will enter into these measurements, they are termed pseudo - 
range and pseudo-range-rate measurements. Pseudo-range mea- 
surements to four different satellites will allow the user to 
solve for his three position coordinates and his clock phase 
synchronization error. Similarly, pseudo-range-rate mea- 
surements to four different satellites will allow the user to 
solve for his three velocity coordinates and his clock freq- 
uency synchronization error. Measurements to more than four 
satellites can reduce the uncertainty of the results. 

••lOMU 



Figure 2.1-1 Elements of the NAVSTAR Global 

Positioning Svstem (GPS) 


*Provided the four satellites do not all lie on a cone with apex 
at the user location. See Appendix E of Ref. 4 for details. 
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Deployment of GPS is to be conducted in three phases 
(Table 2.1-1). The final configuration, Phase III, will con- 
sist of 24 satellites with cesium standard clocks and will 
provide full global navigation coverage for both ground- 
based and orbital users. Each satellite will have an orbital 
period of one-half a sidereal day (11 hours and 58 minutes). 
Thus, the satellites will nominally repeat the same ground 
track every two orbits, or once a day. Table 2.1-2 present.*" 
the orbital parameters of the Phase III satellites. A sub- 
set of these satellites that give good coverage over the 
Western Test Range (WTR) at the time of the Space Shuttle 
launch was used in the ascent phase navigation analysis. 

The eight satellites included in the subset have been 
denoted by asterisks in the tables. 

TABLE 2.1-1 

GPS DEVELOPMENT SCHEDULE 


n42«09 



OPERATIONAL 

DATE 

NTCTER OF 
SATELLITES 

% TI'JE FOR which 4 * SATELLITES 
.ARE VISIBLE .AT ORBITER ALTITUDE 

rKASE I 

LATE 1977 

6 

25 

PHASE II 

1981’" 

9* 

90 

PHASE III 1 

1984 

24 

100 


»A PROPOSAL TO ACCELERATE PHASE II DE\-ELC?VEN*T ANT DEPLOY 12 SATELLITES 
IS U:-:DER CONSIDERATION. 


The satellites will be tracked by four ground track- 
ing stations located in United States territory. Pseudo- 
range measurements taken at each tracking station will be 
forwarded to the master tracking station at Vandenberg AFB 
for real-time satellite ephemeris estimation and clock 
calibration. The locations of the tracking stations are given 
in Table 2.1-3. 
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TABLE 2.1-2 

PHASE III SATELLITE INITIAL CONDITIONS 


Satellite 

No. 

Semimajor 

Axis 

Eccentricity 

Inclination 

Right 
Ascension 
of the 
Ascending 
Node 

Argument 

of 

Perigee 

Mean 

Anomaly 


a 


1 

6 

b. 

u 


(earth radii) 

e 

(dog) 

(deg ) 

(deg ) 

(deg ) 

1 

4.16449 

0 

63 

0 

0 

15 

2 

4.16449 

0 

63 

0 

0 

60 

3 

4.16449 

0 

63 

0 

0 

105 

4* 

4.16149 

0 

63 

0 

0 

150 

6* 

4.16449 

0 

63 

0 

0 

195 

6 

4.16449 

0 

63 

0 

0 

240 

7 

4.16449 

0 

63 

0 

0 

285 

8 

4.16449 

0 

63 

0 

0 

330 

9 

4.16449 

0 

63 

120 

0 

345 

10 

4.16449 

0 

63 

120 

0 

30 

11* 

4.16449 

0 

63 

120 

0 

75 

12 

4.16449 

0 

63 

120 

0 

120 

• 

13 

4.16449 

0 

63 

120 

0 

165 

14 

4.16449 

0 

63 

120 j 

0 

210 

15 

4.16449 

0 

63 

120 

0 1 

255 

16 

4.16149 

0 

63 

120 

0 

300 

17* 

4.16449 

0 

63 

240 

0 

0 

is’ 

4.16449 

0 

63 

240 

0 

45 

19 

4.16449 

0 

63 

240 

0 

90 

20 

4.16449 

0 

63 

240 

0 

135 

21 

4.16449 

0 

63 

240 

0 

ISO 

22 

4.16449 

0 

63 

240 

0 

225 

23 

4.16449 

0 

63 

240 

0 

270 

24» 

4.16449 

0 

63 

240 

0 

315 


TABLE 2.1-3 

GROUND TRACKING STATION LOCATIONS 


Station 

No. 

Name 

Latitude 

(deg) 

Longitude 
(deg ) 

1 

Vandenberg AFB 

34.7 

-120.5 

2 

Elmendorf AFB 

61.3 

-149.5 

3 

Hawaii 

20.0 

-158.0 

4 

Guam 

13.0 

145.0 
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GPS requires periodic uploading of the satellite ephe- 
meris and clock calibration parameters in order to maintain 
navigation accuracy. Uploading of a satellite with its 
currently estimated ephemeris will occur once a day when 
that satellite makes its closest approach to Vandenberg AFB. 

The time of the uploading to each satellite is given in 
Table 2.1-4. These times are relative to the initial time 
assumed in Table 2.1-2. 


TABLE 2.1-4 

« 

SATELLITE UPLOADING TIMES 


Satellite 

Ko. 

Time of 
Uploading 
min 

Satellite 
No. . 

Time of 
Uploadirg 
min 

1 

870 

13* 

1260 

2 

750 

14 

1140 

3 

615 

15 

1020 

4* 

495 

16 

915 

5* 

405 

17* 

90 

6 

1230 

18* 

1410 

7 

1140 

19 

825 

8 

1005 

20 

720 

9 

315 

21 

615 

10 

210 

22 

480 

11* 

105 

23 

360 

12* 

1395 

24* 

225 


The net satellite contribution to pseudo-range and 
range-rate measurement error is not simply related to the 
satellite-transmitted ephemeris errors, because satellite 
clock errors are also involved, and the clock errors and 
ephemeris errors are highly correlated. In the case of 
pseudo-range measurements, the net error to the user is 
typically much smaller than the satellite position error 
itself. This happens because the combination of clock 
errors and position errors important to the user is nearly 
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the same combination measured by the ground tracking stations. 
While the ground network lacks good information to correctly 
allocate satellite position errors and clock errors indepen- 
dently, it has good information on the combination of these 
errors important to the user. 

An example of the satellite contribution to user 
measurement errors and the importance of the ephemeris 
upload in maintaining measurement accuracy is provided in 
Fig. 2.1-2. Prior to the upload, the user measurement 
errors due to satellite errors alone are on the order of 
5 ft and 0.003 fps. Immediately following the upload the 
errors are reduced to less than 2 ft and 0.001 fps. These 
errors define the fundamental limit of GPS measurement 
accuracy if all user-related errors (receiver errors, propa- 
gation errors, clock errors) could be eliminated. 



TIME (hr) 

Figure 2.1-2 Typical User Measurement Error Statistics 

Due to Satellite Ephemeris and Clock Errors 
(Phase III GPS) 
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The user-furnished equipment for GPS measur^ent 
processing consists of an L-band antenna, a GPS receiver and 
a clock. In addition to GPS satellite eph^eris errors, GPS 
measurenent accuracy is a function of the particular receiver 
(and clock) selected by the user. Several classes of GPS re- 
ceiver have been defined. The principal difference from a 
performance vantage point is whether the receiver can process 
the P-code (high accuracy) or C/A-code (medium accuracy) GPS 
signals — the multichannel X-receiver can process both signal 
types while the single channel Z-receiver can process only 
the C/A-code. The receiver single-measurement performance 
specifications for a 0.1 sec integration time are summarized 
in Table 2.1-5. Greater accuracy can be achieved by averaging 
successive measurements . 

The design specifications for the GPS X set require 
a user clock of comparable accuracy to the Shuttle Master 
Timing Unit (MTU) — short-term accuracy of 1 part in 10^® 
and an aging rate of 1 part in 10® per day. The user clock 
error model presented in Section 2.4 is based upon MTU test 
data. 

TABLE 2.1-5 


GPS RECEI\T:R ?1EASUREMENT ACCURACY" SPECIFICATION (REF. 8) 


CODE TYPE 

Icj RANGE ERROR 
(ft) 

lo RANGE-RATE ERROR 
(fps) 

Military (P) Code 

5 

0.2 

Cl ear /Acquis it ion 

50 

0.2 

(C/A) Code 
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2.2 TRAJECTORY AND MEASUREMENT SCHEDULE 

The GPS Navigation Error Sensitivity Analysis (NESA) 
computer program performs a statistical covariance analysis 
of the satellite system interacting with the ground tracking 
stations. Beginning with some initial condition statistics, 
to which this process is quite insensitive, the NESA program 
was used to simulate two days of GPS operation, achieving a 
close approximation to steady state conditions. The simula- 
tion correctly accounts for the fact that each satellite has 
its own "idea" of its ephemeris, clock synchronization, etc., 
based on the daily uploading from the master tracking station 
at Vandenberg AFB. Because each satellite's ground track and 
the time of its data uploading repeat every sidereal day, so 
do the statistics of the associated error processes, once the 
transient effects of the initial statistical conditions have 
faded. The steady state situation referred to above thus 
exhibits time variations within the day. 

After the two day settling period, a 90 min seg- 
ment of the statistical simulation was run with output at 
one minute intervals. This 90 min segment of GPS simulation 
was the input to the filter covariance program and to the 
truth model program for the GPS-aided Space Shuttle ascent 
navigation study. In order to do this, it was necessary 
to make a time connection between the Shuttle ascent trajec- 
tory and the NESA output. 

The Shuttle trajectory used for this study was 
reference mission 3B with a launch from Vandenberg AFB. The 
ground track for the ascent trajectory, along with the 
location of the eight GPS satellites used in the simulation 
(assuming a launch at t=0) is shown in Fig. 2.2-1. Its 
key events are suiranarized in Table 2.2-1 From liftoff to 


Jli 
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Figure 2.2-1 Baseline GPS Satellite Coverage 
(Space Shuttle Ascent Phase) 

TABLE 2.2-1 

TRAJECTORY PARAMETERS AT KEY EVENTS DURING ASCENT 


T-ooes 


Event 

GET 

(sec) 

Inertial 

Velocity 

(fps) 

Relative 
Flight-Path 
Angle (deg) 

Altitude 

(it) 

Geodetic 
Lat itude 
(deg) 

Longitude 

(deg) 

Lift-off 

0.0 

1,251.3 

0.00 

-..25 

35.000 

239.500 

Initiate 

Pitchover 

5.0 

1,257.1 

90.00 

297.0 

35.000 

239.500 

Begin Gravity 
Tarn Phase 

20.0 

1,350.5 

66.18 

5,159.0 

34 . 999 

239.500 

Staging 

116.4 

4,163.2 

65.61 

169.910.0 

34.850 

239.432 

Begin Con- 
stant 3g 
Accelera- 
tion Phase 

363.4 

j 

15,454.4 

-1.40 

584,412.0 

29.873 

237.488 

Begin Coast 
Phase 

486.5 

25,250.2 

0.06 

562,667.0 

24.440 

235,467 

Begin OMS*to 
Oi^it Phase 

5X6.5 

25,250.5 

1 

0.01 

1 

561,474.0 

22.482 < 

234.762 

Orbit 

Insertion 

771.0 

25,597.1 

0.12 

551.798 0 

5,673 1 

1 

1 229.225 


*Orbit ManeuveriDg System 
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o:?bit insertion, the ascent takes only 771.0 sec, or less 
than 13 min; and the measurement taking portion of the 
trajectory, scheduled between 630.0 and 768.0 sec after 
liitoff, takes less than 138 sec, or 2.3 min. This period 
covers the second half of the final OMS insertion maneuver 
and should permit sufficient time for the aided navigation 
filter to improve the state estimates and for the guidance 
system to refine the insertion parameters. Thus, no apparent 
advantage would be obtained by starting the measurement 
processing earlier. 

The measurement schedule calls for sequential process- 
ing of GPS measurements at the rate of one each 6 sec* starting 
with Satellite No. 4. Thus, over the 138 sec measurement inter- 
val, three measurements from each of the eight visible satellites 
are used to update the navigation filter state estimates. The 
measurement interval is shown relative to the ascent accelera- 
tion profile in Fig. 2.2-2. It is worth noting that the mea- 
surements were scheduled during a low acceleration portion of 
the mission. This was done to ensure that acceleration-sensi- 
tive errors in the Space Shuttle Master Timing Unit (MTU) would 
not significantly degrade the GPS measurement accuracy; how- 
ever, the results in Section 2.5 suggest that the navigation 
accuracy is reasonably independent of the acceleration profile 
during the measurement sequence. 


In order to connect the Orbiter trajectory to the NESA 
output , it is only necessary that the measurement interval of 


*GPS receivers are capable of generating measurements at 
a much higher data rate. The slower measurement rate was 
selected to reduce computer requirements. 


3 » 
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Figure 2.2-2 Ascent Tra.iectory Acceleration Profile 


2.3 min fall somewhere within the 90 min of NESA output. The 
Phase III GPS constellation provides good satellite cover- 
age on a continuous basis. Because the NESA run was made with 
a subset of eight satellites, however, the results are repre- 
sentative of Phase III GPS coverage in the WTR for only a short 
period each day. In order to determine where the measurement 
interval should fall within the NESA output for simulation pur- 
poses, it was necessary to investigate satellite visibility from 
the Orbiter. Figure 2.2-3 indicates satellite visibility over 
the output period from four points. The four points are 

V) on the ground at Vandenberg AFB; 

B) on the Orbiter at 630.0 sec after 

liftoff, the beginning of the measure- 
ment interval ; 

M) on the Orbiter at 702.0 sec after 

liftoff, about midway into the measure- 
ment interval ; 

E) on the Orbiter at 768.0 sec after 

liftoff, the end of the measurement interval. 
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Figure 2.2-4 GPS Satellite Visibility Criterion 


The criterion used for satellite visibility was the 
same in each case . To be visible , it was assumed that the 
satellite must be 10 deg or more above the horizon as viewed 
from the user location (Fig. 2.2-4). On the ground, the 
horizon is 90 deg below the zenith, so the half cone angle for 
satellite visibility is 80 deg. Above the earth on the Orbiter, 
however, the horizon is over 100 deg below the zenith, so the 
half cone angle for satellite visibility is increased to more 
than 90 deg. 

As Fig. 2.2-0 shows, it makes little difference, in 
terms of visibility, where the Orbiter measurement interval 
interfaces with the NESA output. Satellites No. 13 and 18 
disappear from view of the Orbiter near the end of the 90 
min window, but the remaining six satellites should be capable 
of providing good coverage. Based on Fig. 2.2-3, the Shuttle 
launch time was synchronized with 0 hrs in the GPS time frame. 
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2.3 FILTER-INDICATED PERFORMANCE 

The 11-state Kalman filter used during the GPS-aided 
ascent phase predicts pseudo-range and range-rate measurements 

/V 

(Piu* Pj^^) based on the before-measurement estimates of Orbiter 
position, velocity and clock states, combined with the satellite 
delivered estimates of the corresponding satellite states. Dif- 
ferences between the actual measurements P^^) ^be 

predicted measurements drive the filter and cause post-measure- 
ment updates in the Orbiter state estimates vsee Fig. 2.3-1). 

The filter gains used for the updates are based upon a filter- 
generated covariance matrix which reflects the filter's confi- 
dence in its before-measurement state estimates. The filter- 
indicated performance implied by this covariance matrix is 
summarized in this section. 


R-20621 



Figure 2.3-1 Typical GPS User Mechanization 

The general structure of the 11-state filter used 
in the ascent phase study is presented in Table 2.3-1. 
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TABLE 2.3-1 


NAVIGATION FILTER STATES 


State 

Number 

State Definition 

1-3 

Position* 

4-6 

Velocity* 

7-9 

IMU Misalignment Angles* 

10 

Receiver Clock Phase Error 

11 

Receiver Clock Frequency Error 


The initial system covariance matrix, Pj,(tQ), and continuous 
process noise matrix, Qf^» both diagonal. Initial 

standard deviations for the filter states and spectral 
densities for their associated process noises are presented 
in Table 2.3-2. 

The initial velocity error uncertainty specified in 
Table 2.3-2 is based on the preflight vehicle motion model 
given in Section 3.2 of Ref. 1. The initial misalignment 
and the velocity and misalignment process noise values are 
based on the IMU performance specification of Ref. 28. The 
values for the clock states (States 10 and 11) were selected 
under the assumption that the oscillator in the onboard re- 
ceiver will have operating characteristics similar to those 
of the Shuttle MTU and that it will be calibrated prior to 
launch. The MTU errors are discussed in Section 2.4. 

The measurement covariance matrix, R, for the GPS 
measurement is diagonal. RMS values of 12 ft and 0.05 fps 
are assumed for the position and velocity measurements — the 
12 ft value assumes that the P-code (high accuracy) GPS trans- 
missions are used on the Space Shuttle. Additional measurement 
_ 

These states are expressed in the I (Inertial) coordinate 
frame of Appendix A of Ref. 1. 


3S 
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TABLE 2.3-2 

ASCENT NAVIGATION FILTER INITIAL CONDITIONS 


T-OO-i-l 


State 

No. 

Description 

Initial Filter 
Standard Deviations 

Filter Process Noise 
Spectral Densities 

Value 

Units 

Value 

Units 

1 

X Position Error 

10 

ft 

— 

f t^/sec 

2 

Y Position Error 

10 

ft 

— 

ft 2/sec 

3 

Z Position Error 

10 

ft 

— 

f t^/sec 

4 

X Velocity Error 

0.01 

fps 

0.001 

( fps )2/sec 

5 

y Velocity Error 

0.01 

fps 

0.001 

(fps)2/sec 

6 

Z Velocity Error 

0.01 

fps 

0.001 

( ips)2/sec 

7 

X Platform Misalignment 

60 

sec 

0.035 

s^c“/sec 

8 

Y Platform Misalignment 

60 

sec 

0.035 

sec2/sec 

9 

Z Platform Misalignment 

60 

sec 

0.035 

sec^/sec 

10 

Receiver Clock Phase Error 

100 

ft 

0.1 

ft^/sec 

n 

Receiver Clock Frequency Error 

0.1 

fps 

3 10"^ 

(fps)2/sec 


errors, such as ionospheric propagation uncertainties, plus 
unmodeled satellite errors are assumed to be included in the 
12 ft and 0.05 fps values. 

As a preliminary study of filter performance, the 
filter was simulated twice, once as an optimal Kalman filter, 
and once with measurement underweighting. * The purpose of 
underweighting is to prevent the filter from making sudden 
large changes in the state estimates when a new sensor is 
first acquired. This reduces the possibility of filter 
divergence as well as reducing the sensitivity of the filter 
to a single bad measurement. In an optimal Kalman filter, 
the gain matrix for a measurement is computed as: 


K = PH^ (2.3-1) 

With underweighting, this equation is modified to read: 


♦Underweight ing is discussed in Ref. 1 and is included in 
the Shuttle baseline for entry (Ref. 45). 
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K = Ph’^Iu^qHPh'^+R]"^ (2.3-2) 

where is the underweighting factor, a scalar larger than 
unity. For the results presented here, an underweighting 
factor of 1.2 was used until the root sum of the position 
variances [the square root of the sum of P(l,l), P(2,2) and 
P(3,3)i was less than 1000 ft. 

When the filter was run with measurement underweight- 
ing, the underweighting was turned off after 3 measurements, 
but its influence continued to be significant until about 
the eighth measurement. After the 8 measurements, the 
optimal and underweighted filters had nearly identical per- 
formance. Since measurement underweighting did not appear 
to sacrifice any navigation accuracy at orbit insertion, it 
was used for the baseline filter. 

The filter covariance results are presented in 
Fig. 2.3-2 for the first 8 measurements. No measurements 
are made until 630 sec after liftoff; the filter covariance 
propagates in an open loop manner until that time. The 
first measurement is from Satellite No. 4, which has a 
line-of-sight approximately perpendicular to the Shuttle 
trajectory (Fig. 2.2-1). Consequently, the first update re- 
sults in a significant decrease in the filter-indicated rms 
crossrange position and velocity error. Subsequent measure- 
ments provide significant decreases in the rms downrange and 
radial errors. After eight updates, the filter-indicated 
rms errors have decreased from approximately 2400 ft to 20 ft 
per axis in position and from appro.ximately 7 fps to 0.2 fps per 
axis in velocity. Although measurements are taken until orbit 
insertion (771.0 sec, 24 updates), there is no significant re- 
duction in the filter-indicated covariance after the eighth 
update , 
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Figure 2.3-2 Filter-Indicated Ascent Navigation 

Performance for First Eight GPS 
Measurements 


The filter covariance data for all 11 filter states 
is summarized in Table 2.3-3. The data is presented at four 
times : 

• Liftoff (t = 0 sec) 

• prior to the first GPS measurement 
(t = 630" sec) 

• after the eighth GPS measurement 
(t = 672+ r'ic) 

• at orbit insertion (t * 771 sec) 

In addition to improving the position and velocity estimates, 
it is interesting to note that the filter-indicated results 
suggest that GPS-aiding during ascent can estimate IMU mis- 
alignments to an accuracy better than that obtained from 
ground alignment. Although the filter is somewhat optimistic 
in its estimates, this capability is confirmed in Section 2.5 


Jf 
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TABLE 2.3-3 

FILTER- INDICATED PERFORMANCE DURING ASCENT 


Event 

T 

AMS Position Errors 
(ft) 

RMS Velocity 
( fps) 

Errors 

RMS Misalicmnerit Errors 

(O) 

RMS Ciock Errors 

(sec) 

DR 

CR 

R 

DR 

CR 

R 

DR 

CR 

R 

Phas<^ 

Freq . 
(Jps) 

Launch 

0. 

10.0 

10.0 

30.0 

0.01 

0.01 

0.01 

60.0 

60.0 

60.0 

100.0 

0.10 

Before First 
UfasurcAent 

630. (-) 

2184.6 

2667 . 9 

1803,3 

4.96 

7.46 

7.19 

60.2 

60 2 

60.2 

128.4 

0.17 

After Eighth 
Ueasurement 

672.0 

e.e 

le.o 

28.6 

0.07 

0.16 

0.17 

12.0 

€.6 

12.4 

19 9 

0.10 

Orbit 

Insertion 

771.0 

5.0 

6.7 

12.5 

o.oe 

0.16 

0.13 

11.0 

6.8 

13.4 

7.5 

0.06 


2.4 ASCENT PHASE TRUTH MODEL 

The ascent phase truth model provides a detailed de- 
scription of the environment in which the navigation filter 
must operate. As such, it is the basic element of the SYSTEM 
module (Fig. 2-1) which is used to analyze the performance 
of the candidate filter. The error sources included in the 
truth model fall into five general classes: 

• Receiver measurement noise 

• IMU errors 

• Shuttle MTU errors 

• Gravity modeling errors 

• GPS satellite position, velocity and 
clock errors 

The mathematical details of the models for these error sources 
were presented in Ref. 1. A brief description of the models 
and a summary of the numerical values used for the performance 
analysis are presented in this section. 
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2.4.1 Truth Model Definitions 


With the exception of receiver measurement noise, each 
of the error models included in the ascent phase truth model 
requires that a set of error states be defined. A total of 128 
error states are included in the truth model; the states are 
identified in Table 2.4-1. Except for the eleven filter states, 
the error states are grouped in terms of the error source 
classes defined above. The total of 88 GPS-related states 
includes eleven for each of the eight visible GPS satellites. 


TABLE 2.4-1 


ASCENT PHASE TRUTH MODEL STATES 


Error Source 

Individual 

Error Sources 

To1 al 

Error Sources 

1. Filter States 

11 

11 

II. IMU Errors 


20 

Acceleromo; er 



Biases 

3 


Scale Factors 

3 


Asymmetries 

3 


Non-Orthoponal ities 

3 


Gyroscope 



Bias Dril’t 

2 


Mass Unba.iances 

5 


III. Shuttle MTU Errors 

6 

6 

IV. Gravity Model Errors 

3 

3 

V. GPS Satellite Errors 


8H 

Posi t ion 

3* 


Velocity 

3^ 


Clock 

3* 


Solar Pre.si^ure' 

1 ^ 


Gravi ta t i 1 Uncertai nt v 



Total 


12^ 


♦Number of brror Stales Per SatcHilt 


The eleven filter states are identical to those in 
the navigation filter mechanization; however, the position, 
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velocity and misalignment states are not driv'en by process 
noise as was implied in Table 2.3-1. Instead, these states 
are driven by the appropriate error sources in the truth 
model, e.g. , IMU and gravity model errors. (The filter treats 
these unmodeled errors as process noise.) Also, whereas 
states 10 and 11 reflect the filter estimates of the BiTU 
phase and frequency errors, the actual MTU errors are speci- 
fied in Group III of Table 2.4-1. 

The states defined for the IMU error model are a 
subset of those included in the detailed error model used 
for the preflight calibration and alignment analysis (Ref. 1). 
Only those IMU error sources which can make an appreciable 
contribution to the navigation errors at orbit insertion 
were retained. Thus, the original model was reduced to 20 
key states for this application. 

The gravity error model is based on Ref. 44 and con- 
sists of two vertical deflection states (downrange and cross- 
range) and one gravity anomaly state. The error magnitudes 
are typical of those expected for an onboard gravity model 
which does not account for high frequency gravitational ef- 
fects. The model parameters are determined as a function of 
altitude and velocity by table look-up. As altitude increases . 
the error magnitude decreases. 

The Shuttle MTU error model is the 5-state model 
developed in Ref. 1, except that an acceleration-sensitive 
error state has been added to reflect the most recent 
test data available (Ref. 50). The resulting 6-state model 
conforms with a NASA generated specification for normalized 




2-23 



the analytic sciences corporation 


rms phase error, and maximum random aging rate 

for the MTU. The specifications are: 


• Normalized rms phase error of 10“^® over 
the interval 0.5 sec i t £ 10 sec. 

-9 

• A maximum random aging rate of 10 per 

day. This is equivalent to an rms aging 
rate of 3.3 10”10 per day. 


The NASA specification and the 6-state approximation to this 
specification (assuming zero acceleration) are shown in 
Fig. 2.4-1. The model was developed based on the approach 
presented in Ref. 10. Figure 2.4-1 also demonstrates the 
normalized rms phase error for the 2-state model developed 
for the navigation filter. 



AVERAGING TIME - t (»ec) 

Figure 2.4-1 Normalized RMS Phase Error 

for Quartz Crystal Clock Model 
(Assuming Zero Acceleration) 


♦The MTU specification (Ref. 11) does not clearly indicate 
whether the oscillator performance is to be evaluated using 
Allan variances or normalized rms phase error. For a given 
numerical value, the latter is a less stringent specification 
and was therefore assumed for this study. 
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The linear system block diagram for the 6-state clock 
model is presented in Fig. 4.2-2. The parameters K. and Kn 

-1 -4 X ^ 

assume values of 1.01 x 10 and 1 x 10 , respectively, and a 

is the Shuttle acceleration along the major thrust axis. 

Table 2.4-2 contains a description of the six clock model states. 



Figure 2.4-2 Quartz Crystal Clock Error Model 

TABLE 2.4-2 

SHUTTLE MTU ERROR STATE DEFINITIONS 



0 
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The (g>S satellite error model was implemented in the 
NESA program (see Section 2.2), which produced a data tape 
with an 88 x 88 covariance matrix and an associated state 
transition matrix for each measurement time. Although the 
NESA program used more than 88 states in its error model 
(tracking station clock errors, survey errors, etc.), the 
88 states defined in Table 2.4-1 are sufficient to character- 
ize the GPS satellite contributions to user measurement errors. 

2.4.2 Truth Model Data Base 

The data base for the ascent phase truth model con- 
sists of the numerical data required to evaluate the filter 
performance in an environment defined by the truth model. 

In most cases this consists of initial covariances, process 
noise spectral densities, and measurement variances, but in 
certain instances it also includes time constants for corre- 
lated error sources. The data base for all error sources 
except the GPS satellite errors is presented in Table 2.4-3. 
The GPS satellite errors are characterized in the text of 
this section. 

The measurement errors specified in the table assume 
that a multichannel X-receiver is used to process the P-code 
(high accuracy) on the Space Shuttle. Additional measurement 
errors, such as ionospheric propagation uncertainties, plus 
unmodeled satellite errors are assumed to be included in the 
12 ft and 0.05 fps values. 

The standard deviations assigned to the IMU errors 
(including IMU misalignment) are based on the IMU performance 
specification (Ref. 28). The position and velocity uncertainty 
values were derived from a model of wind-induced oscillations 
of the Shuttle on the launch pad (Ref. 1). It is important to 
note that these error sources are modeled as independent 
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TABLE 2.4-3 


ASCENT NAVIGATION TRUTH MODEL DATA BASE 

T-0268 


Error Source 

Standard Deviation 

Time 

Constant 

Process Koise 
Spectral Density 

0. Uncorrelated C3>S Beceiver 




Keasureaeat Errors 




Pseudo^Baugtt 

12 

0 

0 

Pseudo-Bange-Rate 

0.05 fps 

0 

0 

I . Filter States 




Position 

10 ft 

m 

0 

Velocity 

0.01 fps 

m 

0 

Eisaligomeots 

60 gei 

• 

0 

11 . IWJ Errors 




Acceleroneter 




True Biases 

50 ug 

m 

0 

Scale Factors 

40 ppa 

m 

0 

Asyanetries 

40 ppo 

m 

0 

lloDortbogoaalities 

15 ^ 

m 

0 

Gyro 




True Bias Drifts 

0.015 deg/hr 



Bass Unbalance 

0.025 deg/hr /g 



111. Shuttle UTU Errors 




Phase Error 

100 ft 

— 


Frequency Error 

0.01 fps ^ 


/ 

Aging 

3.56*10-^ ft/sec^ 

«• 

> # #V Jl dh \ 

Rudott Frequency Error 

20.2 fps 


/ £q.(2.4«2) 

Randoa Frequency Bate Error 

1.37 ft/sec« 

i 

) 

Acceleration Sensitivity 

0.1 ft/sec*/g 



IV. Gravitational Deflections 

8 ^ec 

Function of 

Function of 

and Anomalies 

67 mgal 

Altitude and 

Altitude and 

! 


Vehicle Speed 

Vehicle Speed 

V. GPS Satellite Errors 

Computed Using 

See text 

See text 

j 

1 

KESA Prograas 




rather than correlated errors. The IMU errors resulting from 
the preflight calibration will be correlated, but the effort 
presently underway at TASC to estimate those correlations 
is not complete. When they become available, these correla- 
tions can be added to the truth model. 


Except for acceleration sensitivity, the values 

shown for the MTU are those used in generating Fig. 2.4-1. 

The MTU has a preferred axis along which the acceleration 

sensitivity is significantly smaller than for the other two 

2 10 

axes. In selecting the 0.1 ft/sec /g (1 part in 10 per g) 
value in Table 2.4-3, it was assumed that this axis will be 
oriented along the major thrust axis for the Shuttle, but 
the results of Section 2.5 indicate that, for the selected 
measurement schedule, acceleration sensitivity would not be 
a dominant error source even if the MTU orientation were 


arbitrarily chosen. 
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The process noises Wg and for the MTU model 
(Fig. 2.4-2) are linearly dependent. The process noise 


spectral density used 

4.56x10"^ 0 

0 0 

Q « 0 0 

4.14x10"^ 0 

4.51x10"^ 0 

0 0 

with units specified 


to generate Fig. 

0 4.14x10"^ 

0 0 

0 0 

0 3.77x10^ 

0 4.10x10"^ 

0 0 

in Table 2.4-2. 


.4-1 is 


4.51xl0‘^ 0 

0 0 

0 0 

4.10x10'^ 0 

4.47x10"^ 0 

0 0 


(2.4-2) 


The covariance matrix for the GPS satellite error 
sources was generated by the NESA program. The satellite 
contributions to measurement error vary as a function of 
Orbiter/satellite geometry, time since last ephemeris update, 
etc., but Fig. 2.1-2 provides a general indication of the 
magnitude of the errors involved. For satellites which were 
uploaded immediately prior to the launch time, e.g., satellite 
No. 18, the satellite-contributed measurement errors would 
be on the order of 2 ft and 0.001 fps. For satellites with 
less recent updates, the errors could be a factor of 2 or 3 
greater. 


2.5 GPS-AIDED ASCENT NAVIGATION PERFORMANCE 

The previous sections of this chapter have defined 
the proposed navigation filter for GPS-aided ascent navigation 


A/ 
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and a stochastic model (truth model) of the real world 
environment in which the filter is designed to operate. The 
performance of the navigation filter in this environment 
is analyzed via covariance analysis in this section. The 
results of the analysis are presented in the form of both a 
time history of overall navigation performance during the 
GPS measurement interval and as a detailed error budget at 
orbit insertion. 

The first GPS measurement is taken at 630 sec, midway 
through the final QMS insertion burn. Prior to that time the 
navigation mode is pure-inertial. The rms navigation errors 
at 630 sec are approximately 2000 ft and 7 fps along each 
axis, which roughly corresponds to the current specification 
for orbit insertion accuracy (Ref.17). The dominant contri- 
butor to this error is initial IMU misalignment, but various 
gyro and accelerometer error sources are also significant. 

The overall performance of GPS-aided navigation with 
the 11-state navigation filter is indicated in Fig. 2.5-1. At 
the time of the first measurement, there is an immediate 
decrease in both position and velocity errors. Because of the 
relative geometry of the Shuttle and the GPS satellite (Sat- 
ellite No. 4) for this measurement, the greatest initial 
reduction is in crossrange errors. As subsequent measure- 
ments from different satellites are processed, however, the 
improvement in all three components of position and velocity 
is both uniform and dramatic. After measurements from all 
eight visible satellites have been processed (at 672 sec), 
the rms navigation errors are on the order of 20-30 ft and 
0.2 fps along each axis, with the vertical errors being 
slightly larger because of the lack of a GPS satellite 
directly overhead of the Orbiter. Over the remainder of 
the measurement interval, the navigation errors slowly 
decrease so that they are less than 20 ft and 0.2 fps rms 
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a) RMS POSITION ERROR 


b) RMS VELOCITY ERROR 
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Figure 2.5-1 GPS-Aided Ascent Navigation: Overall Performance 


along each axis at orbit insertion. These errors are some- 
what larger than those of the GPS receiver and suggest that 
error sources other than the receiver errors (e.g., satellite 
ephemeris and clock) remain important. 


The conclusion to be drawn from Fig. 2,5-1 is 
that the proposed navigation filter can effectively incor- 
porate GPS measurements into the navigation solution. The 
degree to which GPS-aiding improves navigation performance 
is further illustrated by Table 2,5-1 which compares the 
navigation errors prior to the first GPS measurement with 
those at orbit insertion. The table indicates that GPS- 
aiding has decreased the rros position errors by better than 
a factor of 100 and the rms velocity errors by better than 
a factor of 30 . Furthermore, the table indicates that the 
filter has used GPS measurements to align the IMU during 
powered flight . Whereas the rms IMU misalignments were 
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TABLE 2.5-1 

EFFECT OF GPS IN REDUCING ORBIT INSERTION ERRORS 


T-0269 



RMS Positl 
Errors (ft 

.on 

> 

RMS Velocity 
Errors (fps) 

RMS 

Misalignment 
Errors ( ) 

RMS Clock 
Errors 

Implementation 

DR 

CR 

R 

DR 

CR 

R 

DR 

CR 

R 

Phase 

(ft) 

Fret}. 

(fps) 

Pure Inertial 
Navigation* 

2260 

2730 

1864 

5.2 

7.6 

7.4 

64 

62 

61 

132 

0.20 

GPS-Aided 

Navigation 

7 

9 

16 

0.1 

0.1 

0.2 

25 

13 

15 

4 

0.14 


^Errors at mid^oiot of final Orbital Uaneuverlug Syst^n (OHS) Insertion bum, 142 sec 
prior to Insertion 


slightly greater than 60 sec prior to the first measurement, 
they are approximately 13 s^ about the orthogonal (cross- 
range and vertical) axes and 25 s^ about the downrange axis. 
This capability for accurate inflight alignment could signi- 
ficantly reduce the importance of the current system specifi- 
cations for IMU alignment prior to launch and the deorbit burn . 

One measure of navigation filter performance is the 
degree to which the filter covariance reflects the actual 
navigation accuracy. Comparison of Table 2.5-1 with 
Table 2.3-2 indicates that the filter was somewhat optimistic 
relative to the true position and velocity accuracy, but the 
difference was not significant. The filter was considerably 
over-optimistic in its estimate of IMU misalignments, however, 
which suggests that it might be desirable to increase the 
process noise which the filter uses to account for gyro errors. 

The detailed error budget at orbit insertion is 
presented in Table 2.5-2. The most prominent feature of the 
table is that the navigation performance is completely domin- 


ated by GPS-related errors . The IMU is the major contributor 
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TABLE 2.5-2 


ASCENT ERROR BUDGET AT ORBIT INSERTION 


T-0270 



Position C 

ft) 

Velocity (fps) 

Misalignment {Seh) 

Error Sourco 

DR 

CR 

R 

DR 

CR 

R 

DR 

CR 

R 

Icltial Conditions 

0.0 

0.0 

0.1 

0.00 

0.00 

0.00 

1.6 

0.3 

2.7 

Measuresient Koise 



<so> 

0.04 

0.05 

0.06 

0.9 

0.4 

1.5 

Accelerometer : 










Quantization 

0.1 

0.1 

0.2 

0.00 

0.01 

0.01 

0.6 

0.3 

0.8 

Bias 

0.1 

0.6 

0.5 

0.01 

0.04 

0.03 

4.1 



Scale Factor 

0.0 

0.1 

0.2 

0.00 

0.00 

0.01 

0.3 

3.6 

0.2 

A83rmnetr7 

0.0 

0.1 

0.2 

0.00 

0.00 

0.01 

0.3 

3.7 

0.2 

Nonort bogonal i ty 

0.0 

0.1 

0.4 

0.00 

0.00 

0.00 


3.8 


Gyroscope : 









ceTto 

Bias 

0.0 

0.0 

0.1 

0.00 

0.00 

0.00 

Mass Unbalance 









cgTX^ 

Input Axis 

0.0 

0.0 

0.1 

0.00 

0.00 . 

0.00 

3.9 

Spin Axis 

0.0 

0.0 

0.0 

i 

0.00 

0.00 

0.00 

4.0 ^ 

1.1 

2.6 

MTU (Clock): | 



1 



Cg~T& 




Five State Model 

2.5 



CD . f > g > 

1.7 

1.1 

3.1 

Acceleration Term 

0.0 

0.0 

0.2 

0.00 

0.00 

0.01 

0.0 

0.1 

0.1 

Gravity Model 

0.0 

0.0 

0.0 

0.00 

0.00 

0.00 

0.8 

0.5 

0.3 

GPS Satellite Errors 


CC2> 


0.00 

o.co 

0.01 

0.2 

0.2 

0.3 

RSS Total 

6.9 

8.7 

16.0 

0.06 

0.12 

0.19 

25.7 

12.8 

14.7 


to the remaining misalignment uncertainty, but the velocity 
estimation accuracy is dominated by the frequency errors 
in the MTU and the position estimation accuracy is dominated 
by the combination of GPS receiver measurement noise and GPS 
satellite ephemeris (and clock phase) errors. This dominance 
of GPS-related error sources implies that navigation perform- 
ance is insensitive to an increase in the magnitude of any 
other error source as long as GPS is available . 

The most important conclusion to be drawn from 
Table 2.5-2 is that, when GPS-aiding is used, ascent navigation 
accuracy may not be sensitive to the quality of the IMU . The 
role of the IMU is to provide an inertial reference with short 
term stability so that GPS measurements taken at different 
times can be "tied together" to generate an accurate naviga- 
tion estimate. The IMU need not have good long term stability 
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performance characteristics and, more important in a Space 
Shuttle context , it need not be accurately calibrated or 
aligned . 


A second important observation is that the accelera- 
tion-sensitive MTU error had a negligible effect on the navi- 
gation accuracy. This is partly due to the fact that the GPS 
measurements were processed during a low acceleration phase of 
the mission (approximately 0.1 g) , but the error budget sug- 
gests that even during 3g acceleration the velocity error due 
to this error source would be less than 0.3 fps. Thus, the 
possibility exists for utilizing GPS-aided inertial navigation 
during the entire ascent mission phase . 

The importance of the MTU frequency errors in limit- 
ing velocity accuracy is reflected in Table 2.5-1. Although 
GPS-aiding improved the position and velocity estimates, it 
could not improve the clock frequency error estimates signi- 
ficantly, primarily because of short-term instabilities in the 
MTU frequency (see Table 2.4-3). In effect, the residual 
MTU frequency estimation error establish the lower limit for 
the velocity errors. Any attempt to further decrease the 
velocity errors would require an onboard clock with improved 
frequency stability characteristics. 


2.6 ASCENT SENSITIVITY ANALYSIS 

The error budget in Section 2.5 was generated under 
the assumption that the truth model presented in Section 2.4 
is an accurate statistical representation of the real world 
error sources which would affect navigation system performance. 
Since many of the models are based upon performance projections 
for hardware which has not yet been built, it is necessary 


J/ 
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to evaluate the sensitivity of the navigation filter to 
unknown variations in the rms values of the dominant error 
sources. Sensitivity curves for the ascent navigation filter 
are presented in this section. 

The error budget in Table 2.5-1 indicates that 
position errors at orbit insertion are dominated by GPS- 
related errors, primarily receiver noise and GPS satellite 
errors. The performance sensitivities to variations in these 
error sources are presented in Figs. 2.6-1 and 2.6-2, respec- 
tively. As an example of how to interpret the results, 
consider Fig. 2.6-1. The error budget indicates that, for 
the baseline error model, the radial position error at orbit 
insertion is 16 ft. If the receiver noise were a factor of 
5 worse than predicted by the error model, however, and if 
the filter parameters were not adjusted to reflect the 
increase, then the radial error woulo increase to 60 ft. 

Both Figs. 2.6-1 and 2.6-2 indicate a strong perfor- 
mance sensitivity to the GPS error sources. However, the 
error magnitudes anticipateu in the model are not likely to 
be significantly in error for Phase III GPS. If the Space 
Shuttle were to process the C/A-code (low accuracy) signal 
rather than the P-code (high accuracy) signal, the receiver 
errors would be a factor of 10 greater than indicated in 
the error model, but this would be a known degradation for 
which the filter could be adjusted. The resulting position 
errors would be greater than those indicated in the error 
budget, but quite likely would be considerably less than 
the xlO numbers indicated in Fig. 2.6-1. 

The only ma,,or contributor to velocity errors is 
the MTU frequency stability. As a result, the velocity error 
will exhibit a nearly linear dependence on unknown variations 
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Figure 2.6-1 RMS Position Error Sensitivity at 

Orbit Insertion to GPS Receiver 
Measurement Noise 



RELATIVE GPS SATELLITE ERRORS 


Figure 2.6-2 RMS Position Error Sensitivity at 

Orbit Insertion to GPS Satellite 
Tracking Errors 
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in the MTU errors. Figure 2.6-3 indicates that, if the MTU 
were replaced by a more accurate clock, the velocity errors 
could be significantly decreased even without adjusting the 
filter parameters. If velocity accuracies of 0.2 fps are 
adequate, however, the MTU will suffice as a time reference 
as long as it meets the established performance specifications. 



Figure 2.6-3 RMS Velocitv Error Sensitivity at 

Orbit Insertion to IITU Errors 

The principal contribution of IMU-related error 
sources to navigation performance at orbit insertion is to 
the residual misalignment error. This error will be important 
only if reentry is to be performed without GPS measurements 
and without realigning the IMU — constraints which potentially 
could be applied to short duration missions such as 3B. 

Because ■'f geometric effects, the largest residual misalign- 
ment is about the downrange axis. The results in Fig. 2 6-4 
indicate that, even if the IMU calibration were a factor of 
5 worse than expected, the downrange alignment error (also 
crossrange and radial) would still be less than 100 sec at 
orbit insertion. Since initial conditions are a negligible 
contributor to the error budget, the performance at orbit in- 
s'^ition would be even less sensitive to large initial azimuth 
misalignments at launch. 


I Sf 
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Figure 2.6-4 RMS Downrange Misalignment Error Sensitivity 

at Orbit Insertion to Gyro Errors 


2.7 GPS-AIDED ORBITAL AND DEORBIT NAVIGATION 
PERFORMANCE PROJECTION 


All major changes to the Space Shuttle orbit, in- 
cluding both the orbit insertion and deorbit maneuvers, will 
be accomplished with the OMS (Orbital Maneuvering System) 
engines. The thrust magnitudes as well as the burn durations 
will be similar for the final orbit insertion and the deor- 
bxt maneuvers (Table 2.7-1). Consequently, the orbit inser- 
tion maneuver and the techniques discussed in the previous 
sections can be used to gain insight into potential GPS-aided 
navigation performance for the deorbit maneuver. 

As discussed in Section 2.5-2, the position and 
velocity errors for GPS-aided navigation at orbit insertion 
are entirely due to GPS-related error sources. Since the 
GPS-related errors are not acceleration sensitive*, the 

* 

The .MTU does have an acceleration sensitive error compon- 
ent, but it was shown to be a negligible error source for the 
acceleration magnitudes of interest. 
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TABLE 2.7-1 

COMPARISON OF ORBIT INSERTION AND DEORBIT MANEUVERS 



MAXIMUM 

MANEUVER 

EVENT 

ACCELERATION 

DURATION 


(ft/sec2) 

(sec) 

Orbit Insertion 
Maneuver 
(Mission 3B) 

1.5 

255 

Deorbit 
Maneuver 
(Mission OFT-1) 

2.1 

212 


position and velocity estimation accuracy is essentially in- 
dependent of the presence (or absence) of OMS thrusting. Thus, 
the orbit insertion and deorbit maneuver performance evalua- 
tions can also be used to predict GPS-aided orbital navigation 
accuracy with two different sets of initial condition errors. 

Performance projections for both GPS-aided orbital 
and deorbit navigation are presented in this section. The 
necessary changes to the navigation filter and the truth 
model for deorbit are discussed in Section 2.7-1 and the 
performance estimates are presented in Section 2.7-2. 

2.7.1 Leorbit Maneuv^^er Study Baseline 

For the most part, the baseline used for the de- 
orbit study is the same as that used for the ascent study . 

The principal differences are in the duration of the mea- 
surement sequence, the parameter values in the navigation 
filter, and the initial conditions (both filter and truth 
model). These differences are summarized in this section. 
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The measurement schedule for the ascent study 
specified a total of 24 GPS measurements at 6 sec intervals 
over the final 138 sec of the insertion maneuver. Because 
the initial navigation errors will generally be worse for 
deorbit than for ascent, and because the navigation filter 
is suboptimal, the measurement schedule for deorbit was ex- 
tended to 36 measurements over 210 sec — approximately the 
duration of the deorbit maneuver. 

The initial conditions for orbit insertion naviga- 
tion are a consequence of IMU cal /align and MTU errors pro- 
pagated through ascent. By contrast, the initial conditions 
for deorbit navigation are a function of the most recent 
navigation fix. For the deorbit study, the initial position, 
velocity and alignment errors assumed for the truth model 
are the same as those used for entry in Chapter 4 (Ref. 47). 
These errors correspond to an IMU alignment with the onboard 
star tracker and an uploaded navigation fix, both two hours 
prior to the deorbit maneuver, and then open loop propagation 
to the time of the maneuver. The resulting navigation errors 
are primarily in-plane errors — approximately 12,000 ft and 
13 fps rss. The misalignments are 285 s'ec per axis. 

The JITU initial conditions assume an accurate 
calibration at the time of the star tracker alignment and 
subsequent open loop propagation. The calibration errors as- 
sumed are those existing at orbit insertion for the GPS- 
aided ascent study (Table 2.5-1). At the time of the deorbit 
maneuver, the resulting MTU errors are 1129 ft in phase and 
0.25 fps in phase rate. 

The initial filter covariance matrix and the filter 
parameter values are shown in Table 2.7-2. The covariance 
is diagonal except for correlations between downrange and 

2-39 


■^7 



THE ANALYTIC SCIENCES CORPORATION 


TABLE 2.7-2 

NAVIGATION FILTER INITIAL CONDITIONS 


T-0M4fc 


State 

Bo. 

I>&seriptio& 

Initial Filter 
Standard Deeiationa 

Filter Proeeas Boise 
Spectral Densities 

Value 

Units 

Value 

Orlts 

1 

X Posltlott Error 

16.700 

ft 

— 

ft^/sec 

2 

T Position Error 

16.700 

ft 

— 

ft^/sec 

3 

Z Position Error 

16,700 

ft 

— 

ft^/sec 

4 

X Velocity Error 

18.7 

fps 

0.001 

(fps)^/sec 

5 

T Velocity Error 

18.7 

fps 

0.001 

(fps)^/8ec 

6 

Z Velocity Error 

18.7 

fps 

0.001 

(fps)^/sec 

7 

X Plat fora ICisaligsment 

900 

sec 

0.035 

sec^/sec 

8 

T Platform Ifisallgnment 

900 

sec 

0.035 

sec^/sec 

9 

Z Platform Itisalignment 

900 

sec 

0.035 

sec^/aee 

10 

Beceiver Clock Phase Error 

1130 

ft 

0.1 

ft^/sec 

11 

Beceiver Clock Frequency Error 

0.25 

fps 

3 * 10*^ 

(fps)^/sec 


radial errors*. Except for the addition of the clock errors, 
the filter data base is the same as that used in Chapter 4. 
These values are somewhat pessimistic relative to the truth 
model errors, particularly in crossrange position and velo- 
city, in order to account for the effect of unmodeled errors. 
When coupled with measurement underweighting, the net effect 
of this pessimism is to require a longer measurement sequence 
than would be required if a more optimistic initial filter 
covariance were selected. 


2.7.2 Performance Projections 


The error budget in Section 2.5 indicated that posi- 
tion and velocity errors for GPS-aided nav'igation were due 
almost exclusively to : 


*Downrange and radial navigation errors will be coupled due 
to the effect of orbital dynamics — a dov/nrange velocity 
error at orbital perigee will generate a radial position 
error at orbital apogee, etc. The filter accounts for this 
effect by assuming -1 correlations between downrange veloc- 
ity and vertical position and between downrange position and 
vertical velocity. 
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• measurement noise 

• MTU errors 

• GPS satellite errors 

In addition, if "initial conditions" are redefined to mean 
the position, velocity and misalignment errors at the begin- 
ning of the measurement sequence, including the effect of IMU 
errors up to that time, then the misalignment errors at orbit 
insertion are dominated by: 

• initial conditions 

Because of the degree to which they dominate the ascent study- 
results, these fo’ur error groups are sufficient to character- 
ize navigation performance for the orbital and deorbit navi- 
gation studies. The performance projections summarized in 
this section are based on these error groups. 

The results of the deorbit study are presented in 
Table 2.7-3. Even though the true initial errors were 
considerably larger than for the ascent study, and the filter 
initial covariance was not representative of the true errors, 
the position and velocity errors at the end of the deorbit 
measurement sequence were virtually identical to those at 
orbit insertion. This suggests that, at the very least , 
the 11-state filter can yield good navigation performance 
over the range of operating conditions expected for exoatmos- 
pheric mission phases . A future effort will investigate the 
filter performance during atmospheric flight. 

The disappointing aspect of the deorbit results 
is the inability to estimate IMU misalignments. Whereas the 
misalignments were estimated to better than 25 s'ec for the 
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TABLE 2.7-3 

EFFECT OF GPS IN REDUCING DEORBIT ERRORS 
T-0462 


EVENT 

RMS POSITION 
ERRORS (ft) 

RMS VELOCITY 
ERRORS (fps) 

RMS 

MISALIGNMENT 
ERRORS (^) 


DR 

CR 

R 

DR 

CR 

R 

DR 

DR 

R 

Begin Deorblt 

11265 

866 

1036 

1.1 

0.9 

12.4 

285 

285 

285 

End Deorbit 

7 

9 

17 

0.1 

0.1 

0.2 

283 

256 

257 


ascent study, there was minimal improvement in the estimates 
during the deorbit measurement sequence. The principal dif- 
ference is that the misalignments were highly correlated with 
the position and velocity errors during ascent, and hence 
readily recoverable from range and range-rate measurements. 

By contrast, the navigation errors during orbital operations 
are independent of IMU misalignments. Fven though the deorbit 
acceleration will introduce a small correlation between mis- 
alignment and velocity errors, the table suggests that the 
acceleration level w'ill not be large enough to make accurate 
(60 sec) misalignment estimates unless the duration of the 
deorbit maneuver (and measurement sequence) is increased.* 

Although IMU misalignment errors apparently cannot 
be significantly reduced during deorbit, it is possible that 
the misalignments can be estimated with sufficient accuracy 
to provide a consistency check of the IMU alignments. Th is 
capability could eliminate the redundancy management problem 
of detecting the occurrence of major shifts in gyro parameters 
in the two hours between IMU alignment and deorbit. A more 
thorough investigation is required to determine the perfor- 
mance capability of such a consistency check. 

*The true acceleration level for deorbit will be somewhat 
greater than that used for this study (2.2 vs 1.5 ft/sec^), 
but this would not significantly improve misalignment esti- 
mation accuracy. 


(,a 
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The present baseline for operational Space Shuttle 
missions specifies the Tracking and Data Relay Satellite 
System (TDRSS) as the principal navigation aid for orbit 
determination. TDRSS can provide both range and range-rate 
information, but it suffers a number of operational disadvan- 
tages relative to GPS. Because TDRSS consists of only two 
satellites, TDRSS/user geometry is a more important factor 
in orbit determination accuracy than it would be for Phase III 
GPS. Also, utilization of only two satellites means that 
longer measurement sequences must be processed in order to 
generate a navigation solution. Finally, orbit determination 
using TDRSS requires ground processing and subsequent upload- 
ing of the navigation solution. This will tend to decrease 
the operational flexibility of TDRSS-aided navigation rela- 
tive to an autonomous GPS navigation system. 


Table 2.7-4 compares orbit determination accuracy 
for TDRSS and GPS. The TDRSS data was taken from Ref. 15;* 
the GPS perfonnance projections are based on the orbit in- 
sertion and deorbit results after two min of data processing. 
The table clearly implies that GPS can prov'ide significantly 


better navigation accuracy than TDRSS from shorter measurement 


sequences . In contrast to TDRSS, which requires a relatively 
long trajectory arc to estimate orbital parameters (5 n.in was 
the shortest arc considered in Ref. 15), GPS could attain the 
indicated accuracy with only a few seconds of data if the 
navigation filter were suitably optimized. 


♦The error analysis presented in Ref. 15 considered only re- 
ceiver noise. It did not consider the effects of receiver 
oscillator instability, user dynamics, etc. As a consequence, 
the results in Ref. 15 are perhaps overly optimistic. 
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TABLE 2.7-4 

PROJECTED ORBIT DETERMINATION ACCURACIES FOR 
TDRSS AND GPS 


NAVIGATION 

AID 

LENGTH OF 
MEASUREMENT 
SEQUENCE 
(»nin ) 

RSS POSITION 
ERRORS (ft) 

RSS VELOCITY 
ERRORS (fps) 

TDRSS 

5 

200* 

Not Specified 

GPS 

2 

20 

0.2 


*See footnote on previous page. 


\ 
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3. CALIBRATION AND ALIGNMENT COVARIANCE ANALYSIS 


The performance analysis for the Space Shuttle cali- 
bration and alignment algorithm is being conducted in two 
phases. The first phase, completed in the previous contract 
period (Ref. 1), used monte carlo techniques to generate an 
overall performance projection for the Approach and Landing 
Test (ALT) algorithm (Ref. 51) and'to identify potential risk 
areas. The second phase, currently in progress, uses covari- 
ance analysis techniques to develop a detailed error budget 
for the Orbital Flight Test (OFT) algorithm (Ref. 18 ) and to 
identify the major error mechanisms which limit calibration 
and alignment performance. For those error parameters for 
which the algorithm cannot provide the required calibration 
accuracy, the results of this phase will provide a clear in- 
dication as to the modifications which must be made in order 
to achieve the desired performance. 

An overview of the calibration and alignment mechani- 
zation for OFT is presented in Section 3.1. Section 3.2 con- 
tains a general description of the mathematical structure 
used in the covariance simulations and Section 3.3 summarizes 
the truth model states and error sources used in the covari- 
ance analysis. A partial error budget for the cal/align 
software is presented in Section 3.4. 


3.1 CALIBRATION AND ALIGNMENT OVERVIEW 

The calibration and alignment of the Space Shuttle 
inertial measurement units (IMUs) consists of two parts; 
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the hangar calibration and the preflight calibration and 

alignment. The hangar calibration will be performed in a 

low-vibration environment within 5 to 14 days prior to 

launch, and will have a total duration of 8.6* hr. The pre- 

flight calibration and alignment will be initiated at the 

launch pad within 3 to 15 hr before launch, and will have 

* 

a total duration of 2.8 hr. The calibration and alignment 
schedule during both the hangar and the preflight mission 
phases is illustrated in Fig. 3.1-1. 

The cal/align procedure consists of steering the 
platform t'' a series of discrete positions at which torquing 
rates are applied to compensate for earth rate, and in some 
cases ’’excess" rates are applied to torque the platform 
through a prescribed trajectory relative to a local vertical 
(North, West, Up) coordinate system. At each of these posi- 
tions, the accelerometer outputs ( in a few cases also the 
resolver outputs) are sampled and a least squares fit (LSF) 
filter is used to estimate the acceleration along the up axis 
and the platform tilts and tilt rates about the north and 
west axes. Once the data (vertical accelerations, tilts and 
tilt rates from the various positions) necessary to evaluate 
the required IMU error parameters has been collected, these 
are computed and the compensation pareimeters are updated. 
Between discrete positions, the platform is slewed at a high 
rate and no data is collected during this motion. 

The cal/align mechanization comprises six different 
sequences. Each of these consists of several positions and 
is used for the calibration, alignment or verification ol 
several IMU errors as outlined in Table 3.1-1. 


♦These times are approximate and subject to modification. 




3-2 



THE ANALYTIC SCIENCES CORPORATION 




CAL/ALtGN SEQUENCES 


DURATION (hours) 




I. HANGAR CALIBRATION 

INITIALIZATION. WARM-UP 
COARSI ALIGNMENT 
HANGAR CAL A 
HANGAR CAL B 
HANGAR CAL C 

IL PREFLIGHT CALIBRATION 

INITIALIZATION. WARM-UP 
COARSE ALIGNMENT 
PREFLIGKT CAL A 

HI. PREFLIGHT ALIGNMENT 

GVROCOMPABSING 
VELOCITY, tilt INITIALIZATION 


0.6 


6.5 


l.3€ 


as 


1 72 


05 




6-14 DAYS BEFORE LAUNCH 


'A 
3-16 HOURS 
BEFORE LAUNCH 


O 

LAUNCH 


Figure 3.1-1 Calibration anii Alignment Schedule 


TABLE 3.1-1 

ERRORS CONSIDERED DURING THE CAL/ALIGN PHASE 


T-004C 


1 IVV Errors 


gar Phase 

Pre: light Pnase 


!\o . of 

Cal A 

Cal B 

Cal C 

* 

Gyrc» 

ITI 

Nar.e 

Errors 





compass 


Atcel. Bias High Gain 

3 







Accel. Scale Factor High Gain 

3 

1 






Acccl. Bias Lou- Gale 

3 

1 



1 



Accel. Scale Factor Low Gain 

3 

1 



1 



Gyro Bias Drift Rate 

3 

1 



1 



Gyro Mass Inba lance 

5 

1 



4 



Accel. A&yr.nietry 

3 

1 






Accel. Noncrthogonality 

*i 

1 






Gyro Tcrcuer Scale Factor 

3 

1 






Gyro Misalignrter.t 

6 

1 






Resolver Offset 

3 


I 





Global Xonortncgonality 

1 


1 





Accel ar.d Gyro to Platform 








Misalignments 

2 


1 





Relative Cluster Attitude 

• 



3 


4 


Tilt Misalignments 

2 





5 

5 

Azimuth Misaligxunent 

1 





5 


Key: 1 - Compete and calibrate 






2 - Compute and, if commanded, calibrate 



3 - Co.Tpjte and store 





4 - Compute ana check for 

reasonabl 

eness 



5 - Compute and align 
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3.2 COVARIANCE ANALYSIS PROGRAM DESCRIPTION 
3.2.1 Program Overview 

The structure of the covariance analysis program is 
the same as for the OFT IMU software. The program is orga- 
nized into seven separate sequences with an executive program 
which defines the interface between those sequences (Fig. 3.2-1). 
The input to the Coarse Alignment sequence is a data file 
which defines the IMU status (values for all error coefficients, 
gyro torquing rates, etc.), a priori calibration and alignment 
estimates, and the transformation matrices between the plat- 
form coordinate system and the vehicle-, earth- and inertially- 
fixed coordinate systems. The output of the Coarse Alignment 
sequence is a new data file updated by the Coarse Alignment 
estimates. This file is then the input to Hangar Cal A and 
so on. 


The program preserves the correlation between esti- 
mation errors in different sequences and has the capability 
of simulating a complete calibration and alignment procedure 





Figure 3.2-1 Structure of Covariance Simulation 

for Performance Analysis 
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from beginning to end. Alternatively, the initial condition 
data for any sequence can be used (and modified if desired) 
for performing special studies for a given sequence. 

3.2.2 State and Covariance Equations 

The LSF filter used in the cal/align software is a 
recursive filter, i.e., the tilt and tilt rate estimates 
generated at a particular time are a functio- only of the 
current accelerometer outputs and the estimates generated at 
the previous time. This recursive property permits the esti- 
mation procedure to be modeled as a (suboptimal) state esti- 
mator for the IMU error parameters based on a linearized model 
of the system dynamics. Because of the fact that the IMU 
parameter estimates are updated only at the end of each cal/ 
align sequence, however, it is necessary to introduce tem- 
porary storage states into the filter equations. The form of 
these equations and th approach used in the program to gen- 
erate the accompanying covariance equations are discussed in 
this section. 


The state and covariance error equations are sum- 
marized in Table 3.2-1. Equation (3.2-1) found in the table 
models the error propagation for the linearized IMU system. 
The equation represents the fundamental operational mode of 
the IMU including its associated software during the hangar 
and preflight phases (the Ground Operation Sequence, Ref. 18). 
The state vector x is composed of three parts: 


X.. - the least squares fit (LSF) filter states 

used to estimate the platform accelerations 
(tilts) and acceleration rates (tilt rates) 
by measuring the platform velocity outputs, 
and/or to estimate the platform attitude and 
attitude rate by measuring the gimbal angles 
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TABLE 3.2-1 

SYSTEM AND COVARIANCE ERROR EQUATION SUMMARY 


1 


CAL/ALIGN STEP 

STATE EQUATION 

COVARIANCE EQUATION 

ESTIMATE TILTS AND DRIFTS 
PROPAGATE IMU ERRORS 

X • F X ♦ w (3.2-1) 

P « r P 4 p F*^ 4 Q (3.2-2) 

TORQUE OUT TILTS 
RESET FILTER 

X* - L x' (3.2-3) 

P* •> L P" (3.2-4) 

STOKE F.ST1 MATES 
RESin* FILTER 

x"^ • C X* (3.2-5) 

P* - C P' (3.2-C) 

UPDATE IMU PARAMETERS 1 

RESirr STORAGE STATES 

x**^ • D x' (3.2-7) 

p* . D P* (3.2-f.) 


whfre the vectors art defined by 





LSF Filter States 


Si 


Measurement Noise 

X * 

-2 

m 

Storage States 

W «i 

0 

e 

- 


-3 


IMV Error States 


-3 


Process Noise 


and the matrices are given by 


JKI 

ERROR COVARIANCE 


0 0 KH 


P P P 

11 12 *13 

y- 

0 0 0 

P «• 

P P P 

*^21 22 ^23 


0 0 F3 


P P P 

*^31 32 ’ 33 

PLATFORM LKVELING 

DATA STORAGE 


NOISE COVAHIANCE 

k 


11 

0 


<*‘33 


0 0 
PARAMETER I'PDATL 


(3.2-fl) 



0 

0 

0 


0 

0 

0 


0 

0 

0 

L ® 

0 

I 

0 

c - 

s 

1 

s 

I> - 

0 

'>2 

0 



0 

I 


0 

0 

1 


0 

-»3 

1 


(3 2-10) 


x« - the dummy or storage states, used to store 
some of the filter estimates, , which 
are required to update IMU error parameters 
at the present or at a future time, 

x_ - the system states, including the platform 
velocity errors, misalignments and all IMI' 
error states being considered in the present 
cal/align evaluation. 
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Equation (3.2-1) is driven by white process noise repre- 
sented by the vector w. The components of w are resolver 
measurement error and the IMU process noise w^. There 
arc no errors associated with the transfer of filter states 
into storage states Xg. The propagation of the statistics 
associated with the state vector x is governed by Eq. (3.2-2). 
The details about the elements in and H are given in 
Appendix B of Ref. 7. Details about K, Lg, Cg, Dg and Dg 
are given in Appendix D of Ref. 1 and in Appendix B of this 
report, and about Q and P(0) in Section 3.3.2 of this report. 

At the end of each estimation period some of the 
filter estimates x^^ are stored in the temporary memory states 
Xg and, simultaneously, the filter states are reset to zero. 
This storage and reset is represented by Eqs. (3.2-5) and 
(3.2-6). Depending on the particular point in the cal/aiign 
process, the next step will be one of the following: 


• Updat-^ one or more of the IMU errors, X 3 , 

using ;ne or more of the estimates stored 
in X 2 > reset the storage states in X 2 

which are not required for future IMU error 
updates. This step is represented by Eqs. 
(3.2-7) and (3.2-8). 

• Rotate the platform about the horizontal 
axes by an amount equal and opposite to the 
filter estimates of the platform m'.salign- 
ments from level. This seep is represented 
by Eqs. (3.2-3) and (3.2-4). 

• Initiate the next estimation of platform 
tilts and drifts and IMU error propagation, 
as described above and represented by Eqs. 
(3.2-1) and (3.2-2). 


Ttie cal/align evaluation performed by TASC is based 
on computation of the covariance Eqs. (3.2-2), (3.2-4), 
(3.2-6) and (3.2-8). Since the corresponding state Eqs. 
(3.2-1), (3.2-3), (3.2-5) and (3.2-7) are all linear, it is 
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possible to initialize the matrices P and Q with only those 
error sources of interest at a particular time and to obtain, 
at the end of the run, the sensitivity of the cal/align pro- 
cess to that particular error source or group of errors. 


3.3 TRUTH MODEL DESCRIPTION 

Generation of an error budget for the cal /align 

* 

software requires that a "truth model" describing the real 
world error sources for the cal/align process be defined. 

The description must include a complete list of states and 
error sources, and a data base. The truth model to be used 
to evaluate the Space Shuttle cal/align software is presented 
in this section. 

3.3.1 States and Error Sources 

The truth model states and error sources used in 
evaluating the IMU cal/align mechanization are listed in 
Table 3.3-1, which divides them into three major categories: 

• Category I - Platfonn misalignments and 
calibrated IMU error states 

• Category II - Noncalibrated IMU error states 

• Category III - Random and quantization errors 

These three categories contain the error sources 
which will be used tc generate the baseline error budget . 

The first category corresponds to the platform misalignments 


♦The "truth model" is a mathematical model of all potentially 
significant error sources and the way they affect the cal/ 
align performance in the real world. 


76 
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TABLE 3.3-1 

I MU TRLTH MODEL STATES AND ERROR SOURCES 


T-of>n 


ERROR SOURCE NAME 
AND 

group NUVOEH 

NTMHKR 

or 

STATES 

NUMBER 
OP ERROR 
SOUPCLR 

ERROR 

M.NKM0MC 

1. PLATFORM! MISALlGNIfFNTS 
AND CALIBRATED ERRORS 




1. 

Platfonn Misalignmcots 

3 

0 


2. 

Gyro Bias Drifts 

3 

3 

DFK,Y.Z 

3. 

Gyro Mass Unha lances 

5 

5 

DlX.Z.DSX.Z.lXlZ 

4. 

Gyro Torquer Scale Factors 

3 

3 

KTX,Y,Z 

5. 

Gyro. Ilisalignments 

6 

6 

BXY.XZ ZY 

6. 

Accel. Biases - Lov Gain 

3 

3 

KOX.Y.Z 

7, 

Accel. Scale Factors - Low Gain 

3 

3 

KlX.Y.Z 

a. 

Accel. Biases - High Gain 

3 

3 

K0ilX,V.Z 

9. 

Accel. Scale Factors - High Gain 

3 

3 

ElHX.i ,Z 

10. 

Accel Scale Factor Asyraactries 

3 

3 

ESX,Y.7 

11 

Acccl . Honor t hogona lit ies 

3 

3 

DELTYX , Z.X , 

12. 

Accel, and Gyro Misalignments 

2 

2 

MXZ.MYZ 

13. 

Resolver Offsets, Giiabal 
Honor thogonal Ities 

4 

6 

lRO.PO,AZO.DP,a,b 

11. NOKCALIBKATED ISfU ERRORS I 




14. 

Gyro Anisoelasticities 

n 

11 

D12X,Z.n.S2\,Y .2 
D1SX,Z.D0S\.Z,D10X,2 

15. 

Gyro Output Axis Mass Unbalances 

2 

2 

DOX.Y 

16. 

Gyro Heading Sensitive Drifts 

3 

3 

RSN,« ,U 

17. 

Gyro Thermal Transient Drift 

3 

3 

TTX.Y.Z 

la. 

Accel. Monlinearit ies - 2°^ Order 

3 

6 

K2X,Y,2.X2,yz,7.X 

19. 

Accel. Konllnearit ies - 3**^ Order 

3 

3 

K3X,Y,Z 

20. 

Acccl. Heading Sensitive Biases 

3 

3 

KORSN.h.U 

21. 

Acccl. Heading Sensitive Scale Factors 

3 

3 

K1HSN,%,U 

22. 

Outer Bell Off. set and Misalignments 

3 

3 

ORO.AY.&Z 

23. 

Resolver Harmonics 

24 

24 

jSt^^th 

111. RAKDOV AND QUANTIZATION ERROR.S 




24. 

Gvro Randomness and Quantisation 

3 

6 


25. 

Accel. Randomness and Quantization 


6 


26. 

Resolver Randomness and Quantization 


7 


27. 

Vehicle Motions 

7 

2 


T0TAI5 

115 

124 



and those IMU error sources which are calibrated during one 
or more calibration sequences. Category II corresponds to 
additional IMU error sources which are not estimated or cali- 
brated during either the hangar or the preflight calibrations. 
Category III contains those states which are driven by random 


7 / 
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errors, including instrument randomness, quantization effects 
and vehicle motions induced by wind gusts. 

All error source categories are divided into smaller 
sets, each of which is associated with a group number. These 
group numbers were quite useful in defining the detailed truth 
model equations in Appendix B of Ref. 7. In addition, each 
line of the error budget that will be generated in this study 
will correspond to the contribution from all error sources in 
one particular group rather than each error taken one at a 
time. This simplifies the error budget table and places the 
various error source groups in better perspective. 

3.3.2 Truth Model Data Base 

Generation of a detailed error budget requires the 
numerical values for all truth mod-'l e’’ r sources. Two sets 
of errors must be specified: the initial calibration error 
(instrument error minus software compensation term) at the 
beginning of the hangar calibration phase, and the changes 
of the instrument errors during the time elapsed between the 
hangar calibration and the preflight calibration due to the 
turn-off turn-on errors. This data is summarized in Tables 
3.3-2 and 3.3-3. 

The numerical values described above are required to 
construct F, P and Q in Eq. (3.2-2). The specific elements 
of F, P and Q which must be assigned are: 

P33(0) - the initial IMU error covariance matrix 

P„ 3 (l) - the IMU turn-off turn-on instability 
error covariance matrix 

Qgg - the process noise covariance matrix 
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TABLE 3.3-2 

IMU TRUTH MODEL DATA BASE FOR CALIBRATED ERRORS 


T-0042 


Error Source Name 
and 

Group Number 

Standard Deviation 

Data 

Source 

Units 

At Beginning 
of Hangar Cal 

Turn-off Turn-on 
Instabilities 

2. Gyro Bias Drifts 

sec /sec 

0.035 

0.016 

Ref. 1 

3. Gyro Uass UabaXances: 





Input Axis Acceleration 

s"^/sec/g 

0.3 

0.015 

Ref. 1 

Spin Axis Acceleration 

s'^/sec/g 

0.02 

0.005 

Ref. 22 

Output Axis Acceleration 

j TOC / sec /g 

0.025 

0.01 

Ref. 22 

4. Gyro Torquer Scale Factors 

ppm 

400 

100 

Ref. 1 

5. Gyro liisalignments 

sec 

120 

15 

Ref. 1 

6. Accel. Biases • Low Gain: 





X.Y 


100 

50 

Ref. 1 

Z 

ug 

200 

50 

Ref. 1 

7. Accel. Scale Factors - Low Gain 

ppm 

200 

100 

Ref. 22 

8. Accel. Biases • High Gain: 





x.y 


100 

50 

Ref. 1 

z 


200 

50 

Ref. 1 

9. Accel. Scale Factors - High Gain 

ppm 

200 

100 

Ref. 22 

10. Accel. Scale Factor Asymmetries 

ppm 

200 

0 

Ref. 22 

11. Accel. Nonorthogonalities 

sec 

60 

10 

Ref. 1 

12. Accel, and Gyro Uisalignments 

sec 

60 

1 

10 

Ref. 1 

13. Resolver Offsets 

1 sec 

200 

20 ; 

Ref. 22 

Gifflbal Nonortbogonalities 

s?c 

60 

° 1 

Ref. 22 



the IMU truth model error dynamics matrix 
the measurement error covariance matrix 


Elements of Pgg, Qgg are normally chosen together to 

define random processes with desired properties. For instance, 
a first-order markov process (exponentially correlated random 
process) is modeled with a specified rms value, o^, and cor- 
relation time, T^. In this case the relevant truth model 
matrix elements are; 


P 


33 . . 


1 1 


2 

0 . 

1 


(3.3-1) 


V 
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TABLE 3.3-3 

IB4U TRUTH MODEL DATA BASE FOR NONCALIBRATED ERRORS 


T-0043 



Error Source Name 

Units 

Standard 
Deviat ion 

Data 

Source 

14. 

Gyro Anisoclasticities: 

Input*Input Anisoel. Drift 

^ 2 
sec / sec /g 

0.005 

Ref. 

22 


Spio-Spin Anisoel. Drift 

sec/sec/g^ 

0.005 

Ref. 

22 


Input*Spin Anisoel. Drift 

i— . 2 

sec/sec/g 

0.025 

Ref. 

22 


Output^Spin Anisoel. Drift 

, ,2 
sec/sec/g 

0.025 

Ref. 

22 


Input-Output Anisoel. Drift 

sec/sec/g 

0.005 

Ref. 

22 

15. 

Gyro Output Axis Mass Unbalances 

s^/sec/g 

0.005 

Ref. 

22 

16. 

Gyro Heading Sensitive Drift 

^c/sec 

0.005 

Ref. 22 

17. 

Gyro Thermal Transient Drift 

sec /sec 

0.0064 e-t/feOOO) 

Ref. 

22 

18. 

Accel. Non linearities: 
2"“ Order 

ug/g* 

15 

Ref. 

22 


Product Nonlinearity 


10 

Ref. 

22 

19. 

rd 

Accel. Nonlinearities - 3 Order 

ug/g^ 

5 

Ref. 

22 

20. 

Accel. Attitude Sensitive Biases 


10^2) 

Ref. 

1 

21. 

Acccl. Attitude Sensitive Scale Factors 

ppm 

15<2> 

Kef. 

1 ' 

22. 

Outer Roll Offset and Misalignments 

st:c 

200 

Ref. 

1 

23. 

Resolver Harmonics* 

Harmonic i 

sec 


Ref. 

22 


8^^ Harmonic 

oec 

19 0 

Ref. 

22 


9^*^ Harmonic 

sec : 

4.2<2) 

Ref. 

22 


16^^ Harmonic 

S 

sec 

20.0 

Rof . 

22 

24. 

Oyro Randomness 

/ 3/2 

sec/sec 

0-0001 

See text 

25. 

Accel . Output Noise 

yg 

5 

Ref. 

22 


Velocity Quantization - Low Gain 


0 0003 

See t-^xt 


Velocity Quantization - Hiph Gain 


0. 00003 

See text 

26. 

Resolver Output Noise 

sec 

i 

12 1 

See text 


Resolver Quantization 

sec 

5 8 i 

See text 

27 

Vehicle Motions 


See tf'M 

See text 


(1) t * seconds after comp 1 Mine the 30 second dipat^ii torquinp lollouinc a platform 
s leu . 

(2) Standard deviation of maxintrt amplitude. 


« 33 ,. - (3.3-3) 
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A constant error source can be considered a special case of 


the above with x . = «= and F__ 

1 Jd. . 


%3ii ~ 


The gyro randomness error assumes a random walk 

—6 2 

model with a growth rate of 32 x 10 (deg/hr) /hr. This vari- 

^ 3/p 

ance corresponds to a standard deviation of 0.0001 sec/sec ' 
as given in Table 3.3-3. 


The velocity quantization error with the acceler- 
ometers operating in the low gain mode is determined as follows 

One AV pulse = 0.0343 fps (3.3-4) 


Assuming that each AV pulse count has a random error selected 
from a uniform distribution ranging from -0.0172 fps to 
+0.0172 fps, and one such error is introduced every 0.32 sec 
computational cycle*, the error variance is 

Qg., = (0.0343)^/12/0.32 (3.3-5) 

''ii 

= 0.0003 (fps)^/sec (3.3-6) 


The quantization level for the resolvers is 20 sec. 
Because dithering rates are applied to all axes along which 
successive resolver measurements are made, this quantization 
error may be represented as an uncorrelated measurement noise 
of 


Rll = (20)^/12 = 33.3 (s^c)^ (3.3-7) 


The vehicle motions model was specified in detail in 
Ref. 1. The equivalent state space representation of this 
effect is 

♦NASA is presently considering a longer computational cycle. 
If a change is made, the quantization models must be changed 
accordingly . 
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where 


North wind velocity 
West wind velocity 
West wind acceleration 


X 

— V 


North IMU displacement 
North IMU velocity 
West IMU displacement 
West IMU velocity 


(3.3-8) 


N is a white, gaussian random process and A, B, C and K are 
defined in Ref. 1. 

In the OFT launch configuration, the vertical IMU 
velocity is caused mainly by the vehicle sway motion. Assum- 
ing that this motion is primarily along the vehicle yaw direc- 
tion, then the vertical IMU velocity can be modeled as being 
proportional to the corresponding horizontal motion. The IMU 
velocity then becomes a measurement error source for the LSF 
filter since the filter uses successive velocity measurements 
to derive IMU accelerations and acceleration rates. 


3.4 PARTIAL ERROR BUDGET FOR CALIBRATION AND ALIGN!IENT 

The covariance analysis program for the IMU cal/align 
software was developed and checked out during this contract 
period; however, generation of a detailed error budget is part 


3-14 


76 



THE ANALYTIC SCIENCES CORPORATION 


of a planned future effort. During the program checkout 
process, error budget contributions for a number of spe- 
cific error groups were computed. The purpose for these 
runs was twofold: (1) to illustrate the type of results 

which can be obtained via covariance analysis, and (2) to 
improve understanding of the relative importance of certain 
major error sources. An example of the results obtained 
is presented in this section. 

Detailed error budgets showing the rms contribu- 
tion to all hangar and preflight calibrated parameters 
from two groups of error sources — second-order acceler- 
ometer nonlinearities and resolver harmonics — were genera- 
ted at the end of each cal/align sequence. Tables 3.4-1 
and 3.4-2 summarize the results from the hangar calibra- 
tion and the preflight calibration, respectively; where 
the former includes Hangar Coarse Alignment, Hangar Cal A 
and Hangar Cal B, and the latter Preflight Coarse Align- 
ment and Preflight Cal A. Each entry in these tables 
represents the rms contribution of the error source group 
indicated at the top to the calibrated parameters speci- 
fied in the left-hand column. As a comparison, the contri- 
bution f»*om all error source groups combined, as obtained 
from the monte carlo simulations reported in Ref. 1, are 
also shown. 

As is apparent from Tables 3.4-1 and 3.4-2, the 
second-order accelerometer nonlinearities are significant 
contributors to the calibration errors of accelerometer 
bias and scale factor asymmetries . This result was to be 
expected because the calibration algorithm cannot distin- 
guish accelerometer quadratic nonlinearities from biases 
and scale factor asymmetries. These results also confirm 
the conclusions obtained from the monte carlo simulations 
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TABLE 3.4-1 

HANGAR CALIBRATION ERROR BUDGET FOR 
ERROR SOURCE GROUPS 18 AND 23 


T-0405 


CALIBRATED PARAMETERS 

UNITS 

GROUP 16 

GROUP 23 

AU. GROUPS* 

Level Coarse Alignment 

sec 

0.5 

28.9 

216.0 

Azimuth Coarse Alig&fflent 

sec 

0.0 

28.9 

630.0 

Gyro Bias Drifts 

sec/ 8 

0.000 

0.000 

0.014 

Gyro Mass Unbalance 
Input Axis 

sec/s/g 

0.000 

0.000 

0.059 

Spin Axis 


0.000 

0.000 

0.020 

Output Axis 


0.000 

0.000 

0.007 

Gyro Torquer Scale Factors 

ppm 

8.1 

0.0 

64.0 

Gyro Misalignments 

sec 

0.0 

0.0 

8.0 

Accel. Biases-Low Gain 

vg 


0.0. 

33.0 

Accel. Scale Factor8->Lo« Gain 

ppm 

0.0 

0.0 

15.0 

Accel. Biases-Bigb Gain 

ug 


0.0 

43.0 

Accel. Scale Factors«High Gain 

ppm 

0.0 

0.0 

15.0 

Accel. Scale Factor Asymmetries 

ppm 


0.0 

59.0 

Accel . Konort hogonal i t ies 

sec 

0.0 

0.0 

2.0 

Accel, and Gyro Misalignments 

sec 

0.0 

0.0 

18.0 

Resolver Offsets, Gimbal 
Nonurthogonali ties 

S'*C 

0.0 


37,0 


•Based on Monte Carlo simulations including all error groups reported 
in Ref. 1 


TABLE 3.4-2 

PREFLIGKT CALIBRATION ERROR BUDGET FOR 
ERROR SOURCE GROUPS 18 AND 23 


T-?i9P ■ 


CALIBRATED PARAMETERS 

L*XITS 

GROUP 18 

CROUP 23 

ALL GROUPS* 

Level Coarse Alignment 

sec 

3.2 

29.1 

196.0 

Azimuth Coarse Alignment 

sec 

5.5 

37.5 

546.0 

Gyro Bias Drift 

sec/s 

0.000 

0.000 

0.009 

Accel. Biases-Low Gain 

ug 


0-0 

47.0 

Accel. Scale Factors-Low Gain 

ppm 

0.0 

0.0 

11. 0 

Accel. Biases-High Gain 

Ug 


0-0 

52.0 

Accel. Scale Factors»Klgh Gain 

ppm 

0.0 

1 

0.0 

13.0 


•Based on Uonte Carlo slsmlatloss using all error groups reported 
In Ref. 1 


P 
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which indicated that supplemental errors, primarily acceler- 
ometer heading sensitivities and nonlinearities, are signif- 
icant contributors. 

Table 3.4-1 shows that resolver harmonics are a 
significant contributor to the resolver offset and gimbal 
nonorthogonality calibration errors . The monte carlo 
sensitivity studies in Ref. 1 also support this conclusion. 
Tables 3.4-1 and 3.4-2 also indicate that resolver har- 
monics contribute to both Hangar and Preflight Coarse 
Alignment errors, however, these contributions are not 
significant when compared to other error sources. Another 
result exhibited by the error budgets is that the calibra- 
tion of the remaining IMU error parameters is insensitive 
to resolver harmonic errors . 

Each of the two error source groups causes larger 
coarse alignment errors during the preflight phase than 
during the hangar phase. This is so because both error 
source groups propagate into other accelerometer and 
resolver errors during hangar calibration. These errors, 
in turn, degrade the Preflight Coarse Alignment accuracy. 
This does not imply, however, that preflight coarse align- 
ment is less accurate than hangar coarse alignment. 
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4. ENTRY. PRELAND AND LANDING NAVIGATION 


Previous TASC reports (Refs. 1 through 3) have eval- 
uated entry, preland and landing navigation performance for a 
number of different hardware/ software configurations as 
part of NASA's development program for the Space Shuttle 
navigation system. As a consequence of performance limita- 
tions encountered with the previous filter design, NASA has 
modified the baseline filter and the recommended filter 
parameter values in an attempt to improve navigation accu- 
racy (Refs. 45 and 47 respectively). The objective of the 
current study is to generate performance projections for 
the revised filter for the OFT-1 mission in 1979. The level 
of detail provided in this analysis should be useful for 
inferring the error mechanisms which limit navigation sys- 
tem performance for the current mechanization and for 
evaluating the impact of proposed system modifications. 

The trajectory and measurement schedule for the per- 
formance study are presented in Section 4.1 and the baseline 
filter design is reviewed in Section 4.2. Section 4.3 con- 
tains the truth model which defines the environment in which 
the navigation performance is evaluated. The results of the 
navigation performance analysis and performance sensitivity 
to the dominant error sources are discussed in Sections 4.4 
and 4.5 respectively. 


4.1 TRAJECTORY AND MEASUREMENT SCHEDULE 

The reference trajectory for the entry, preland and 
landing phases is the nominal trajectory for the OFT-1 mission. 
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The OFT-1 mission entails a launch from Cape Kennedy into a 
38 degree inclination orbit and a landing at Edwards Air Force 
Base. The approach path is straight except for a dog-leg 
maneuver when the Terminal Area Energy Management (TAEM) align- 
ment cylinder is intercepted at an altitude of approximately 
28,000 ft. Figure 4.1-1 shows the ground track for the last 
segment of the entry phase and for the landing phase. The 
total elapsed time between entry interface at 400,000 ft and 
touchdown is 1630 sec. 





Figure 4.1-1 Approach and l,anding Ground Track 

for OFT-1 Mission 


The OFT-1 trajectory differs from previous trajec- 
tories used by TASC to evaluate the Shuttle entry and landing 
navigation systems in two important respects — it is low 
latitude rather than polar and the reentry flight path is 
much steeper. The low latitude approach will avoid the large 
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atmospheric density modeling errors associated with polar 
regions. Such modeling errors can significantly affect drag 
update accuracy. On the other hand, the steeper flight path 
angle will yield higher vehicle acceleration which, in turn, 
will make navigation system performance more sensitive to 
error sources such as IMU misalignments and accelerometer 
scale factor errors. 

The measurement schedule assumed for the OFT-1 mis- 
sion is indicated In Table 4.1-1. Processing of the first 
drag update measurement occurs at an altitude of 282,000 ft 
when the vehicle deceleration has reached approximately 0.025g. 
Acquisition of TACAN at 145,000 ft ends the entry navigation 
phase and signals the start of preland navigation. Drag up- 
date measurements are processed down to 100,000 ft; the baro 
altimeter is switched on at an altit-..ue of 85,000 ft. MLS 
acquisition occurs at 20,000 ft wLen the Shuttle is on the 
TAEM alignment cylinder and the range-to-g .s approximately 
51,000 ft. Through the entry and preland navigation phases 
the measurements are processed at 3.84 sec interv'als. 

The switch to the landing navigation phase (t=Tj^^.) 
occurs at 13,580 ft (approximately 12,000 ft above the Shuttle 
touchdown point), at which time the update interval is reduced 
to 1.92 sec. Finally, at an altitude of approximately 100 ft 
above the runway (near the runway threshold), the MLS eleva- 
tion measurement is replaced by radar altimeter measurements 
at an update interval of 0.16 sec. 

The antenna locations for the radio navigation aids 
are given in Table 4.1-2. The TACAN antenna is located 
34,320 ft short of the nominal touchdown point on the runway 
center line. Lo< ation of the TACAN antenna on the Shuttle 
approach path rather than at the far end of tho runway (as 
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TABLE 4.1-1 

KEY EVENTS FOR OFT-1 MISSION 


XAVICATIOX 

PHASE 


E>THY r 

IXAVICATIO.V i 
PHASE V 


PREULVD 
IXAVICATION X 
PHASE 


, LAKDIKC I 
I NAVIGATION i 
PHASE 


T-0489 


ZVZST 

ELAPSED 

TIME 

(sec) 

1 TIME TO 
1 TOUCHDOWN 
(sec) 

ALTirJDE 
1 ABOVE GEOID 

RELATIVE 

VELOCITY 

(fps) 

HANCE-TO- 

GO 

(it) 

Entry Interlace. 

0-0 

1630.1 

398,837. 

24295. 

19.600.000. 

First Drag Measuree^at . 

188.2 

1441-9 

282.174. 

24399. 

15.400.000. 

First TACA.K Measurenent, 

1052.2 

577.9 

144.780. 

7391. 

1,425.000. 

Last Drag Ifeasureaent 

1248 . 0 

382.1 

100,263. 

3446. 

479.000. 

First Bare Altimeter Measurement. 
^BA 

1309.4 

320.7 

84.257, 

2447. 

282.000. 

First MLS Aziautb. Elevation, & DUE 
Measurements . 

(No TaCAN Measurement) 

(No Bare Altimeter) 

1536 . 0 

94.1 

20.370. 

662. 

51.000. 

Reduce Filter Dimen-ion. 

Switch from 3.84 sec Update Cycle 
to 1.92 sec Update Cycle 

1562.9 

67.2 

13,580. 

622. 

35,000. 

First Radar Altimeter Measurement. 
^RA Elevation Xeasurtsment ) 

Switch from 1.92 sec Update Cycle 
to 0-16 sec Update Cycle 

1622.4 

7.7 

2,387. 

362. 

2.700. 

Touchdown 

1630.1 

c.o 

2.293. 

315. 

0. 


TABLE 4.1-2 

SENSOR LOCATIONS WITH RESPECT TO TOUCHDOWN POINT 


SENSOR 

DOWNRANGE 

(ft) 

CROSS RANGE 
(ft) 

TACAN 

34,320 

0 

MLS Elevation 

850 

300 

MLS Azimuth 

13,800 

300 

and DM£ 




had been assumed in Ref. 3) should yield improved geometry 
during the TAEM alignment maneuver. The MLS antennas are 
offset 300 ft from the cent r line, approximately where they 
had been located in previous TASC studies. 
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4.2 NAVIGATION FILTER 


The navigation filter for the entry through landing 
phases is the baseline filter specified in Ref. 45*. Before 


SW 


the filter is a variable dimension Kalman filter with up 


to 12 states. After Tg^ it is a 6 state complementary filter. 
A schedule of filter states is presented in Fig. 4.2-1. 


Implementation of the filter covariance equations 
is detailed in Refs. 1 and 3. The data base for the filter 
includes an initial covariance matrix and parameter values 
associated with each of the navigation sensors, including 
both the correlated error states and the measurement noise. 
The initial covariance matrix at initiation of the deorbit 
burn was provided by NASA (Ref. 47) and assumes that the 
most recent navigation fix and IB4U alignment were two hours 


n-ii305b 


AF^flOXtiiATe 

ALTiTuOC 

AfrOVE RkJMMAV 


38 Gl 000 rT 

I 


XWXMfT 


XA&JtBD Ft 
t 


3 QJ 0 CCFT 

I 


12JOOO FT 1W FT TOUCHOCWM 

I I I 


1 


FosiTiOA AfiO vtLOcrry 


! 


MtSALlGJ«M£»{T$ 


FlLTfA 

STATES 


10 


OAAC UFOaTE 


I- 


•AAO 

ALTumCTEK 


I 


MLS ELEVaTHM 


It 


TACAM OmE 


MLS OmE 


TACAM VOA 


«kS A 21 MUTH 


Figure 4.2-1 Schedule for Entry, Preland and Landi^r 

Phase Navigation Filler States 


♦The navigation filter has undergone a major redesign since 
this effort was completed. Navigation performance with the 
new filter will be evaluated in a subsequent report. 
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prior to deorbit. The covariance matrix is diagonal except for 
correlations between downrange and crossrange errors . The 
remaining parameter values were also provided by NASA. The 
filter data base is summarized in Table 4.2-1. 


Prior to the landing navigation phase the measure- 
ment matrices for all external navigation aids are of the form 

= jh 1 0 ! 0 1 h.l (4.2-1) 

F {_ p I I I bj 

where h is the 1x3 partial derivative of the measurement 
P 

with respect to the nominal position and h^ is the partial 


TABLE 4.2-1 

FILTER STATE AND MEASUREMENT ERROR STATISTICS FOR 
ENTRY AND PREL4ND NAVIGATION FILTER 


T-0267 



Correlated Error 

Filter State 

Measurement 

Noise 

Heasur^ent 

Standard 

Deviation 

Correlation 

Time 

(sec) 

Standard 

Deviation 

Position 

16,700 ft 

oo 

— 

Velocity 

18.7 fps 

CO 

— 

Misaligirnents 

4.4 mrad 

00 

— 

Drag-Altitude 

9000 ft 

120 

5148 ft 

Baro Altimeter 

443 ft + 3.5% 
of altitude 

60 

60-750 ft (altitude 
dependent) 

TACAK DME 

524 ft 

80 

131 ft 

TACAN VOR 

29 mrad 

20 

14 mrad 

MLS DME 

80 ft 

. 1 

1 

26 ft 

ULS Azimuth 

0.96 mrad 

00 

0.3 mrad 

MLS Elevation 

0.96 mrad 

1 

0.3 mrad 


♦Downrange and radial navigation errors will be coupled due to 
the effect of orbital dynamics — a downrange velocity error 
at orbital perigee will generate a radial position error at 
orbital apogee, etc. The filter accounts for this effect by 
assuming -1 correlations between downrange velocity and ver- 
tical position and between downrange position and vertical 
velocity . 
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derivative with respect to the bias states. The drag update 

measurement is implemented as a pseudo-altitude measurement . 

During landing navigation H„ is of the form; 

r 


H 


F 




h I 
P I 



(4.2-2) 


For all external navigation aids, the first three elements 
of Hp are identical to those of Hjg in the truth model. 


The gain computation algorithm in the Kalman filter 
mode (i.e., prior to Tg^) includes measurement underweighting 
in order to minimize the possibility of filter divergence: 

1-1 


K = PH' 


u.^qHPh'^ + R 


(4.2-3) 


where the measurement underweighting factor u.^q is 1.2 when 
the rss (filter- indicated) position error is greater than 
3280 ft and 1.0 otherwise. For the measurement schedule and 
filter parameter values used in this study, underweighting 
is switched off (u^q is set to unity) at approximately 
49,000 ft. 


In the complementary filter mode (i.e., after Tg^) 
the filter covariance is not computed and all gains are based 
upon pre- stored gain schedules. In this mode the >!LS DME and 
radar altimeter gains are padloaded constants scaled by the 
upaate cycle, anJ the MLS azimuth and elevation gains are 
padloaded constants scaled by the update cycle and the range- 
to-go. The gains used for the study are summarized in Table 
4.2-2 and are based on Ref. 45. 


*In previous efforts the drag update measurement was imple- 
mented as a pseudo-drag measurement (Rets. 1 and 3). The 
measurement matrix Hp in that case included a partial deriva- 
tive with respect to the nominal velocity, h , which was non- 
zero . 
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TABLE 4.2-2 

LANDING NAVIGATION GAINS IN RUNWAY COORDINATES 


T-M93 


MEASUBEUDH 

POSITION gains’ 

VELOCITY GAINS* 

R 

PR 

CR 

R 

DR 

CR 

ULS Azimuth 

0 

0 

0.2 R^6 

0 

0 

0.02 Rgoi 

MLS Range 

• 0 

0.26 

0 

0 

0.026 

0 

IILS Elevation 

0.1 Bgo6 

0 

0 

0.01 Rgo6 

0 

0 

Radar Altimeter 

? .56256 

0 

0 

0 

0 

0 


•Bqq Is raaee-t'v.go. 6 is measureiMnt ictervsl (1.92 s«c for HLS 
measuremeDte, 0.16 sec for the radar altimeter ffleaeureiaent ) . 

The filter outlined above has been exercised over the 
trajectory and measurement schedule presented in Section 4.1. 
The time history of the filter-indicated performance between 
the deorbit burn and initiation of the landing phase is sum- 
marized in Table 4.2-3 in the form of rms values of position 
and velocity errors in the UVW coordinate frame (vertical, 
downrange, crossrange). The filter does not compute a per- 
formance measure (covariance) during the landing phase. 

TABLE 4.2-3 

FILTER-INDICATED PERFORMANCE 


T-04M 


EVENT 

ELAPSED 

TIME 

RMS POSITION ERROR 
(ft) 

RXS VEUXriTY 
(fps) 

ERROR 


(sec) 

R 

DR 

CR 

R 

DR 

CR 

iDii.atioQ of Deorbit 
Burn 

— 

16700 

16700 

16700 

18.7 

18.7 

18.7 

Entry Interface 

0.0 

7016 

35391 

15942 

28.6 

6.7 

19.7 


163.2 

4633 

36917 

15828 

29.9 

3.8 

19.9 

First Drag lfe&.suremeot 

188.2* 

3966 

35639 

15813 

29.2 

3.1 

19-9 


1029.1* 

7959 

27189 

32753 

35.9 

34.9 

79.8 

First Tacan Measurements 

1052.2* 

7847 

12575 

28088 

25.5 

26.0 

73. 


1305.6* 

6693 * 

2275 

4163 

22.6 

14.5 

30.6 

First Baro Altimeter 
Measurement 

1309.4* 

3717 

1553 

4047 

18.2 

11.5 

29-6 


1536. 0" 

828 

964 

459 

10.7 

6.4 

5.6 

Ids Acqulsltior 

1536. 0"*“ 

63 

69 

74 

9.7 

6.8 

4.6 

Initiation of Landing 
Navigation Phase 

1562.9* 

1 

54 

78 

44 

1.0 

1.2 

i 

0.9 
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The filter- indicated performance down to the first 
TACAN measurement at 1052.2 sec is considerably worse than 
that determined in Ref. 3, primarily because the initial I!1U 
misalignment and crossrange navigation uncertainties are much 
more severe than in the previous study. This pessimism rela- 
tive to the previous study continues down to MLS acquisition 
because of larger uncertainties assumed for the TACAN and 
baro altimeter measurements. At the switch to landing naviga- 
tion, the filter-indicated errors are approximately twice as 
large as those in Ref. 3. 

4.3 ENTRY, PRELAND AND LANDING PHASE TRUTH MODEL 

The "real world" environment in which the filter is 
evaluated is characterized by the truth model . i.e., a mathe- 
matical model of the various error sources which occur in 
the operational environment for the navigation system. The 
truth model for the entry phase is an extended version of the 
model first presented in Ref. 3. The truth model for the 
preland and landing phases is a modification of that detailed 
in Ref. 1. 

4.3.1 Truth Model Definitions 

The truth model for the entry phase consists of a 
total of 47 error states plus random errors associated with 
accelerometer quantization and the altitude pseudo-measure- 
ment. The error states are listed in Table 4.3-1 which gioups 
them into four general categories; 

• Filter states and measurement noises 

• IMU-related states 

• Gravity model errors 

• Drag-related states 
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TABLE 4.3-1 

ENTRY PHASE TRUTH MODEL STATES AND ERROR SOURCES 


T-0495 



NUMBER 

NUIBER 

ERROR SOURCE 

OF 

OF ERROR 


STATES 

SOURCES 

I. FILTER STATES AND lIEASUREflEKT NOISES 



1. Position Errors 

3 

3 

Velocity Errors 

3 

3 

Platform Misalignments 

3 

3 

Correlated Measurement Errors 



Altitude 

1 

1 

Uncorrelated Measurement Noise 



Altitude 

— 

1 

INS Quantization Noise 

— 

3 

II. IMU ERRORS 



Accelerometer 



2. Biases 

3 

3 

3. Scale Factors 

3 

3 

4. Asymmetries 

3 

3 

5. Non-orthogonalit ies 

3 

3 

Gyro 



7. Bias Drifts 

3 

3 

8. Mass Unbalances 

5 

5 

III. GRAVITY MODEL ERRORS 



6. Gravitational Deflections 

3 

3 

and Anomalies 



IV. NON^ESTI MATED. DRAG- RELATED STATES 



21. Non-Standard Density 



1962 Standard Atmosphere Error 

1 

3 

Time-Varying Bias 

1 

1 

First-Order Markov 

1 

1 

22. Non-Standard Wind 



Westerly (Time-Varying Bias) 

1 

1 

Headwind (First-Order Markov) 

1 

1 

(''•Dssw'ind (Second-Order Markov) 

2 

2 

Turbulence (First- and Second- 

5 

5 

Order Markovs) 



23. Non-Standard Aerodynamics 



Time-Varying Bias 

1 

1 

First-Order Markov 

1 

1 

TOTALS 



51 
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The principal differences between the entry phase truth model 
used in this study and that of Ref. 3 are improvements in the 
,^IMU error model and in the aerodynamics error model in the drag- 
related states. These changes are discussed in detail in 
Section 4.3.2 where the truth model data base is presented. 

The truth model states and other error sources used 
to evaluate the Shuttle navigation system during the preland 
and landing phases are listed in Table 4.3-2, which also 
divides them into four categories: 

• Filter states and measurement noises 

• IMU-related states 

• Gravity model errors 

• External aid-related states 

The last category includes the drag-related states of Table 
4.3-1. In addition to the changes in the IMU and aerodynamics 
models cited above, the major modifications to the truth 
model for these phases involve a revision of the data base to 
reflect the best a\'ailable sensor performance data, and the 
addition of radar altimeter terrain modeling errors for 
several additional runways. 

4.3.2 Truth Model Data Base 

This section presents the truth model data base for 
the entry through landing phases of the OFT-1 mission. A 
short discussion of the change to the aerodynamics error 
model is included. 

Initial conditions at deorbit will be a function 
of the accuracy of the most recent navigation system update 
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TABLE 4.3-2 

TRUTH MODEL STATES AND ERROR SOURCES FOR THE 
PRELAND AND LANDING PHASES 


T-0496 



NUMBER 

NUMBER 

ERROR SOURCE 

OF 

OF ERROR 


STATES 

SOURCES 

I. FILTER STATES AND MEASUREMENT NOISES 



1. Position Errors 

3 

3 

Velocity Errors 

3 

3 

Platform Misaligoments 

3 

3 

Correlated Sensor Measure- 

7 

7 

ment Errors 



Uncorrelated Measurement Noise 



Navigation Sensors 

— 

8 

INS Quantization Noise 

— 

3 

II. IlfU ERRORS 



Accelerometer 



2. Biases 

3 

3 

3. Scale Factors 

3 

3 

4. Asynroetries 

3 

3 

5. NoD-orthogonalitles 

3 

3 

Gyro 



7. Bias Drifts 

3 

3 

6. Hass Unbalances 

5 

5 

III. GRAVITY MODEL ERRORS 



6. Gravitational Deflections 

3 

3 

and Anomalies 



IV. NON-ESTIMATED. EXTERNAL ATn-KFI.ATED 



STATES 



10. TACAK Range Bias First-Order 

1 

1 

Markox' 



11. TACAN Range Scale Factor 

1 

1 

12. Baro Altimeter Errors 



Bias 

1 

1 

Scale lactor 

1 

1 

First-Order Markov 

1 

1 

Static Defect 

1 

1 

13. UL5 Scal€> Factor and Biases 



Range Bias 

] 

1 

Range Scale Factor 

1 

1 

Azimuth 

1 

1 

Elevation 

1 

1 

14. MLS Second-Order Markov 



Range 

2 

1 

Azimuth 

2 

1 

Elevat ion 

2 

1 

3. Radar Altimeter Errors 



Bias 

1 

1 

Terrain Modeling 

1 

1 

Terrain Averaging 

1 

1 

Antenna Poi* ting 

1 

1 

Timing Delay 

1 

1 

Receiver Noise 

1 

1 

Fading 

1 

1 

16. TACAN Survey Errors 

6 

6 

17. MLS Surx Errors 

6 

6 

19. TAC *■ B(“aring Bias First-Order 

1 

1 

Markov 



20. MLS Tlrtinr Bias 


3 

21 .22 .2'^ Drug-Related .ales 

4 

14 

(Tai 14 3-1) 



TOTAI-S 

92 

t 

100 
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and of the subsequent open loop propagation. The initial co- 
variance for the current study assumes an uploaded navigation 
fix and an IMU alignment via star tracker two hours prior to 
deorbit (Ref. 47). The resulting navigation errors at deorbit 
are primarily in-plane errors - 12,000 ft and 13 fps rss. 

The IMU misalignments are 285 sec rms per axis, which accounts 
for both the star tracker errors and gyro drift during the 
two hour interval. 

The IMU error model has been modified to be con- 
sistent with the IMU performance specification (Ref. 28). 

The IMU states are a subset of those defined for the IMU 
calibration and alignment analysis in Ref. 1. Gyro aniso- 
elasticities , accelerometer nonlinearities, etc. have been 
ignored in this study because of their negligible effect on 
navigation errors (Refs. 1 and 3). Table 4.3-3 presents the 
truth model data base for the IMU-related error sources. 

The major change to the drag-related error models 
is an improvement in the aerodynamics error model to reflect 
the relationship between drag coefficient errors and Shuttle 


TABLE 4.3-3 

TRUTH MODEL DATA BASE 
FOR B!U- RELATED ERROR SOURCES 


T-0499 


ERROR SOURCE 

STA!;DARD DEVIATION 

DATA* 

SOURCE 

1. 

INS Quantization Error 

1 , 0 cm/ sec 

Ref 

. 28 

2. 

Accelerometer Biases 

50 ug 

Ref 

. 28 

3. 

Accelerometer Scale Factors 

100 ppm 

Ref 

. 28 

4. 

Accelerometer Asymmetries 

100 ppm 

Ref 

. 28 

5. 

Accelerometer Nonorthogonalities 

15 sec 

Ref 

. 28 

t . 

Gyro Bias Drifts 

0.035 deg/hr 

Ref 

. 28 

8. 

Gyro Mass Unbalances 

0. 025 deg/hr/g 

Ref 

. 28 
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control surface deflections. The variation (Group 23), 
previously modeled as a stationary markov process, has been 
modeled as the markov process plus a time-varying bias term. 
This bias term is determined as a function of relative velo- 
city, flight-path angle and control surface deflections. 

The drag coefficient computational model used in the 
filter is a quadratic function of the computed angle of attack 

- <=0 * =l“c * Vc 
c 

where Cq, and C 2 are constants. Aerodynamic data presented 
in Ref. 46 gives the true drag coefficient by: 


6gp, 6gg, 5j.) 


Cg (a,M, 6gp, 5gg, 6g) + ACg (a,M,5gp) 
untr BF 


+ ACg (a,M,6gg) + ACg (a,M,6p) 
SB B 


(4.3-2) 


where 


cx(t) is the angle of attack, 

M(t) is the Mach number, and 

Cg is the untrimmed drag coefficient, 

untr 

AC-. , AC^ and AC-. are the increments of the drag 

coefficient due to the body flap, speed 
brake and eleven settings, respectively, 

5gp(t), 5gg(t) and 0 £(t) are the body flap, speed 

brake and eleven deflections, respectively. 

With the assumption that the error in ot^(t) is a negligible 
contributor to the error in the deviation of from 

the true Cp, for the OFT-1 mission is presented in Fig. 4.3-1. 
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Figure 4.3-1 Deviation of Navigation Filter C„ 

Model From True 

The greater deviations observed at high altitudes are a func- 
tion of the poorer accuracy of Cj)^ for the larger angles of 
attack used in the upper atmosphere. The range of this model- 
ing error is 1.1% to -9.5%. 

The Group 21 time-varying density bias which exhibits 
both latitude- and season-dependence, is shown in Figs. 4.3-2a 
and 4.3-2b respectively as percent deviations from the 1962 
Standard Atmosphere. Dashed lines on the profiles denote the 
mean density deviations corresponding to the altitude-latitude 
profile for the OFT-1 mission. Because of the low inclination 
for OFT-1, this error is smaller at high altitudes than it 
would be for higher inclination orbits such as reference 
mission 3B (see Ref. 3). The remaining drag-related errors 
utilize the same data base presented in Ref. 3. Table 4.3-4 
summarizes the parameter ranges for the drag-related error 
models (Group 1, Groups 21 through 23). 
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Figure 4.3-2a Time-Varying Density Bias for December 

January Given as Percent Departure From 
1962 Standard Atmosphere (Ref, 53) 



05PA«TUP€ CPOM t9«2 STANOAPO OCNSiTv 


Figure 4.3-2b Time- Varying Density Bias for June - July 

Given as Percent Departure From 1962 
Standard Atmosphere (Ref. 53) 
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TABLE 4 . 3-4 

TRUTH MODEL DATA BASE 
FOR DRAG- RELATED ERROR SOURCES 


T-0497 


ERROR SOURCE 

STANDARD DEVIATION 

C0RB£LATI0X TltO: 
(Distance ) 

DATA SOURCE 

1. 

Uncorrelated Altitude 
Pseudo-Measurement * Noise 

1250 tt 


Ref. 

3 

21. 

Kon-Standard Density 






1962 Standard Atmosphere 

0 - 21% 

— 

Ref. 

3 


Error 

Time-Varying Bias 

0 - m 


Fig. .i 

.3-2 


First-Order Markov 

2-6% 

(14 - 270 nm) 

Ref. 

3 

22. 

Non-Standard Wind 






Westerly (Time-Varying 

0 - 295 fps 

— 

Ref. 

3 


Bias) 

Headwind (First-Order 

(270 cm) 

Ref. 

3 

40 - 110 fps 


Markov ) 

Crosswind (Second-Order 

(270 nm) 

Ref. 

3 

40 - 110 fps 


Markov ) 





Turbulence (First- and 

13 fps 

(1750 ft) 

Ref. 

3 


Second-Order Markovs) 





23. 

Cp Variations 






Time-Varying Bias 

0-9.5% 


Fig. 4 

.3-1 


First-Order Markov 

5% 

100 sec j 

Ref. 

3 


A detailed discussion of the radar altimeter error 
model can be found in Ref. 1. Terrain modeling errors for 
runways at Andersen Air Force Base, Guam, and Hickam Air 
Force Bata, Hawaii, have been added to the radar altimeter- 
related truth model data base. Both of these locations are 
secondary landing sites for the Space Shuttle. Figure 4.3-3 
shows the terrain modeling errors for both the north and 
south runways at Andersen, and for runways 8 and 26 at Kickam. 

The truth model data base for the external aid- 
related error sources of Group 1 and Groups 10 through 20 
is presented in Table 4.3-5. Although the state definitions 
are the same as developed in Refs. 1 and 3, many of the data 
base v'alues have been revised to reflect more recent sensor 
performance data. The table contains applicable references 
for the data values selected. 

The gravity model errors (Group 6). which have a 
negligible effect on system performance, have been ignored 
in this study. 


9(o 
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R-mS4o 


AAFB- NORTH RUNWAY 
(Northeast Approach) 


AAFB- SOUTH RUNWAY 
(Northeast Approach) 


HAFB-26 


HAFB-6 


a 


A 



RANGE TO- GO (FT) 


28 26 24 22 20 18 16 14 12 10 8 6 4 2 

TIME-TO^O (SEC) J f 


Figure 4.3-3 Terrain Modeling Error for Approach 

^aths to Secondary Space Shuttle 
Landing Sites 
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TABLE 4.3-5 

TRUTH HODEL DATA BASE FOR EXTERNAL 
AID-RELATED ERROR SOURCES 


T-0498 


ERROR SOURCE 

STANDARD 

DEVIATIW 

TIME 

CONSTANT 

DATA 

SOURCE 

1. 

UDCorrelated Ueasuremeot Noise 





TACAK Range 

100 ft 

— 

Ref. 49 


TACAN Bearing 

12 mrad 

— 

Refs. 49.35 


Baro Altimeter 

0.24 ft 


Ref. 33 


ULS Azimuth 

0.2 mrad 


Ref. 34 


MLS Elevation 

0.1 mrad 

— 

Ref. 34 


MLS Range 

40 ft 

— 

Refs. 34.35 


Radar Altimeter 

0.9 ft 

— 

Ref. 1 

10. 

TACAN Range Bias First<-Order 

385 ft 

300 sec 

Ref. 49 


Markov 




11. 

TACAN Range Scale Factor 

100 ppm 

m 

Ref. 49 

12. 

Baro Altimeter Errors 





Bias 

100 ft 

m 

Ref. 2 


Scale Factor 

3% of alt 

m 

Ref. 2 


Static Defect 

1.52 jc 10“^ 
ft/ft^/sec^ 

m 

Ref. 2 


First-Order Markov 

20 ft 

100 sec 

Ref. 2 

13. 

MLS Scale Factor and Biases 





Range Bias 

60 ft 

«D 

Ref. 34 


Range Scale Factor 

400 ppm 

m 

Ref. 34 


Azimuth 

0.6 mrad 

ea 

Ref. 34 


Elevation 

0.4 mrad 

«R 

Ref. 34 

14. 

MLS Second -Order Markovs 





Range 

40 ft 

4 sec 

Refs. 34,35 


Azimuth 

0.3 mrad 

4 sec 

Ref. 34 


Elevation 

0.3 mrad 

4 sec 

Ref. 34 

15. 

Radar Altimeter Errors 


1 



Bias 

1 ft 

__ 



Terrain Modeling 

< 150 ft* 


Ref. 1 


Terrain Averaging 

< 25 ft* 


Pef. 1 


Antenna Point ing 

< 7 ft* 


Ref. 1 


Timing Delay 

< 0.5 ft* ! 


Ref. 1 


Receiver Noise 

negligible 

0.5 sec 

Ref. 1 


Fading 

< 0.3 ft* 

0.5 sec 

Ref. 1 

16. 

TACAN Survey Errors 

1 ft 

00 

- 

17. 

MLS Survey Errors 

1 ft 

CO 

- 

19. 

TACAN Bearing Bias First-Order 

12 mrad 

2 X 10^ 

Refs. 49,35 


Harkov 




20. 

MLS Timing Bias 

50 msec 

OP 

- 


^Absolute range rather than standard deviation. 


<76 
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4.4 ENTRY. PRELAND AND LANDING NAVIGATION PERFORMANCE 

This section presents detailed performance results for 
the Space Shuttle entry, preland and landing navigation system. 
Overall performance curves, showing the time histories of 
position and velocty errors due to all error sources combined, 
are given. "Baseline Error Budget" tables show the contri- 
butions of individual error sources, or small groups of error 
sources, at specific times, to position and velocity errors. 
These detailed error budgets have been computed at five signi- 
ficant time points: 

• The end of the entry navigation phase 

(t » T,j.^cAN’ * 145,000 ft) 

• The initiation of MLS measurement processing 

(t * Alt. = 20,000 ft) 

• The initiation of the landing navigation phase 
(t = *^SW’ Alt.® 12,000 ft above the runway) 

• The initiation of radar altimeter measurement 
processing (t = T„^, Alt. = 100 ft above the 
runway ) 

• Touchdown (t * Tjj)- 

Section 4.4.1 presents the results from entry down to 20,000 ft 
and Section 4.4.2 presents the results for altitudes below 
20,000 ft. 

4.4.1 Navigation System Performance: Entry to 20,000 ft 

Figures 4.4-la and 4.4.1b present overall perfor- 
mance curves for the entry phase and the initial portion of 
the preland phase of the OFT-1 mission. The curves show rms 
position and velocity errors due to the combined effects of 
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14000 



«.24C3e 



Figure 4.4-lb Overall Performance Down 

to 20,000 ft: Velocity 


* 
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REPRODUCIBILITY OF TH£ 
ORIGINAL PAGE IS POOR 


all error source groups in the truth model. They were gen- 
erated by root-sum-squaring individual error source contri- 
butions at specified intervals (as tabulated in Appendix C). 

The initial covariance for the entry filter is quite 
pessimistic relative to the truth model errors and this pessi- 
mism is reflected in the overall performance curves. The 
filter-indicated downrange error at is 35,600 ft rms vs 

11,400 ft rms for the truth model. As a consequence, the 
filter overweights the first drag update measurement to 
yield a net increase in the downrange position error. A simi- 
lar pessimism pervades crossrange position uncertainty and all 
three components of velocity, but the altitude pseudo-measure- 
ment does not contain sufficient information to attempt to 
improve these estimates. 


After the initial transients, the drag update filter 
does an effective job of limiting downrange and radial posi- 
tion errors. The growth in the velocity errors (and in the 
crossrange position error) is attributable primarily to the 
large IMU misalignments, although the errors are not nearly 
as large as anticipated by the filter. The increase in the 


velocity errors at T, 


TACAN 


is disconcerting, but again it is 


attributable to the pessimistic filter-indicated performance . 
(At the filter-indicated rms crossrange velocity 

error is 80 fps vs a truth model uncertainty of 14 fps. ) 
Following the navigation errors exhibit a general 

decreasing trend. At Tg^ vertical position and velocity 
errors show sharp declines, reflecting the acquisition of the 
baro altimeter measurement. The processing of TACAN bearing 
measurements accounts for the large steady decrease in 
crossrange errors between and 


To illustrate the effect that the drag update mea- 
surements have on the vertical position error, several truth 


/b/ 
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model runs without this "external aid" were made. Just prior 
to TACAN acquisition the total vertical position error for this 
no measurement case is approximately 4700 ft as compared to the 
1436 ft observed when the drag measurements are employed. 
Clearly then, the processing of the drag altitude pseudo-mea- 
surements distinctly reduces vertical position errors . 

System error covariances have been computed for each 
group of error sources down to the 20,000 ft altitude of the 
reference trajectory. Detailed error budgets showing rms 
position and velocity errors were generated at the end of the 
entry phase (t * '^tACAN’ “ 145,000 ft) and at 20,000 ft. 

Each entry in the error budgets is the rms contribution of 
the error source or sources indicated in the left-hand column. 
The entries in each coliimn are root -sum-squared to obtain the 
total projected performance of the navigation filter. These 
performance values are compared with the Space Shuttle land- 
ing specifications (from Ref. 43) at the bottom of the error 

* 

budget at 20,000 ft , To focus attention to the major contri- 
butors in each colvimn, every number whose magnitude is greater 
than 20% of the rms total for that column has been circled. 

Table 4.4-1 presents the baseline error budget for 
the entry filter. The major contributors to navigation errors 
at 145,000 ft are: 

• The initial condition errors (Group 1) 

• The altitude pseudo-measurement noise 
(Group 1) 

• The non-standard atmospheric density 
modeling error (Group 21) 

These landing specifications are guidelines, not inviolable 
requirements. The vertical velocity specification, in 
particular, may ultimately be relaxed. 
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TABLE 4.4-1 

BASELINE ERROR BUDGET FOR THE ENTRY NAVIGATION PHASE 


T»0S01 






RMS 

NAVIGATION ERROh AT 14S.000 ft | 



ERROR SOURCE 

VALUE 

POSITION (ft) 

VELOCITY (fps) 1 





R 

DR 

CR 

R 

DR 

CR 

I. 

FILTER STATES AND MEASUREHENT 









NOISE 










Initial Condition Errors 









Positions and Velocities 

Ref. 47 
285 sec 

153. C10490,)C3618.1I 

fl6.S7) 


2.00 


Uisalignments 

^357) 

mr 


C5.63) 

C5.03)C13.501I 


Uncorrelated Measurement Noise 








Altitude 

1250 ft 

( 47d.) 

697. 

306. 

2.01 

0.70 

0.36 


INS Quantization Noise 

0.0328 fps 

27. 

160. 

132. 

0.55 

0.47 

0.13 

II. 

IMD-BELATED STATES 









Accelerometer Errors 





(5.31) 

0.53 

0.41 

0.42 

1.32 

0.12 




2. Biases 

3. Scale Factors 

4. Asynmetries 

5. Nonortbogonalities 
Gyro Errors 

7. Bias Drifts 

8. Mass Unbalances 

50 ng 
100 ppm 
100 
IS ^ 

0.035®/hr 

0.025®/hr/g 

r"32cm 

6. 

45. 

31. 

105. 

6. 

1906. 

161. 

206. 

58. 

257. 

23. 

C1633.) 

53. 

80. 

19. 

428. 

25. 

(EUD 

0,87 

0.95 

0.04 

1.29 

0.14 

1.22 

0.29 

0.38 

0.02 

2.51 

0.20 

IV. 

DRAG- 

RELATED STATES 









21. 

Non-Standard Density 










1962 Standard Atmosphere 
Modeling Error 

4-Term (<9 %) 
Ref. 3 

(515E)(5m7) 

(5^ 

(5.56) 

1.25 



Time-Varying Bias 

Winter (<8 5f) 
Fig. 4.3-2a 

142. 

117. 

248. 

0.52 

0.43 

0.46 



First-Order Markov 

(<6%) Ref. 3 

177. 

500. 

303. 

0.72 

0.55 

0.33 


22. 

Non-Standard Wind 
Westerly 

Winter 
(<295 fps) 
Ref. 3 

j 

259. 

9. 

4. 

0.72 

0.20 

0.04 



Crosswind | 
Beadwlnd | 

(<110 fps) 
Ref. 3 

84. 

89. 

55. 

0.28 

0.12 

0.06 


23. 

Non-Standard 

Aerodynamics 










Time-Varying Bias 

(<9.6%) 
Fig. 4.3-1 

(865.) 

67. 

796. 

2.37 

1.76 

2.41 



First-Ordei Markov 

5% of 

208. 

438. 

278. 

0.77 

0.51 

0.31 

Total Projected Performance 


1436. 

13805. 

5516. 

19.76 

30.06 

14.21 


Except for vertical position and crossrange velocity, which 
are dominated by drag-related errors and initial misalignments 
respectively, each of the navigation errors is attribut- 
able to both initial condition and drag-related errors. 

This relative balance between initial conditions and drag- 
related errors suggests that the filter combines the altitude 
pseudo-measurements with inertial information in an effective 
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manner . However, filter parameter and gain variation studies 
would be required to prove this hypothesis. 

The large error contributions of initial position 
and velocity uncertainties can be explained via a mechanism 
which couples these initial errors to the navigation errors 
at 145,000 ft (and later) through platform misalignments. 

The filter responds to the large initial position and veloc- 
ity uncertainties by adjusting the misalignment estimates 
to account for a portion of these errors - thus, in effect, 
introducing an initial position- and velocity-generated 
misalignment. (This behavior can be observed in the appropri- 
ate table — Group 1 — in Appendix C by noting the immediate 
increase in platform tilt errors following the first drag 
update at f * ) In turn, the misalignments propagate 

via the system dynamics directly into position and velocity 
errors . 


Of the IHU-related sources, only the accelerometer 
bias errors contribute significantly to vertical and down- 
range velocity errors. The minor role of IMU error sources 
in this study vis-a-vis Ref. 3 can be attributed to a differ- 
ence in assumed initial conditions. In Ref. 3 the navigation 
system had been operating in the pure inertial mode since 
launch. Hence, the majority of the navigation errors at de- 
orbit were directly attributable to IMU error sources. By 
contrast, the present study assumes a navigation fix and IMU 
alignment two hours prior to deorbit. This eliminates all 
IMU contributions to deorbit navigation errors except those 
which can accrue (accelerometer bias, gyro bias drift) dur- 
ing the ensuing two hour orbital coast. 

Table 4.4-2 lists the contributions of alternative 
error source models (e.g., season-dependent atmospheric den- 
sity error models for the summer instead of the winter profile). 
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TABLE 4.4-2 

ALTERNATIVE CONTRIBUTIONS FOR THE ENTRY NAVIGATION PHASE 

T-0502 





RMS 

NAVIGATION ERROR AT 145.000 

ft 


ERROR SOURCE 

VALUE 

POSITION 

(ft) 

VELOCITY (fps) 




R 

DR 

CR 

R 

DR 

CR 

21. 

Non-Standard Density 









1962 Standard Atmosphere 
Uodellng Error 

R«f. 3 

3813. 

2830. 

3076. 

7.54 

2.00 

5.26 


Time* Varying Bias 

Sumer ( <16%) 
Fig. 4.3-2b 

938. 

572. 

77. 

2.04 

0.65 

0.95 

22. 

Non-Standard Vlnd 









Westerly 

Sumer 
(<200 fps) 
Ref. 3 

137. 

129. 

41. 

0.25 

0.02 

0.06 


The 1-term density modeling error is seen to be potentially a 
much greater contributor to vertical position error than the 
4-term model because of the larger uncertainties (20% vs 5%) 
associated with the 1-term atmospheric density model near 
145,000 ft. For the filter evaluated in this study, the ver- 
tical position error would more than triple with the 1-term 
model. This performance sensitivity underscores the need for 
an accurate onboard atmospheric density model for drag 
updating . 


Table 4.4-3 presents the baseline error budget for 
the navigation filter at 20,000 ft. The major contributors 
to navigation errors are: 

• The initial condition errors (Group 1) 

• The TACAN range bias error (Group 10) 

• The TACAN bearing bias error (Group 19) 

• The baro altimeter bias, scale factor 
and static defect errors (Group 12). 
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TABLE 4.4-3 

BASELINE ERROR BUDGET AT 20,000 FT 


T>0503 


ERROR SOURCE 


I. FILTER STATES AND MEASUREMENT 

rorsE 

iDltial Cooditions 
Pcsitioos and Velocities 
Uisaligoments 

Uncorrelated Measuremeot Noise 
TACAN Range 
TACAN Bearing 

Baro Altimeter 

INS Quantisation Noise 

II. IMU-REI ATSD STATES 

Accelerometer 

2. Biases 

3. Scale Factors 

4. Asymmetries 
Nonort bogona 1 i t ies 


5. 

Gyro 

7. 

8. 


Bias Drifts 
Mass Unbalanees 


IV. EXTERNAL AID^RELATED STATES 

10. TACAN Range Bias 

11. TACAN Range Scale Factor 

12. Bare Altimeter Errors 

Bias 

Scale Factor 
Static Defect 

First --Order Harkov 
16. TACAN Survey Errors 
19. TACAN Bearing Bias 
1. .21. ,22. ,23. 

Drag«Related Errors 


VALUE 


Ref. 47 
285 

100 ft 
12 mrad 

0.24 

0.0328 fps 


50 Ug 
100 ppm 
100 
15 ^ 

0.035®/hr 

0.025'*/hr/B 


38S ft 
100 ppm 

100 ft 
3% 


1.52 » 10-J 
ft/ftVsec^ 
20 ft 
I ft 
12 mrad 


Table 4.4*1 


RMS NAVIGATION ERROR AT 20,000 ft 


POSITION (ft) 


DR 


CR 


55.3 

49.9 


(3.88) 

C.2.2?> 

(3.7P, 

(ISO) 

110.5 


0.99 

0.63 

0.30 

eTT 

122.9 

dill) 

0.45 

0.54 

0.60 

0.3 

0.2 

0.1 

0.00 

0.00 

0.00 

5.8 

9.4 

1.9 

0.21 

0.14 

0.12 

75.8 

127.9 

20.0 

Cl .923(T7^ 

0.86 

23.8 

52.5 

10.3 

0.35 

0.67 

0.46 

25.7 

56.4 

10.8 

0.26 

0.76 

0.45 

8.0 

21.6 

2.7 

0.26 

0.40 

0.17 

15.7 

20.9 

13.0 

0.95 

0.65 

0.82 

2.9 

5.4 

1.0 

0.17 

0.16 

0.13 


C^44^)fa65TD( 31471) 
7.0 9.8 0.1 


(319.9)7284 70) (SOTT 
10.0 8.0 3.1 


6.5 16.6 5.7 


VELOCITY (fps) 


R 


DR 


CR 


z3Dl 


^ . 

0.22 0.32 0.01 


(4.00) 0.91 0.44 


0.08 0.05 0.03 

0.01 0.00 0.03 

1.16 1.07 I.IC 

0.60 0.31 0.29 


Total Projected Performance 


497.1 698.4 273.4 


8.57 5.83 6.73 


Touchdown Specifications 


3.0 RO.O 


4.7 


0.20 3.00 2 00 


Note that the large initial condition errors con- 
tinue to be major contributors to the velocity errors. The 
TACAN range and bearing bias errors are principal contri- 
butors to all components of position error, and the range 
bias is also a contributor to all components of the velocity 
error. The baro altimeter errors contribute significantly 
to position errors and to the vertical velocity error. Addi- 
tionally, TACAN range and TACAN bearing measurement noises 
show large effects in the vertical and crossrange position 
errors, respectively. Accelerometer bias errors continue 
to propagate into significant vertical and downrange velocity 
errors. 
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The Shuttle is approximately midway through its TAEM 
alignment maneuver at 20,000 ft. TACAN range and bearing 
information collected prior to the initiation of the maneuver 
becomes crossrange and downrange information, respectively, 
when the Shuttle is at or beyond the 20,000 ft mark because 
of the right angle turn. This accounts for the contribution 
of the TACAN range bias to crossrange errors, etc. The 
large contributions of the TACAN range bias to vertical 
channel errors, and the baro altimeter errors to downrange 
channel errors can be attributed to the high correlations 
which exist between the vertical and downrange channels (at 
deorbit, vertical and downrange errors are perfectly corre- 
lated) . 


Appendix C contains detailed computer printouts of 
all error contribution time histories corresponding to 
Tables 4.4-1 through 4.4-3. 

4.4.2 Navigation System Performance: 20,000 ft to 

Touchdown 


Figures 4.4-2a and 4.4-2b present overall performance 
curves from 20,000 ft to touchdown. At T ^^j g. there are large 
decreases in all three error components, reflecting the high 
accuracy of the MLS measurements as compared to TACAN and 
baro altimeter measurements. Subsequent to T ^^ ,, the position 
errors become essentially the navigation accuracy of MLS. 

The initial radar altimeter measurement results in an increase 
in vertical position error due to the large terrain modeling 
error for Vandenberg AFB, as was also observed in the Ref. 1 
study. 


All components of the velocity error decrease with 
MLS measurement processing until the switch to the comple- 
mentary filter at Tg^. The downrange velocity error 


/07 
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R-24«32< 



Figure 4.4-2a Overall Performance from 20,000 ft 

to Touchdown: Position 



Figure 4.4-2b Overall Performance from 20,000 ft 

to Touchdown; Velocity 

RlvPRODUCIBlLITY OF THE 
OHKHNAL PAGE IS POOR 
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increases after Tg^ due to the effects of the MLS range 
second-order markov error. The unstable growth in vertical 
velocity error at observed in Ref. 1 does not occur here 
since the vertical channel is not coupled with the downrange 
or crossrange channels via the MLS gains and since there are no 
vertical velocity updates after (see Table 4.2-2). Hence 
neither MLS range errors nor the increase in the vertical 
position error due to the terrain modeling error drive the 
velocity errors up. Subsequent to there is a slight 
increase in the vertical velocity component due to IMU mis- 
alignments . 


Table 4.4-4 is the baseline error budget at Tgy^. 
The major contributors to navigation errors are: 

• The MLS azimuth Lud DM£ measurement 
noises (Group 1) 

• The MLS scale factors and biases (Group 13) 

• The MLS second -order markovs (Group 14) 

• The MLS timing biases (Group 20). 


The significance of MLS measurements to the filter 
is clearly evident by the fact that nearly all the major con- 
tributors are MLS-related. The MLS azimuth error sources 
dominate both position and velocity crossrange errors, as do 
the MLS range error sources in the downrange channel. The 
steepness of the descent path, approximately a 20-24 deg 
glide slope, causes the downrange error sources (MLS range 
errors) to contribute large vertical channel errors. These 
MLS range-related error sources along with MLS elevation- 
related errors dominate the vertical channel errors. 


The MLS second-order markov errors (Group 14) show 
particularly large effects on the velocity errors. These 
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TABLE 4.4-4 

BASELINE ERROR BUDGET AT 


T-OS<M 






RMS NAVIGATION ERROR AT 

13.S8C ft 1 



ERROR SOURCE 

VALUE 

POSITION (ft: 

VELCX:XTT (fps) 1 





R 

DR 

_CR 

R 

DR 

CTR 

1. 

PILTEB STATES AND NEASUREMENT 



















Initial CoDilitioa& 









Positions and Velocities 

2BS ^ 

1.2 

0.6 

1.3 

0.09 

0.34 

COD 


Misaligmaents 

2.8 

4.8 

7.4 

0.18 

0.18 


Dncorrelated Neasureaent Noise 









TACAN Range 

100 ft 

0.5 

0-7 

0.4 

0.00 

0.01 

0.01 


TACAN Bearing 

12 arad 

0.3 

1.3 

1.4 

0.01 

0.01 

0.02 


Bare Altimeter 

0.S4 fx 

0.0 

0.0 

di) 

0.00 

0.00 

JLJKL 
Co. 46) 


HLS 

Azimuth 

.0.2 mrmd 

1.0 

3.8 

0.03 

0.03 


NLS 

HLS 

Elevation 

DUE 

0.1 arad 
40 ft 

2.9 0.7 

mT?>(25775 

0.2 

2,7 

0.18 0.05 

t0JM)(1.76) 

C.Ol 

0.08 


INS Quantization Noise 

0.0328 fps 

0.1 

0.1 

0.2 

0.01 

0.01 

0.01 

11. 

lUD-BELATED STATES 





- 




Accelerometer 









2. 

Biases 

50 tig 

G.8 

1.7 

3.1 

0.05 

0.07 

0.15 


3. 

Scale Factors 

100 ppa 

0.1 

0.2 

0.5 

0.03 

0.05 

0.01 


4. 

Asymmetries 

100 DOB 
15 ^ 

0.2 

0.2 

0.6 

0.03 

0.05 

0.02 


b. 

Nooort hogonal i t i es 

0.2 

0.3 

0.7 

0.04 

0.01 

0.03 


Gyro 










7. 

Bias Drifts 

0.03S“/hr 

0.8 

1.4 

1.9 

0.08 

0.05 

0.08 


8. 

Hass Unbalances 

0.025^/hr/g 

0.2 

0.3 

0.4 

0.02 

0.01 

0.02 

IV. 

EXTERNAL AID-RELATED STATES 









10. 

TACAN Range Bias 

365 ft 

S.5 

8.6 

7.3 

0.10 

0.05 

0.17 


11. 

TACAN Range Scale Factor 

100 ppm 

0.2 

0.4 

0.6 

0.02 

0.00 

0.01 


12. 

Baro Altimeter Errors 










Bias 

100 ft ) 
3% of alt. f 
1.S2 X 10-*( 









Static Defect 

0.7 

7.5 

4.5 

Co. 383 0.13 

0.05 




tX/U^/aec^ ) 









First-Order HarFov 

20 ft 

0.0 

0.0 

0.0 

0.00 

0.00 

0.00 


13. 

MLS Scale Factor and 
Biases 
Range Bias 










60 ft 

(1872)^0) 

5.6 0 773) 

0.1 

0.05 

0.04 

0.13 



Range Scale Factor 

400 ppa 

(25.7) 

0.10 

0.21 

0-06 



Azimuth Bias 

0.6 arad 

0.3 

0.8 

0.06 

0.08 Ccrrw)\ 



Elevation Bias 

0.4 arad 

03.63 

3.6 

0.4 

0.23 

0.09 



14. 

NLS Second-Order Harkovs 










Range 

40 ft 

C16.6) 


3.2 

C0.90K2.12) 

0.09 



Azimuth 

0.3 mrad 

6.5 cis.n 

Q06 

0.05 (O.Wl 



Elevation 

0.3 arad 

(i3tn 

3,0 

1.0 

Co. 78) 0.24 

0.04 


17. 

MLS Survey Errors 

1 ft 

1.0 

0.2 

1.0 

0.00 

0.00 

0.00 


19. 

TACAN Bearing bias 

12 arad 

0.5 

2.8 

2.3 

0.02 

0.03 

0.03 


20. 

HLS Timing Bias 










Range 

50 msec 

CO) (20; 

0.0 

0.03 

0.10 

0.06 



Azimuth 

50 msec 

0.9 

1.2 

4.3 

0.05 

0.05 (0.65)1 



Elevation 

SO msec 

2.0 

0.5 

0.1 

0.03 

0.02 

0.00 

1 Total Projected Performance 


37.0 

81 .4 

37.9 

1.46 

2.81 

1.30 

Touchdown Specifications 


3.0 

80.0 

4.7 

0.20 

3-00 

2-00 


//^ 
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error sources are primarily due to multipath effects, but 
little data on the magnitudes or frequency characteristics for 
multipath is available. Such large velocity error contribu- 
tions due to multipath effects suggest that an improved error 
model based on more extensive test data is necessary to 
accurately assess landing navigation performance. 

The IHU misalignments, which have been major contri- 
butors all through the higher altitudes, are lesser contri- 
butors at Tg^. The baro altimeter continues to be a major 
contributor to vertical velocity, and the TACAN range bias 
still makes considerable contributions to position errors. 

Several error sources found in the past (Ref. 1) to 
be minor contributors to navigation errors at this and lower 
altitudes have been ignored in Table 4.4-4 and all subsequent 
baseline error budgets. These include the drag-related errors 
and the TACAN survey errors. 

Table 4.4-5 presents the baseline error budget at 
100 ft above the runway touchdown point , just prior to 
the first radar altimeter measurement. The major error con- 
tributors are; 

• The initial conditions (Group 1) 

• HLS errors and measurement noises 
( Groups 1 , 13 , 14 and 20 ) 

• Baro altimeter errors (Group 12). 

After T„„, with the suboptimal gain structure of 
the complementary filter, the misalignment terms experience 
a marked increase in importance. The initial position and 
velocity errors show an increased error contribution after 
Tg^ through their coupling with the misalignments via the 
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TABLE 4 • 4-^5 

BASELINE ERROR BLT)GET AT 100 FT ABOVE RUNWAY 


T^0505 


ERBOR SOURCE 


VALUE 


RMS NAVIGATION ERROR AT 100 ft 


POSITION (ft) 


R 


DR 


CR 


VELOCITY (fps) 


R 


DR 


CR 


II. 


FILTER STATES AND MEASUREMENT 
tolSE 

Initial Conditions 
Positions and Velocities 
Misallgaments 

Oncorrelated Measurement Noise 
TACAN Range 
TACAN Bearing 

Bare Altimeter 
MLS Azimuth 
MLS Elevation 
MLS DME 

IRS Quantization Noise 

IMO^RELATED STATES 
Accelerometer 


Hef. 47 
28S aSt 

100 ft 
12 Brad 
b/24X70 

0.24 e 

0.2 Brad 
0.1 mrad 
40 ft 

0.0328 fps 


ft 



0.01 

0.01 


0.00 

0.00 


0.01 

0.01 



2. 

Biases 

50 ug 

1.2 

1.7 

1.3 

0.12 

0.12 

0.09 

3. 

Scale Factors 

100 ppm 

0.4 

0.1 

0.4 

0.05 

0.06 

0.03 

4. 

AsynsDetries 

100 ppm 

0.4 

0.3 

0.4 

0.05 

0.06 

0.03 

5. 

Honor t hogona 1 i t ies 

IS sec 

0.2 

0.5 

0.3 

0.03 

0.00 

0.02 

Gyro 









7. 

Bias Drifts 

0.03S®/hr 

0.9 

1.7 

0.6 

0.10 

0.08 

0.07 

8. 

Mass Unbalances 

0.025°/hr/g 

0.1 

0.4 

0.1 

0.02 

0.01 

0.02 

IV. EXTERNAL-AID RELATED STATES 








10. 

TACAN Range Bias 

385 ft 

0.3 

10.2 

dlD 

0.15 

0.04 

0.08 

11. 

TACAN Range Scale Factor 

100 ppm 

0.1 

0.6 

0.3 

0.02 

0.01 

0.00 

12. 

Bare Altimeter Errors 









Bias 

ion ft 1 








Scale Factor 

3T of alt ,, 1 

1 (TD 10.3 

1.6 

CO. 52) 

0.39 

0.04 


Static Defect 

i . 62 a 10 ^ j 









ft/ft^/sec^ I 








First-Order Markov 

20 ft 

0.0 

0.0 

0.0 

0.00 

0.00 

0.00 

13. 

MLS Scale Factor and 





1 




Biases 
Range Bias 

60 ft 

0.8 

dO) 

1.3 

Co. 48) 0.03 

0.02 


Range Scale Factor 

400 ppm 

0.5 

5.2 

9.2 

0.08 

0.18 

0.00 


Azimuth Bias 

0.6 mrad 

0.1 

0.7 

(8^ 

0.03 

0.01 

Co. 27) 


Elevation Bias 

0.4 mrad 

1.1 

0.2 

0.1 

0.20 

0.03 

0.00 

14. 

HLS Second-Order Harkovs 









Range 

40 ft 

1 C3^B}C42^43 

1.3 

lC0.41)(2.96) 0.06 


Azimuth 

0.3 mrad 

1 0^ 

7.3 

(ill) 

0.08 

0.02 CO. 45) 


Elevation 

0.3 mrad 

dif) 

0.7 

074 

CO, 23) 

0.05 

0.01 

17. 

MLS Survey Errors 

1 ft 

1.0 

0.0 

1.0 

oM 

0.00 

0.00 

10. 

TACAN Bearing Bias 

12 mrad 

0.1 

3.2 

0.9 

0.02 

0.00 

0.02 

20. 

MLS Timine !:ias 

i 

1 







nange 

50 msec 

0.5 

0.5 

0.20 

0.31 

0.00 


Azimuth 

i 50 msec 

0.1 

2.0 

0.6 

1 0.04 

0.01 

0.00 


Elevation 

50 msec 

1.1 

0.0 

0.0 

(o3D 

0.03 

0.00 

Total Projected Performance 


8.2 

80.5 

12.2 

i 1.11 

3.63 

0.70 

Touchdown Specifications 


3.0 

80.0 

4.7 

0.20 

3.00 

0 

0 
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mechanism described in Section 4.4.1. The dominant velocity 
error contributors continue to be the BttiS second-order 
markov errors. In particular, the 40 ft error of the range 
second-order markov couples directly into velocity errors. 

The MLS range and azimuth measurement noises and 
range and azimuth biases affect the downrange and crossrange 
error components about as expected. Even though the down- 
range velocity error grows following Tg^^, the suboptimal 
filter appears to select reasonable MLS gains in the down- 
range and crossrange channels. In the vertical channel, 
however, the baro altimeter errors contribute at least as 
much as the MLS errors. In addition, MLS range errors remain 
important contributors to vertical velocity even though the 
final approach following the flare maneuver is made with a 
low (1-3 deg) glide slope. This behavior suggests that the 
complementary filter underweights the MLS elevation measurement . 

The baseline error budget at touchdown is given in 
Table 4.4-6. The error in vertical position slightly exceeds 
the touchdown specification — with nearly all of the vertical 
position error due to the terrain modeling error for Van- 
denberg's runway 30 (other landing sites are discussed later). 
The downrange position error is nominally equal to the 
touchdown specification, and the crossrange position error 
is approximately twice the specification — both totals 
primarily due to MLS error sources. The MLS biases and 
second-order markovs are again the dominant error sources 
in the downrange and crossrange channels. 

The vertical velocity error exceeds the touchdown 
specification; however, it is significantly smaller than that 
observed in Ref. 1. Since the radar altimeter measurement is 
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TABLE 4.4-6 

BASELINE ERROR BUDGET — TOUCHDOWN 


T-0506 





RMS NAVIGATION ERROR AT TOUCBDOVN 


HBBOR SOURCE 

VALUE 

POSITION (ft) 

VELOCITY (fps) 1 




R 

DR 

CR 

R 

DR 

CR 

I* 

Initial ConditioBS 
Posltioas and Velocities 

S«f . 

0.1 

2.2 

1.1 

_ 

0.37 

OD 


Misaligmients 

285 SSh 

0.0 

5.2 

CD 

C6~48J 

0.26 

C6 . 22v' 


Uncorrelated Noise 








Radar Altisteter 

0.9 ft 

0.3 

0.0 

0.0 

0.00 

0.00 

0.00 


!2LS Aaimuth 

0.2 mrad 

0.0 

4.2 

2.0 

0.04 

0.01 

0.C2 


MLS Elevation 

0.1 mrad 

0.0 

_0.J2 

0.1 

0.04 

0J)1 

0.00 


MLS DME 

40 ft 

0.0 


0.9 

0.21 

atow 

0.00 


RSS for Other External 


0.0 

1.6 

0.5 

0.05 

0.03 

0.02 


Navigation Aids 








11. 

IMU*Related States 

Fable 4.4-S 

0.0 

2.5 

1.4 


0.17 

0.12 

IV, 

External Aid-^Related States 
13. MLS Scale Factor and Biases 




CUD 





Azinmth Bias 

0.6 mrad 

0.0 

0.7 

0.03 

0.01 

COT) 


Elevation Bias 

0.4 mrad 

0.0 

0.2 

0.1 

0.21 

0.02 

0.00 


Range Bias 

60 ft 

0.0 


1.4 

foT4?) 

0.01 

0.03 


Range Scale Factor 

400 ppm 

0.0 

4.2 

0.3 

0.09 

0.16 

0.00 


14. MLS Second-Order Harkov 




CO) 



(F3?) 


Anisutb 

0.3 mrad 

0.0 

7.3 

0.09 

0.02 


Elevation 

0.3 mrad 

0.0 

0.7 

0.4 

0.02 

0.04 

0.01 


Range 

40 ft 

0.0 

C44 . 4^ 

1.2 

(3745) 

C3T3T) 

0.06 


15. Radar Altimeter Errors 


m 







Bias 

1 ft 

0.0 

0.0 

0.00 

0.00 

0.00 


Terrain Modeling 
Error (VA53-30) 
All Instrument and 

Ref. 1 

0.0 

0.0 

0.00 

0.00 

1 

0.00 



Other Errors 

Ref. 1 

0.0 

0.0 

0.0 

0.00 

0.00 

0.00 


10.-12.. 17,-20. 




CO) 

CoTW) 




RMS for Other External 

Table 

0.1 

(STT) 

0.59 

0.09 


Navigation Aids 

4.4-5 






Total Projected Perloraance 


3.8 

60.6 

10.3 

1.42 

3.73 

0.65 

Touchdown Specifications j 

3.0 

80.0 

4.7 

0.20 

3.00 

2.00 


not used to update vertical velocity, and MLS does not up- 
date the vertical channel after the MLS elevation measure- 
ment is lost, the vertical velocity errors show a slight 
increasing trend after Tj^. On the other hand, these same 
factors prevent the excessive increase in vertical veloc- 
ity error after Tj^^ as observed in Ref. 1. That is to say, 
the large vertical position error due to terrain modeling 
errors at Tj^^ does not propagate into large velocity errors 
in any channel. 


//H 
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The downrange velocity error exceeds the touchdown 
specification slightly. The dominant error source is the 
MLS range second-order markov, which again underlines the 
need for a better understanding of the frequency and magni- 
tude of this error source and may warrant an investigation 
of techniques for minimizing its effects. The crossrange 
velocity is within the touchdown specification. The major 
contributors to this error are the MLS azimuth bias and 
second-order markov, and the initial conditions. 

The terrain modeling error included in Table 4.4-6 
is for runway 30 at Vandenberg AFB. The contributions of the 
terrain modeling errors to vertical position errors for the 
other primary landing sites and for the secondary sites at 
Andersen AFB and Eickam AFB are presented in Table 4.4-7. 

Only Vandenberg 's runway 30 and Andersen's North runway 
(southwest approach) yield vertical position errors at 
touchdown which exceed the runway specification. 


TABLE 4.4-7 


CONTRIBUTION OF TERRAIN MODELING ERROR 
TO VERTICAL POSITION ERROR AT TOUCHDOWN 

T-0535 


LA.VDING SITE 

R.MS VERTICAL 
POSITION ERROR 
AT TOUCHDOWN (ft) 

LA.\T)ING SITE 

RF.S VERTICAL 
POSITION ERROR 
AT TOUCHDOWN (ft) 

Edwards AFB 


Andersen AFB 


Runway ; 4 

0.4 

North Runway: 


17 

0.0 

NE Approach 

0.5 

22 

0.0 

SW Approach 

3.1 

35 

0.9 

South Runway: 


Kennedy Space Center 


r?E Approach 

2.6 

Runway : 15 

1.0 

SW Approach 

2.0 

33 

0.0 



Vandenberg AFB 


Hickar. AFB 


Runway : 12 

2.1 

Runway : 8 

0.0 

30 

3.7 

26 

0.0 
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Appendix D contains detailed computer printouts of 
all error contribution time histories corresponding to 
Tables 4.4-4 through 4 .4-7. Of note are the large misalign- 
ments observed in the tables for Group 1 (Initial positions 
and velocities, and misalignments). Group 2 (accelerometer 
biases). Group 7 (gyro bias drifts) and Group 12 (baro altimeter 
bias, scale factor and static defect). These errors are 
transients due to the filter's attempt to estimate platform 
misalignments based on a combination of MLS measurements 
and the assumed correlations between the misalignments and 
velocity errors at 20,000 ft. In order to demonstrate 
that these errors are in fact transients, a partial error 
budget was generated using a 12 state Kalman filter down 
to the runway threshold. The filter substantially reduced 
the misalignments between the baseline Tg^(alt*12,000 ft) 
and the runway threshold. Elimination of some of the 
filter difficulties discussed in Section 4.4.1, such as an 
overly pessimistic filter covariance above 20,000 ft, 
should improve the filter-predicted correlations between 
velocity errors and misalignments and thus reduce the 
transients. 

In summary, the navigation filter performs reasonably 
well down to touchdown . The transition from a Kalman to a 
complementary filter is made without significant transients 
in the estimates, and the fixed MLS and radar altimeter gains 
maintain reasonable navigation accuracy. The only significant 
departure from the touchdown specification is in vertical 
velocity, although downrange velocity and all three components 
of position are also out-of-spec. There is evidence to 
suggest that the velocity estimates at touchdown and the 
vertical position performance prior to the first radar 
altimeter measurement could be improved with a more judicious 
choice of complementary filter gains, but it is doubtful that 
the vertical velocity specification could be met. 
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4.5 ENTRY, PRELAND AND LANDING SENSITIVITY ANALYSIS 

This section contains several curves indicating the 
sensitivity of navigation system perfonnance to variations 
in truth model error source statistics (while holding filter 
parameters fixed) for the entry, preland and landing naviga- 
tion phases. These "fixed filter" sensitivity curves illus- 
trate the effect of an unknown variation in the rms value 
of an error source or group of error sources. 

The baseline error budget for the entry navigation 
phase. Table 4.4-1, indicates that the initial condition 
errors dominate the vertical and crossrange velocity errors. 
Figure 4.5-1 presents vertical velocity error sensitivity at 
145,000 ft to initial position and velocity errors; Fig. 4.5-2, 
the crossrange velocity error sensitivity to initial misalign- 
ment errors. An increase in the initial condition errors in 
either case results in substantial increases in the velocity 
errors — a factor of two worsening of initial condition errors 
would result in at least a 67% increase in the vertical velocity 
error and a 100% increase in crossrange velocity error. 

Halving the initial position and velocity uncertainty 
would yield only a 29% improvement in vertical velocity 
error. However, similarly reducing the initial misalign- 
ment uncertainty (to 140 sec) would yield a nearly 40% 
improvement in the crossrange velocity. 

At 20,000 ft the vertical velocity channel is domi- 
nated by initial condition errors and the bare altimeter 
bias, scale factor and static defect errors. The sensitivity 
of this error to the baro altimeter errors is presented in 
Fig. 4.5-3. Because of the presence of other major error 
sources, the effects of varying the baro altimeter uncer- 
tainties are not as dramatic as observed for the initial 
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R<24640 



RELATIVE INITIAL POSITION AND VELOCITY ERRORS 


Figure 4.5-1 RMS Vertical Velocity Error Sensitivity 

at 145,000 ft to Initial Position and 
Velocity Errors 


R-24641 



0 12 5 


RELATIVE INITIAL MISALIGNMENT ERRORS 


Figure 4.5-2 RMS Crossrange Velocity Error Sensitivity 

at 145,000 ft to Initial Misalignment 
Errors 
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condition errors in Figs. 4.5-1 and 4.5-2. Since the initial 
conditions errors contribute about as much as the baro altim- 
eter uncertainties to vertical velocity at 20,000 ft, the 
sensitivity curves for these error sources would be quite 
similar. 


Major error sources at Tj^^, 100 ft above the touch- 
down point, are the MLS-related errors. The sensitivity of 
the vertical velocity error to the MLS range bias uncertainty 
is shown in Fig. 4.5-4. Again, because of other error 
sources of equal importance, the vertical velocity error is 
seen to be be basically insensitive to changes in this error 
source. A sensitivity curve for the baro altimeter uncer- 
tainties, also a major contributor to the vertical velocity 
error at Tj^^, would be very similar to Fig. 4.5-4. 

In contrast to the vertical velocity error, the 
downrange velocity error at is dominated by a single 
error source — the MLS range second-order markov error, 'me 
sensitivity curve in Fig. 4.5-5 reveals that doubling the 
marKov uncertainty would result in a 77% increase in down- 
range velocity error; while a markov error of zero would 
yield nearly a 50% improvement in the velocity error. 

The final sensitivity curve. Fig. 4.5-6, presents 
the sensitivity of crossrange position error at touchdown 
to the MLS azimuth bias error, and indicates the appropriate 
touchdown specification. It is clear than even if the MLS 
bias error were zero, the touchdown specification still would 
not be met. This arises because the MLS azimuth second-order 
markov error contribution is itself larger than the specifica- 
tion. Therefore to meet this specification both MLS azimuth 
errors must be significantly reduced. 


Il<^ 
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R*24642 



Figure 4.5-3 RMS Vertical Velocity Error Sensitivity 

at 20,000 ft to Baro Altimeter Bias, 
Scale Factor and Static Defect Errors 


R.24638 



RELATIVE MLS RANGE BIAS ERROR 


Figure 4.5-4 RMS Vertical Velocity Error Sensitivity 

at 100 ft to MLS Range Bias Error 


I aft 
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0 12 5 


RELATIVE MLS RANGE SECOND-ORDER MARKOV ERROR 

Figure 4.5-5 RMS Downrange Velocity Error Sensitivity 

at 100 ft to MLS Range Second-Order 
Markov Error 


R.24637 



RELATIVE MLS AZIMUTH BIAS ERROR 


Figure 4.5-6 RMS Crossrange Position Error 

Sensitivity at Touchdown to MLS 
Azimuth Bias Error 


/•ll 


4-42 




THE ANALYTIC SCIENCES CORPORATION 


SUMMARY AND CONCLUSIONS 


A detailed analysis of Space Shuttle navigation has 
been presented. The study encompassed navigation performance 
during each of the major mission phases — from prelaunch IMU 
calibration and alignment, through ascent and orbital opera- 
tions. to entry and landing . The ascent and orbital opera- 
tions studies evaluated the potential effectiveness of GPS- 
aided navigation in the 1984't' time frame; the remaining 
studies considered navigation performance for the OFT-1 
mission in 1979. In each Instance, the analytical tech- 
niques and navigation error models used for the study were 
described. Detailed navigation error budgets were genera- 
ted for nominal mission profiles and performance sensitivity 
to the dominant error sources was evaluated. A summary of 
the most significant findings of the study is presented in 
the following subsections. 


5.1 ASCENT AND ORBITAL NAVIGATION 

The objective of the ascent navigation study was to 
provide NASA with an assessment of potential navigation 
accuracy in the operational GPS (1984+) time frame. The 
Shuttle was assumed to be equipped with a multichannel 
X-receiver capable of processing the high accuracy (P-code) 
satellite transmissions. An 11-state navigation filter was 
designed and a candidate measurement sequence was defined. 
The measurement sequence provided GPS measurement updates 
to the navigation state vector at a low update rate — range 
and range-rate measurements from one GPS satellite every 
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6 sec. This update rate is significantly below the capa- 
bility of the GPS receiver, but was selected to minimize the 
computational burden on the Space Shuttle computer. No 
significant perfonnance Improvement would be anticipated by 
adopting a higher update rate. 

The results of the ascent study are presented in 
Table 5-1. Relative to pure inertial navigation. GPS-aiding 
can provide a factor of 100 improvement in position esti- 
mation at orbit insertion and a factor of 30 improvement in 
velocity estimation . Furthermore, GPS measurements can be 
used to align the IHU during ascent. Although niU errors 
are the major contributors to the remaining misalignment 
uncertainty, the error budget upon which Table 5-1 was 
based clearly indlcateii that the navigation performance is 
completely dominated by GPS-related errors. This dominance 
of GPS-related error sources Implies that navigation per- 
formance is insensitive to an increase in the magnitude of 
any other error source as long as GPS is available . 

TABLF 5-1 

EFFECT OF GPS IN REDUCING ORBIT INSERTION ERRORS 


T-0269a 


IMPLEMENTATION 

RMS POSITION 
ERRORS (ft) 

RMS VELOCITY 
ERRORS (fps) 

RMS 

MISALIGN!^^ 
ERRORS (lec) 

DR 

CR 

R 

DR 

CR 

R 

DR 

CR 

R 

Pure Inertial 
Navigation* 

2280 

2730 

1864 

3.2 

7.6 

7.4 

64 

62 

61 

GPS-Aided 

Navigation 

7 

9 

16 

O.l 

0.1 

0.? 

25 

13 

15 


*Errors at midpoint of final Orbital Maneuvering System (OMS) 
insertion burn, 142 sec prior to insertion. 
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An important conclusion to be drawn fr<MB the ascent 
navigation study is that, idien <g*S-aidlnig is used, ascent 
navigation accuracy is not sensitive to the quality of the 
IMU or to how accurately the IMU is calibrated and aligned . 
*nie role of the IMD is to provide an inertial reference with 
short term stability so that GPS measuronents taken at dif- 
ferent times can be "tied together" to generate an accurate 
navigation estimate. This requirasent could be met by the 
Space Shuttle IMU without an extensive prelaunch calibration 
and alignment. 

The present baseline for operational Space Shuttle 
missions specifies the Tracking and Data Relay Satellite 
System (TDRSS) as the principal navigation aid for orbit 
determination. Table 5-2 compares orbit determination 
accuracy for TDRSS with the results jof the GPS-aided orbi- 
tal navigation study. The table clearly implies that GPS 
can provide significantly better navigation accuracy than 
TDRSS frcmi shorter measurement sequences. Although the GPS 
results were based on a 2 min measur^ent sequence, compar- 
able accuracy could be obtained with only a few seconds of 
data if the navigation filter were optimized. These 


TABLE 5-2 

PROJECTED ORBIT DETERMINATION ACCURACIES 
FOR TDRSS AND GPS 


NAVIGATION 

AID 

LENGTH OF 
MEASUREMENT 
SEQUENCE 
(min) 

RSS POSITION 
ERRORS 
(ft) 

RSS VELOCITY 
ERRORS 
(fps) 

TDRSS 

5 

>200* 

Not Specified 

GPS 

2 

20 

0.2 


’'‘Data interpretation subject to qualification discussed 
in Section 2.7. 
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performjice reBults, coupled with the operational advantage 
that (gS doeg not require ground processing and subsequent 
uploading of the navigation solution, clearly favor CTS as 
the more desirable orbital navigation aid. 


5.2 PBELAUNCH CALIMtATIQN AND ALIGNMENT 

The prelaunch calibration and alignment study was a 
continuation of the perfomiance phase of the IHD Software 
Verification and Validation Program initially reported in 
Sef. 1. The emphasis was on develoimient of a covariance 
analysis siamlation for the OFT software. This program 
will be used to generate a detailed cal/align error budget 
in a subsequent study. 


5.3 ENTRY, PRELAND AND LANDING NAVIGATION 

The principal finding of the entry, preland and 
landing navigation study is that the navigation filter as 
described in Ref. 45 incorporates measur^ents from the 
navigation aids into the navigation solution in a reason- 
ably effective manner, particularly after the transition 
to landing navigation at 12,000 ft. Previous studies by 
TASC (Ref. 1) had uncovered a tendency for the velocity 
estimates generated by the filter to diverge following the 
switch to the complementary mode because of a 

dynamic coupling between the vertical and downrange esti- 
mates. The filter presented in Ref. 45 was modified by 
NASA to eliminate this coupling. As shown in Fig. 5-1, 
downrange and vertical velocity errors increase somewhat 
during the landing phase, but the overall result is stable 
velocity estimation performance. 
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Figure 5-1 Velocity Performance During Landing — 
Complementary Mode Begins at 

The results of this and previous TASC studies suggest 
that the basic navigation filter design has the potential for 
achieving near-specification performance at touchdown, however . 
careful attention must be given to filter parameter and gain 
selection in order to attain this level of performance . The 
candidate entry and preland parameter values provided by NASA 
are quite pessimistic and result in a number of major naviga- 
tion performance anomalies, e.g., a sudden increase in naviga- 
tion errors following TACAN acquisition. Similarly, although 
the landing filter prevents velocity estimation divergence, 
there is considerable evidence to suggest that appropriately 
selected gains could improve velocity estimation accuracy dur- 
ing landing. A thorough parameter tradeoff study could result 
in an Improvement in entry, preland and landing navigation 
performance over that predicted in this study. 

If subsequent analysis shows that satisfactory 
navigation performance cannot be attained with the present 


)A(, 
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filter/navigation system combination, an alternative mould 
be to replace the navigation aids used prior to MLS acqui- 
sition with GPS. GPS-aiding during entry and preland could 
potentially yield inq»roved navigation performance with a 
simpler navigation filter. A quantitative evaluation of 
this alternative will be ccunpleted by TASC during the next 
contract period. 
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APPENDIX A 

ANALYTIC EVALUATION OF TWO MECHANIZATIONS 
OF THE DRAG UPDATE FILTER 


The drag update filter used in the Space Shuttle 
entry navigation system during blackout employs a type of 
pseudo-measurement in order to improve the on-board estimate 
of position and velocity. To date, an atmospheric drag 
pseudo-measurement has been assumed for the filter; however, 
the alternative of using an altitude pseudo-measurement 
exists. The Kalman filter used in the computer navigation 
program requires the definition of a measurement matrix Hp 
for the pseudo-measurement. Additionally, if the per- 
formance of the drag update filter is to be evaluated using 
a "truth model," the measurement matrix Hg for the truth 
model must also be defined. Two schemes for finding Hg , 
one employing drag pseudo-measurements and the other employing 
altitude pseudo-measurements, are described and compared in 
this appendix. 


A.l 


DESCRIPTION OF MECHANIZATIONS 


The atmospheric drag pseudo-measurement Q^(i) is 

defined to be a drag acceleration measurement constructed from 

the IMU accelerometer outputs. The expected (i.e., predicted) 

drag acceleration q (t) is computed from the current estimate 

P 

of position and velocity, a nominal atmospheric density model 
and estimates of the Shuttle's aerodynamic coefficients. The 
difference between residual 


'AS 


A-1 


(A. 1-1) 
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which appears in the update equation for the truth model 
error state vector 

6xg(t+) - 6Xg(t-) - Kg (t) 6qg(t) (A. 1-2) 

q 

where 6xg(t~) is the estimation error prior to the update 
and 6xg(f*') is the estimation error after the update. This 
mechanization is depicted in Fig. A. 1-1. 



RESOLUTION ALONG |[ <»m 

REUTIVE 

velocity 


R*20S39 


STATE 1 

ESTIMATE 


atmospheric 

1 

DENSITY 


MODEL 

1 



SHUTTLE ' 

\ 

AERODYNAAAIC 


model 



Y* 


3 ^ 

Figure A. 1-1 Simplified Conceptual View of Drag 

Pseudo-Measurement Mechanization 


The matrices of first partial derivatives Hm^Ct) and 
Hpq(t) relate 8-nd qp(t), respectively, to the truth 


♦Figures A. 1-1 and A,i-2 are intended to merely indicate 
conceptually what functional blocks are involved in gen- 
erating the pseudo-measurements versus the predicted 
values. These figures in no way attempt to present the 
detailed steps necessary to develop any of the variables 
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model state vector. The measurement matrix for the truth 
model is then defined by 

H (t) - H (t) - H (t) (A. 1-3) 

q q 

An alternative approach employing altitude pseudo- 
measurements hg^(t) and expected altitude h^(t) can be used 
to generate the truth model measurement matrix. In this 
scheme the measured altitude is computed from IMU acceler- 
ometer outputs, current position and velocity estimates, a 
nominal atmospheric density model and estimates of the 
Shuttle’s aerodynamic coefficients; while the expected altitude 
is interpreted from the current estimate of Shuttle position. 
Similar to Eqs. (A. 1-1) and (A. 1-2) the residual and error 
update equations are 

5hg(t) « hjjj(t) - hp(t) (A. 1-4) 

53^(t+) = 6Xg(t“) - Kg (t) 6hg(t) (A. 1-5) 

h 

Figure A. 1-2 depicts this scheme. 



Figure A. 1-2 Simplified Conceptual View of Altitude 

Pseudo-Measurement Mechanization 
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Hnj^(t) and Hp^(t) relate and hp(t), respectively, to 

the truth model state vector, and Hsjj(t) is the truth model 
measurement matrix given by 

(t) = (t) - (t) (A. 1-6) 

®h Ph “h 

This appendix presents a comparison of these schemes, 
and answers the question of whether or not the system error 
sources have the same effect on the performance for both 
mechanizations. Table A. 1-1 lists the truth model states 
used in evaluating the drag update filter. 

TABLE A. 1-1 

DRAG UPDATE TRUTH MODEL STATES 


Error Source Names 

State 

I. ESTIWATED STATES AND UNCORRELATED NOISES 


Group 1. 


Position Errors 

1-3 

Velocity Errors 

4-6 

Platform Itisallgoments 

7-9 

Accelerometer 'Biases' (First-Order Barkovs) 

10-12 

Correlated Ueasurement Error (First Order Barkov) 


Drag 

13 

II. NON-ESTI MATED. IMU-RELATED STATES 


Group 2. Accelerometer True Biases 

14-16 

Group 3. Accelerometer Scale Factor Errors 

17-19 

Croup 4. Accelerometer Uisaiignments 

20-25 

Group 5. Accelerometer Nonlinearities 

26-28 

Group 6. Gravitational Deflections and Anomalies 

29-31 

Croup 7. Gyro True Bias Drifts 

32-34 

Group 8. Gyro Mass Unbalances 

35-40 

Group 9. Gyro Anisoelasticities 

41-43 

111. NON-ESTIMATED. DRAG RELATED STATES 


Group 21. Non-Standard Density 


1962 Standard Atmosphere Error 

44 

Time-Varying Bias ' 

45 

First-Order Markov 

46 

Group 22. Non-Standard Vind 


Westerly (Time-Varying Bias) 

47 

Headwind (First-Order Markov) 

48 

Crosswind (Second-Order Markov) 

49-50 

Turbulence (First- and Second-Order Markovs) 

Sl-55 

Group 23. Non-Standard Aerodynamics 


First-Order Markov 

56 
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A. 2 DRAG MEASUREMENT MECHANIZATION 

This section is essentially an outline of Appendix 
C of Ref. 2 , in which the truth model measurement matrix 
for the drag acceleration pseudo-measurement case is de- 
veloped. The text here presents only the highlights of 
that development. Detailed derivations and definitions can 
be found in the noted appendix. 

The drag acceleration pseudo-measurement extracted 
from the accelerometer outputs is 

= livR (A. 2-1) 

c 

where ivRc^^^ ^ unit vector in the direction of the com- 
puted relative velocity and f„(t) is the measured specific 
force. The difference between the true drag acceleration 
and that measured is 


= q(t) - q„,(t) (A. 2-2) 

Combining Eqs. (A. 2-1) and (A. 2-2) the perturbation 5q^(t) can 
be defined by 


Sq^jCt) = 6i^j^(t) • f(t) + iyjj(t) * 6f(t) (A. 2-3) 


where 


^ivR ^ ^ ^ ~ ivR ^ ^ ^ ~ ^ ^ ^ 

6f(t) = f(t) - f^(t) 


and iyjj(f) is ^ unit vector in the direction of the true 
relative velocity and ^(t) is the true specific force. The 
error 6^(t) in measured acceleration can be attributed to 
the standard IMU acceleration error sources defined in the 
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drag update truth model, namely, platform misalignments 
(Group 1) and accelerometer errors (Groups 2 through 5). 

Thus, the quantity ^^(t) • 6£(t) in Eq. (A. 2-3) introduces 
terms in Hpjq(t) relating to truth model states 7-9 and 14-28. 
The remaining term in Eq. (A. 2-3) is given by 


“ -lL(t)(£j. X 6r(t)) + Fj^(t) 6v(t) 

- iL^vt) V^.(t) (A. 2-4) 


where 


T 

f (t) 

F (t) = 


I ivR i-VR 

c c 

iv^(t) - Ee 


(A. 2-5) 


Er is the earth rate, V (t) is the atmospheric wind and R«(t) 

and V (t) are the computed position and velocity vectors, 
c 

respectively. Adding in these terms, relating to position, 
velocity and non-standard wind error states, completes (t) 

q 



(t) 


position velocity platform 

errors errors misalignments 


X) 


4 


I 

6 


T 

ivR T 


I/P 


FP 


I 

13 


accelerometer accelerometer 
biases scale factor errors 


I 

I 

I 

I 


T * 

4 * T » 

-VR ^I/P I 

I 

16 



I 

I 

1 


19 


accelerometer 

misalignments 


\ T 
4 T 
1 iVR ^ 

I 


I/P ^4 


accelerometer 

nonlinearities 


t i ^ T F ' 

I -VR M/P ^5 ; 
25 28 


winds 


! -L 


46 56 


(A. 2-6) 
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where Tj^p is the transformation matrix from platform to 
inertial reference coordinates; and F^, F„, F . , F- are 

Ir O H O 

matrices whose elements are the components of the specific 
force f(t) coordinatized in platform coordinates. 


The computation of Hc„(t) requires that (t) also 

q ^q 

be determined; the equation for 6qp(t) must therefore be 
found. The expected drag acceleration is 


* I 't) i ''r 

c c 

where Cj) (t) is the computed drag coefficient for the Space 
Shuttle, A is its cross-sectional area, and m is its mass, 
p (t) and Vo (t) are the computed atmospheric density and 

w c 

computed relative velocity, respectively. The difference 
between the true drag and that computed is 


•5qp(t) = q(t) - qp(t) (A. 2-8) 

From Eqs. (A. 2-7) and (A. 2-8) the relationship defining the 
perturbation 6q (t) is 

+ Cj)(t) |p(t) Vj(t) • «V^(t) (A. 2-9) 

which can be written in terms of the true drag as 


qp(t) = 


q(t) 


6Cp(t) 


6p 


ill 


2ivR 


(t) 




p(t ) 


Vj^(t) 


(t) ) 

> (A. 2-10) 


yay 
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where 


6p(t) = p(t) - P^(t) 

6Cjj(t) - Cjj(t) - Cjj (t) 

c 

5Vp(t) = Vjj(t) - Vjj (t) - 6v(t) - ilE X 6r(t) - V^(t) 


(A. 2-11) 

Considering the drag coefficient term in Eq. (A. 2-10), 6Cjj(t) 
can be approximated by 


6Cp(t) = Cp^(a,M,t) (A. 2-12) 

where Cj)j(a,M,t) is the difference between the true Space 
Shuttle arag coefficient and that predicted by the drag 
coefficient computational model. Assuming that the drag 
correlated error and the non-standard aerodynamics errors 
are defined as percent deviations from their nominal values, 
the terms of Hpg(t) corresponding to error states 13 (drag 
correlated measurement error of Group 1) and 56 (Group 23) 
can easily be written. The 6p(t)/p(t) term in Eq. (A. 2-10) 
is considered by first making the assumptions that the 
atmospheric density model from which P^,(t) is computed is 
a simple exponential function of computed altitude h^(t) 

hc(t) 

Pc(t^c) = hsc (A. 2-13) 

and that the true atmospheric density satisfies 


_ h(t ) 

P(R) = Pq^ hgc + P^Cr.) (A. 2-14) 
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where h(t) is the true altitude, p and h are constants and 

o sc 

Pj(R) is the iifference between the true atmospheric density 
and that predicted by the atmospheric density computational 
model. A first order expansion for 6p(t) yields 

ip(t) • 6r(t) 

6p(R) « “Pc<t) h + Pf(R) (A. 2-15) 

sc 

where ijj('t) is the unit vector along R(t). The dependence 
of the error in expected drag 6qp(t) on position and non- 
standard density errors is evident from Eqs. (A, 2-15) and 
(A. 2-10). The position error related terms in Hp^(t) follow 
directly from the first term on the right hand side of 
Eq. (A. 2-15). With the assumption that the non-standard 
density errors are defined as percent deviations from their 
nominal values, the terms in Hp^Ct) related to the Group 21 
error states follow from the p^(R) quantity in Eq. (A. 2-15) 
and from Eq. (A. 2-10). Substituting Eq. (A. 2-11) into the 
last term in Eq. (A. 2-10) completes the derivation of Hpq(t), 
adding terms related to the position, velocity and non- 
standard wind error states. 


H_ (t) - q(t) 


position 
errors 

■ 


sc 


velocity 

errors 



3 6 


0 


drag 

correlated 

error 


1 


0 


12 


13 43 


DOD-standard 

density error non-standard 
model wind 




1 


1 1 




I 

0 ! 


non-standard 

aerodynamics 

errors 


1 


46 


55 


(A 2-16) 
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A. 3 ALTITUDE MEASUREMENT MECHANIZATION 


Consider the following pair of equations in wnich 
the non-standard density error model term p^(R) has been 
ignored . 

_ hit! 

P(t) = p^e hgc (A. 3-1) 

q(t) = llvR(t) • i(t)l = m (A. 3-2) 

where all variables have been previously defined. Equations 
(A. 3-1) and (A. 3-2) represent the situation in which altitude 
and velocity are known precisely, the accelerometers are 
perfect and the computational models for the drag coefficient 
and atmospheric density are exact. Solving for h(t) yields 
the altitude pseudo-measurement equation 


h(t) = h_ log 


I V*) s 


V^(t) 


SC 


q(t) 


(A. 3-3) 


from which the relationship for 6h (t) can be easily 

P 

established 


6hp(t) 


(6C^(t 
■ •’sc ( Cp(t 




£Vjj(t) 


lytT 


6q(t) 

q(t) 


(A. 3-4) 


The first two terms on the right hand side of Eq . (A. 3-4') 
also appear in Eq. (A.i-10), and so affeci (t) as they 


affect Hp^( t ) . 


That is to say, these terms introduce into 


Hm.-(t) terms related to the position and velocity, non- 
"'h 

standard wind, non-standard aerodynamics and drag corre- 
lated error states. The last tenn in Eq . (A. 3-4) can be 
considered by first noting that q(t) = ,^^,p(t) • f_( t ) ; . 
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Then 6q(t) is given by Eq. (A. 2-3), which equation completely 
determined Thus, 6q(t) inserts into Hmj^(t) all the 

terms in Hm (t). 

The final item of concern in developing H|Q^(t) is 
to now consider a non-standard density error model term in 
Eq. (A. 3-1), namely 

_ h(t) 

p(R) = p^e hg^, + Pj(R) (A. 3-5) 

where again, as before, p^(R) is the difference between the 
true atmospheric density and that predicted by the atmospheric 
density computational model. Eq. (A. 2-15) gives 6p(R) as a 
function of 6hm(t) 

6h (t) 

6p(R) *= -P^,(t) + Pj(R) (A. 3-6) 

“sc 

where 6h(t) has been substituted in for ip(i) * 6r(t). This 
equation can be solved for 


6h„(t) * -h^^ ^ 

m' ' sc p^(t) 


Pf(R) 
sc P^(t) 


(A. 3-7) 


The P^(R)/P (t) term in Eq. (A. 3-7) indicates the effect of 

the non-standard density error model on ^hjjj(t). Note that 

this effect is the same as that exhibited by pj(R) on 6q (t) 

(from Eqs. (A. 2-15) and (A. 2-10)) in deriv’^ing Hj, (t). 

t'q 

Therefore the matrix completely established. 


/3i 
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“h 


sc 


position 
err rs 


^ivR 

'r * 


velocity 

errors 


platform 

misalignments 


^-VR 


I T 

^ 1 -VR 

9 I 


**’l/P 


drag accelerometer accelerometer 

correlated error biases scale factor errors 


1 1 

1 t 

1 n 1 1 

' ivR ^I/P 

I ivR^ ^I/P ^3 

1 

1 

1 

! U * X 

1 1 

t 1 

12 

1 4 

13 

1 ^ 
16 

t 

[ 

19 

accelerometer 

misalignments 

accelerometer 
non 11 near! ties 

non-standard 
density error 
model 


1 ivR ^1/P ^4 
Q 

IhlZuzhl 0 
! ! 

i > 

1 i 

1 1 1 1 1 
1 1 



25 


28 


43 


46 


non-standard 

wind 


2 ivR 


non-standard 

aerodynamics 

errors 


.h oi 

q ‘ 

55 


(A. 3-8) 


To determine Hsj^(t), the matrix of first partials 
Hp^(t) must also be developed. The altitude measurement, 
taken from the current estimate of Shuttle position, is 
given by 

hp(t) = h^(t) = ijj (t) • ^(t) (A. 3-9) 

c 

where ift^(t) is a unit vector in the direction of computed 
position and ^(t) is the computed position vector. The 
perturbation ^n the expected altitude 6hp(t) can be determined 
by considering small perturbations about the nominal values 
in Eq. (A. 3-9). These perturbations are illustrated in Fig. 

A. 3-1. 


|•>‘f 
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Figure A. 3-1 Definition of Perturbations 

in Position Vector and Unit 
Relative Position Vector 

The difference between the true altitude and that predicted 
is 

6hp(t) = h(t) - h^(t) 

= iftCt) * R(t) - ij^ (t) • ^(t) (A. 3-10) 

c 

Using the definitions 

" 1r 

c 

6r(t) « R(t) - ^(t) 

in Eq. (A. 3-10) yields 

6hp(t) = ij^(t) • R(t) - [ij^(t) - filgCt)] • [fi(t) - 6r(t)] 

(A. 3-11) 

6ij^(t) and R(t) are perpendicular, so 6^p,(t) • R(t) = 0. To 
first order, an expansion of the last term in Eq. (A . 3-11 ) yields 
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&h(t) = i_(t) • 6r(t) 
p -R - 


The measur^ent matrix Hp^(t) is simply 


H„ (t) 


position errors 


I 0 

I 

3 


A. 4 COMPARISON OF MECHANIZATIONS 

The truth model measurement matrix using the drag 
measurement scheme is defined by 


Hg (t) - Hp (*) - 
q 9 <j 


which from Eqs. (A. 2-16) and (A. 2-6) is 


2i^p(g x) 

(t) = q(t) -TjS 

I “sc ''r 
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(A, 4-2) 
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If the truth model measur^ent matrix associated with the 
altitude measurement mechanization is similarly defined, i.e., 

Hg (t) = H (t) - H (t) (A. 4-3) 

h ^h h 

then Eqs. (A. 3-13) and (A. 3-8) yield 


Hg (t) 


V. 


^ j 

^ivR 

.fei 
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q { 
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(A. 4-4) 


Now, in the absence of measurement noise the 
following two equations hold: 

6q_(t) = H (t) 6x(t) (A. 4-5) 

6h (t) = H (t) 6x(t) (A. 4-6) 

o ©h 

/V 

where 6x(t) has been previously defined to be the estimate 
of the truth model state vector. In order to examine the 
effects of employing drag acceleration measurements versus 
altitude measurement the relationship between 5q (t) and 
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6hg(t) must be determined and then used to establish the 
proper relationship between Hg^Ct) and Hgj^(t). From Eq. 

(A. 3-3) the relationship between 5hg(t) and 6qg(t) is 
easily derived 

h 

6hg(t) * 6q (t) (A. 4-7) 

Substituting Eq. (A. 4-5) into Eq. (A. 4-7) gives 

h 

6hg(t) = - Hg (t) 6x(t) (A. 4-8) 

from which it is clear from Eq. (3.4-6) that if 

then the effects of employing drag acceleration measurements 
are no different than the effects of using altitude measure- 
ments. It is easily seen that Eos. (A. 4-2) and (A. 4-4) 
satisfy the relationship of Eq. (A. 4-9). Thus, to first 
order, the two mechanizations persented here yield the 
same effect in the performance evaluation of the drag update 
filter , 
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APPENDIX B 

TRUTH MODEL DETAILS FOR THE 
CAL/ ALIGN COVARIANCE ANALYSIS 


An error covariance simulation of the IMU calibration 
and alignment procedure consists basically of two operations: 
(1) propagation of IMU errors during a ground mode sequence, 
including least squares fit (LSF) estimation sequences and time 
delays for thennal stabixization , and (2) updates of IMU 
errors before and/or after an error propagation, including 
setting the initial platform misalignment after a platform 
slewing, torquing out platform tilts, storing LSF estimates 
for later use and updating IMU error compensation terms. 

This appendix contains some of the details required for de- 
fining the IMU error propagations and updates. Section B.l 
gives an explicit derivation of Coarse Alignment linearized 
error equations and Section B.2 contains an analysis of the 
platform misalignments at the end of a platform reposition- 
ing operation. 


B.l COARSE ALIGNMENT LINEAR ERROR ANALYSIS 

To implement the Coarse Alignment sequence of the 
covariance simulation program, it is necessary to express 
each measurement or computation performed during that sequence 
in terms of the IMU errors via a linear relationship. This 
linear relationship is then used to generate an update matrix 
as needed for Eqs. (3.2-5) and (3.2-6). 
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Coarse Alignment starts with a platform position 
given by caged gimbals. To keep the subsequent analysis simple 
assume, without loss of generality, that the initial platform 
orientation is such that the mechanical gimbal angles are 
zero. The first step then is to measure the gimbal angles: 


* 9. + PE 

m t 

'^m * ’^'t 

+ IRE (B.1-1) 


where 


4 , 9 '4/ . a are the measured outer roll, 

’ pitch, azimuth and inner roll 

gimbal angles, 

i = 6. = ® = 0 are the mechanical gimbal 

^ ^ ^ ^ angles, 

ORE, PE, AZE, IRE are errors incurred in mea- 
suring the corresponding gim- 
bal angles (resolver offsets, 
resolver harmonic errors, 
resolver randomness and quan- 
tization, etc.) 


The measured angles are then used to compute a trans- 
formation from navigation base to platform, Tp^g 
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(B.1-2) 
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Since the measured angles 4 ^, 9^^^, 0 ^^^, and are all small, 
Tp/B becomes 


where 



(B.1-3) 



(B.1-4) 


To determine the error in T. 


P/B 


first define the 


navigation base axes corresponding to a perfect IMU. Let 
these axes be represented by the subindex 0 (see Section A. 3 
of Ref. 1) and defined as follows: 


Xq = outer roll axis 

j’q * perpendicular to x^ and in the plane of the 
outer roll gimbal when the mechanical angle 
of this gimbal is zero, 

^0 = ’'o yo 


With these definitions TNBP can now be written as 


( 1 ) 

T = T T 
P/B P /0 0 /B 


(B.1-5) 







= 

I + 

a 

X 




(b-Dpj 






(B.1-6) 


(1) When the mechanical gimbal angles are zero then Tp.g is 
a function of the three small gimbal nonorthogonai iTies 
only . 
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where 


ipQ = Misalignment error of Tp. expressed in "0" 
coordinates 

a » Azimuth gimbal nonorthogonality 

b * Pitch gimbal nonorthogonality 

DP = Outer roll gimbal nonorthogonality calibration 
error 


Substituting Eq. (B.1-3) into Eq. (B.1-6) and solving for 
results in 




ORE + IRE 
PE - a 

,AZE - b + Dpy 


(B.1-7) 


The next step in Coarse Alignment is to estimate the 
transformation from reference to platform, , via a least 
squares fit to accelerometer data. Let be the actual 

transformation matrix generated, then 


"^P/LCA “ '^P/l'^L/LCA 


Tp/l(i * 




) 


(B.1-8) 


where are the misalignment errors of , expressed in 
reference (level) coordinates, and T^.. is the true transfor- 
mat ion at the beginning of Coarse Alignment. 

Tp/LCA 3.nd Tp^g are now combined to give the trans- 
formation from navigation base to reference, TNBR: 


* t^/lca^p/b 


(B.1-9) 
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Now, if the error in TNBR is represented by the misalignment 
vector 'i'jjgjj. then TNBR can also be expressed in terms of the 
transformation from outer roll coordinates to reference, 

TNBR - [l - (B.1-10) 

Combining Eqs. (B.1-5) to (B.1-10) to eliminate Tp^g, 
and TNBR one obtains 

" '^ca''] '^L/p'^P/0 "^0/L 

(B.1-11) 


or 


L 

^NBR 




L 

CA 


h/o*c 


(B.1-12) 


where is defined below. 


L 


(B.1-12) represents the error incurred 


In Eq 

in estimating the transformation ixom navigation base axes 
to reference axes, is a function of IMU errors as defined 

VJ 


by Eq. (B.1-7), and represents the error in estimating 

the transformation from reference to platform. The latter is 
related to the LSF filter output errors as follows 




L 

CA 



(B.1-13) 
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where 


6ft. is the LSF estimation error of the north 
acceleration 

602 is the LSF estimation error of the west 
acceleration 

69^ is the LSF estimation error of the north 
acceleration rate 

G is the local gravity acceleration 

is the north component of the earth's 
rotation 

At the end of Coarse Alignment the misalignments are 

stored in certain dummy states and saved for subsequent 
use. These misalignments would be used, for instance, when 
simulating a platform repositioning operation. The right 
hand side of Eq. (B.1-12) is a linear combination of several 
IMU errors, and this relationship will be used to generate 
the storage matrix (as required in Eqs. (3.2-5) and (3.2-6)). 

B.2 PLATFORM REPOSITIONING LINEAR ERROR ANALYSIS 

When completing a platform repositioning operation, 
the final platform position is such that the relationship 

Tp/B = TRD TNBR (B.2-1) 

is satisfied, were 

Tp,_ = transformation from navigation base to plat- 
' form computed from measured gimbal angles 
(see Eq . B.1-2) 

TRD = dea^red transformation from reference to 
platform, 

TNBR = transformation from navigation base to refer- 
ence determined during Coarse Alignment. 
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When Eq. (B.2-1) is satisfied, then the gimbal angles for a 
perfect gimbal train would be li* , 6. , and ip . The transforma- 
tion from outer roll axes (see definition in Section B.l) to 
platform axes is therefore 


T 

P/0 




(B.2-2) 


or 



1 

O 

O 

m 

'"t 

*sin 

'"t 

o’ 


cos 

®t 

0 

sin 




0 


0 


^P/0 " 

sin 

h 

cos 

♦t 

0 


0 


1 

0 



0 

cos 


-sin 



0 


0 


0 


-sin 

«t 

0 

cos 



0 

sin 


cos 

1 


(B.2-3) 


and the measured gimbal angles corresponding to the imperfect 
gimbal train are related to the ideal gimbal angles as follows: 




IRE 

ORE + DP tan 
PE-a 

lAZE - b + DP/cos e. 


(B.2-4) 


where 


IRE ,ORE,PE , AZE are the inner roll, outer roll, pitch, 
and azimuth resolver errors - including 
resolver offsets, resolver harmonics 
and resolver measurement errors and 
quant izat ion , 

a,b,DP are the azimuth, pitch, adjP outer 

roll gimbal nonorthogonality errors. 


The measured gimbal angles are then used to compute 
the transformation from navigation base to platform, resulting 
in 
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^P/B ' * “®’‘] \ [’ * e°x] 


or 


"^P/B “ "^P/O '^0/h J- + 6 )xj 




where 


= T. T. ,. a® 


•L/P *11', 


6^ * Tl/o®° 


IRE 

PE - a 
,AZE - b + DP/cos e. 


,0 


ORE + DP tan 6^ 
0 
0 


(B.2-5) 


(B.2-6) 


(B.2-7) 

(B.2-8) 


(B.2-9) 


(B.2-10) 


Now solving Eq . (B.2-1) for the desired platform 
orientation, TRD, and subsequently substituting Eqs. (B.2-5) 
and (B.1-10) into the result yields 

1 ■) = Tp^g INBr"^ (B.2-11) 

* "^P/L ' ^^]'^L/0 "^O/L '*NBR^] 

(B.2-12) 
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In Eq. (B.2-12),TRD and are the desired and the resulting 

transformation matrices from reference to platform, respectively. 
Therefore, the initial platform misalignment resulting from a 
cluster repositioning is 


♦a - (B.2-13) 

or, after substituting Eqs. (B.1-12), (B.1-7), (B.2-9), and 
(B.2-10) 


'^'CR " 


IRE 

PE - a 
AZE - b + DP /cos 6, 


- T 


L/Ol 


(• 

y AZE - b + DP 


RE - DP tan 6 \ 

».\ T 

+ ^ 


PE - a 


(B.2-14) 


Note that when ^ ^ then the error in repositioning 

the platform is exactly equal to the error in estimating the 
transformation from reference to platform during Coarse Align- 
ment, which is what one would expect. 


In the covariance simulation program a platform 
repositioning can therefore be simulated by an instantaneous 
change in platform orientation followed by an assignment of 
initial platform misalignments according to Eq. (B.2-14). 


/5A 
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APPENDIX C 

ERROR CONTRIBUTION TI!iE HISTORIES: 
DEORBIT TO 20.000 FT 


This appendix presents the time histories associated 
with the error budget results summarized in Section 4,4.1. 

The data is presented in tabular form with each table indi- 
cating the rms errors in position and velocity, and the 
rms value of the platform alignment estimates, in both the 
R (runway) and V (relative velocity) coordinate frames which 
result from a specific error source or group of errors. 


The time points in each table begin at deorbit 
burn and are printed at appropriate time intervals down 
to MLS acquisition. Time points are also included just 
before and after each external navigation aid is activated. 
Thus there are two rows of data for each of the following 
times (see Table 4,1-1): 


"^DRAG 

"^TACAN 


= 1375.19 sec 
= 2239.19 sec 
2496.47 sec 


The final row in each table presents the rms errors just 
prior to the first MLS measurements at an altitude of 
20,000 ft: 


Tjjl 3 “ 2723.03 sec 


In addition, the rms value of each of states 10 through 
12 is given at 20,000 ft. The magnitudes and mathematical 
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description of the error sources are given in Section 4.3. 
Units are in feet, fps, and radians. 
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APPENDIX D 

ERROR CONTRIBUTION TIME HISTORIES: 
20.000 FT TO TOUCHDOWN 


This appendix presents the time histories associ- 
ated with the error budget results summarized in Section 
4.4.2. The data is presented in tabular form with each 
table Indicating the rms errors in position and velocity, 
and the rms value of the platform alignment estimates, in 
both the R (jmnway) and V (relative velocity) coordinate 
frames, which result from a specific error source or group 
of errors . 

The time points in each table begin at MLS acquisi- 
tion and are printed at appropriate time intervals there- 
after Time points are also included just before and 
after processing of the first radar altimeter measurement. 
Thus there are two rows of data for each of the following 
times ( e Table 4.1-1): 

Tj^Ls “ 2723.03 sec 

T^ = 2809.43 sec 

A single row of each table corresponds to the initiation 
of the landing phase: 

Tg^ = 2749.91 sec 

The magnitudes and mathematical description of the error 
sources are given in Section 4.3. Units are in feet, fps, 
and radians. 
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